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SECTION  1 
INTRODUCTION 


In  this  report  we  shall  use  simple  models  to  study  the  phenome¬ 
non  which  has  come  to  be  known  as  nuclear  winter1.  The  models  are  time- 
independent,  one-dimensional  and  allow  solution  in  terms  of  closed  form, 
analytic  expressions.  In  view  of  the  fact  that  time-dependent,  multi¬ 
dimensional  calculations  of  the  phenomenon  have  been  presented1,  it  is 
perhaps  appropriate  to  motivate  the  use  of  the  much  simpler  models.  The 
models  have  several  uses:  1)  they  allow  one  to  develop  an  understanding 
of  the  phenomenon  at  a  more  fundamental  level  than  may  be  possible  with 
the  use  of  large,  complex  computer  simulations;  2)  they  allow  one  to 
examine  the  parameter  space  of  the  problen  or  to  examine  the  importance  of 
including  a  particular  type  of  phenomenon  in  the  calculations.  Usually  a 
simple  inspection  of  the  formulas  will  reveal  the  sensitivity  to  a  param¬ 
eter  value  or  type  of  effect,  at  most  a  simple  calculation  will  suffice. 
In  this  way  they  can  point  to  interesting  regions  of  the  parameter  space 
where  more  detailed  calculations  should  be  done  or  to  the  most  important 
additions  to  make  in  improving  the  more  complicated  calculations;  3)  the 
models  and  the  physical  intuition  they  allow  one  to  develop  are  useful  in 
developing  "sub-grid"  models:  formulas  which  specify  in  large  scale  simu¬ 
lation  codes  the  results  of  processes  at  scales  smaller  than  the  code  can 
resolve. 

The  main  use  we  shall  nake  of  the  models  in  this  report  -  the 
purpose  for  which  we  developed  them  -  is  to  study  the  effects  of  infrared 
radiation,  I R ,  on  the  nuclear  winter  phenomenon.  In  particular  we  shall 
consider  whether  or  not  the  presence  of  such  IR  radiators  as  hydrocarbons. 
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N02,  C02,  water  vapor  or  condensed  water  (clouds)  can  be  of  importance  to 
the  phenomenon  and  to  consider  the  sensitivity  of  the  phenomenon  to  rea¬ 
sonable  variation  in  the  IR  properties  of  soot  itself.  We  find  that 
hydrocarbons,  NO2  and  CO2  are  not  of  importance;  this  result  is  in  agree¬ 
ment  with  the  conclusions  of  other  workers.  Water  vapor  is  unlikely  to  be 
of  more  than  marginal  importance  to  the  temperature  structure  if  one 
ignores  convective  instabilities;  it  may  be  of  importance  in  determining 
the  stability  characteristics  at  the  top  of  the  soot  layer.  Variation  in 
the  IR  properties  of  the  soot  itself  across  the  range  of  reasonable  expec¬ 
tation  makes  a  significant  difference  in  the  results.  Condensed  water,  if 
it  is  present,  can  play  a  significant,  indeed  a  dominant,  role  in  the  phe¬ 
nomenon. 

A  water  cloud  mixed  with  the  soot  will  alter  the  entire  charac¬ 
ter  of  the  temperature  profile,  changing  it  from  one  in  which  the  tempera¬ 
ture  is  hot  at  the  top  of  the  layer  and  cold  at  the  earth's  surface  -  as 
shown  in  previous  calculations  of  the  nuclear  winter  phenomenon  and  as 
shown  in  our  calculations  in  the  absence  of  condensed  water  -  to  one  in 
which  the  temperature  is  cold  at  the  top  of  the  cloud-soot  layer  and  hot 
at  the  earth's  surface.  Furthermore  the  clouds  can  be  created  by  the 
amounts  of  water  which  are  likely  to  be  present.  The  question  of  whether 
or  not  there  will  be  clouds  is  complicated.  There  exists  a  strong  hyster¬ 
esis  effect:  if  the  cloud  is  there  it  keeps  the  top  of  the  layer  cold 
which  keeps  the  cloud  there;  if  the  cloud  is  not  there  the  top  of  the 
layer  is  hot  which  keeps  a  cloud  from  forming.  Most  fire  plumes  are 
initially  capped  with  condensation  clouds2. 

If  the  clouds  persist  for  long  periods  of  time  they  will  almost 
certainly  have  a  profound  effect  on  the  nuclear  winter  phenomenon.  They 
will  also  affect  the  dynamics  of  the  phenomenon  at  early  times  when  the 
various  individual  plumes,  some  of  which  will  probably  have  condensation 
clouds,  some  of  which  may  not  have  condensation  clouds,  combine  to  form 
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the  hemisphere-wide  soot  layer  of  nuclear  winter.  Even  if  the  condensa¬ 
tion  clouds  disappear  before  the  late  times  their  effects  on  the  dynamics 
of  the  combining  phase  may  have  implications  for  removal  rates  of  the  soot 
and  for  the  later  location  of  the  material.  In  turn,  it  is  the  outcome  of 
these  mixing  processes,  which  the  clouds  will  influence,  which  will  deter¬ 
mine  whether  or  not  there  will  be  clouds  at  late  times.  Until  careful 
calculations  which  include  the  effects  of  the  clouds  are  made  we  do  not 
believe  that  any  calculation  of  the  nuclear  winter  phenomenon  can  be  com¬ 
pletely  credible:  the  presence  of  condensation  clouds,  which  we  have  not 
been  able  to  rule  out,  would  alter  the  entire  character  of  the  event. 

We  also  point  out  that  the  nuclear  winter  phenomenon  is  always 
wanning  over  snowpack.  This  fact  may  influence  global  circulation  pat¬ 
terns  or  may  be  important  in  other  regards.  The  analysis  of  whether  the 
warming  effect  over  snowpack  is  of  importance  to  the  total  phenomenon  is 
beyond  the  scope  of  the  present  report. 

The  report  is  organized  as  follows:  In  Section  2  we  document 
the  data  we  found  necessary  to  accumulate  in  order  to  specify  the  param¬ 
eters  or  ranges  for  parameters  in  the  model  calculations.  Much  of  this 
data  was  gleaned  from  other  sources;  some  represent  calculations  we  made 
ourselves.  In  Section  3  we  present  the  model  calculations  -  the  main 
results  of  the  study.  The  models  are  described  in  order  of  increasing 
complexity,  first  ignoring  convective  instabilities  then  including  their 
effects  in  an  approximate  way.  Concluding  renarks  are  given  in  Section  4 
and  in  Appendix  E  we  make  some  remarks  on  calculating  the  develoDment  of 
the  unstable  interface  at  the  top  of  the  soot  layer  when  no  condensation 
cloud  is  present. 
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SECTION  2 

GASEOUS  AND  PARTICULATE  MATTER  INJECTED  INTO  ATMOSPHERE  BY  MASSIVE 
BURNING,  AND  THEIR  OPTICAL  PROPERTIES 


PRELIMINARY 

One  purpose  of  the  present  study  is  to  investigate  the  possible 
amelioration  of  the  predicted  "nuclear  winter"  effects  at  the  earth's  sur¬ 
face  by  polyatomic  molecules,  injected  into  the  atmosphere  following 
massive  burning,  that  might  inhibit  the  escape  to  space  of  thermal  radia¬ 
tion.  For  this  purpose,  we  first  performed  a  literature  survey  to  deter¬ 
mine  the  principal  gaseous  species  produced  by  the  burning,  and  to  obtain 
estimates  of  the  guantity  of  each  injected.  We  then  sought  to  determine, 
using  band  models  and  parameters  from  the  literature  where  possible,  the 
broad-band  infrared  (IR)  emissivity  and  transmissivity  for  each  of  the 
species. 


As  a  further  aid  in  performing  the  work  of  this  contract,  we 
undertook  calculations  to  estimate  the  broad-band  IR  emissivity,  transmis¬ 
sivity,  and  reflectivity  for  particulate  layers.  These  include  water 
clouds,  dust  layers,  soot  (carbon)  layers,  and  layers  containing  a  mixture 
of  water  droplets  and  soot.  We  also  performed  calculations  to  determine 
the  corresponding  optical  properties  in  the  visible  (0.F5  um)  for  the  same 
particulate  layers. 

This  section  describes  the  methodology  employed  in  the  foregoinq 
studies  and  presents  the  results  in  graphical  and  tabular  form.  The 
results  of  the  literature  survey  of  injected  species  and  their  emission 
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This  is  followed  by  the  optical  properties  of 
finally,  by  the  optical  properties  of  the 


MOLECULAR  SPECIES  INJECTED 


A  survey  of  the  literature,  including  References  1  through  5, 
shows  that  the  principal  species  injected  into  the  atmosphere,  as  a  conse¬ 
quence  of  forest  and  urban  fires,  are  C02 ,  CO,  H20,  N0x(N02+N0),  and 
hydrocarbons.  A  key  item  for  each  of  the  species  is  the  "emission 
factor,"  which  is  the  number  of  grams  of  the  substance  emitted  per  gram  of 
fuel  burned.  These  numbers  are  by  no  means  unique  for  each  of  the 
species,  as  different  measurements  have  revealed  a  range  of  values. 

Table  2-1  shows  values  for  the  emission  factors,  the  associated  refer¬ 
ences,  and  values  adopted  in  this  report. 


There  are  many  hydrocarbons,  but  we  have  assumed  that  they  are 
dominated  by  six  species  and  proportioned  among  them  as  shown  in 
Table  2-2.  The  corresponding  emission  factors  adopted,  based  on  the  total 
hydrocarbon  emission  factor  and  the  relative  abundances,  are  also  shown  in 
Table  2-2. 


In  Table  2-3  we  compare,  data  permitting,  the  species  content  in 
the  ambient  atmosphere  (30°N  -  70°N;  0.33  tt  steradians)  with  the  addi¬ 
tional  amount  added  from  the  burning  of  4500  Tg  fuel  (1  Tq  =  IQ17  q) . 

This  is  the  "baseline  case"  considered  in  Reference  2-1  and  corresponds  to 
3750  Tq  consumed  in  urban  fires  and  750  Tq  in  forest  fires.  The  amounts 
added  from  the  burning  are  based  on  the  emission  factors  given  in 
Tables  2-1  and  9-2.  In  both  cases  the  amounts  are  given  in  terms  of  the 
total  number  of  molecules  and  the  vertical  colunn  number  in  units  of  atm 
cm  at  NTP.  For  the  ambient,  atmosohore  we  have  used  the  followinu:  C/V,  and 
CO  mixing  ratios  of  3  *  Io“4  and  i.u  <  l’)-7,  respective! y,  witn  tne 
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Table  2-1.  Emission  factors  for  products  of  combustion 


Species 

Emission  Factor 

Reference 

Emission  Factor 
Adopted 

C02 

1  —  1.75 

2-2 

1.75 

1.9  (Ponderosa  Pine) 

2-4 

CO 

.05 

2-1 

0.18 

.01  _  .25 

2-2 

.01  (Ponderosa  Pine) 

2-4 

H2n 

.25  —  .75 

2-2 

0.75 

NOx 

1.5  *  10-3 

2-1 

5  x  io-3 

5  x  10“4  _  4.5  x  io-3 

2-2 

4  x  m-4  _  3.3  *  io-2 

2-3 

hydrocarbons 

,.02 

2-1 

(total ) 

5  x  io-3  _  2  x  10-2 

2-2 

4.5  x  io-3  __  6.5  x  io-3 

2-4 

Table  2-2.  Assumed  distribution  of  hydrocarbons  (Reference  2-4). 


Species 

Relative  Abundance  [%) 

Emission  Factor 
Adopted 

CH4  (methane) 

30.8 

3.1  *  IO"3 

C2H4  (ethylene) 

25.8 

2.6  *  IO'3 

C2H2  (acetylene) 

22.3 

2.2  x  io-3 

C3H6  (oropylene) 

7.1 

1 

c 
*— 1 

X 

C3H8  (propane) 

7.0 

7.0  *  10-4 

C2H6  (ethane) 

7.0 

7.0  *  10-4 
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atmospheric  content  and  additions  thereto  in  "baseline  case 


total  concentration  from  Reference  2-6;  the  mid-latitude  mean  values  for 
H2n  from  Reference  2-6;  N0X  and  CH4  values  from  References  2-7  and  2-8. 

The  ambient  atmospheric  content  of  the  hydrocarbons,  other  than  methane, 
is  not  readily  available  to  us  and  so  is  not  included  in  Table  2-3. 

BROAD-BAND  IR  PROPERTIES  OF  MOLECULAR  SPECIES 

Since  thermal  radiation  from  the  earth  approximates  blackbody 
radiation,  we  are  interested  in  determining  the  emission  and  absorption 
properties  of  the  injected  products  of  combustion  averaged  over  a  Planck 
spectrum.  These  averaqe  values  are  called  "broad-band"  IR  properties  of 
the  molecules. 

The  effective  radiatinq  temperature  of  the  earth  lies  between 
about  200  and  300  °K.  At  300  °K  the  blackbody  emission  peaks  at  10  urn;  at 
200  °K  it  peaks  at  about  IS  urn.  So  we  are  interested  mainly  in  molecules 
that  absorb  in  the  lonq-wavelenqth  infrared  (LWIR)  beyond  about  6  urn.  All 
of  the  molecules  listed  in  Table  2-3  satisfy  this  requirement  except  for 
CO  whose  fundamental  band  lies  at  about  4.7  urn.  Thus,  althouqh  Table  2-3 
indicates  that  massive  hurninq  may  inject  more  CO  into  the  atmosphere  than 
is  there  under  ambient  conditions,  its  presence  should  not  alter  siqnifi- 
cantlv  the  transport,  of  thermal  radiation  through  it,  and  so  we  delete  it 
from  further  consideration.  We  shall  also  iqnore  C02.  This  is  because, 
from  Table  2-3,  the  amount  injected  bv  burninq  is  only  about  1  percent  of 
the  ambient  content.  With  the  exception  of  propylene  (C3H&),  for  which  we 
could  find  no  hand-model  data,  estimates  of  the  broad-band  IR  properties 
of  the  remaininq  molecules  in  Table  2-3  are  given  below.  For  NOx  we 
will  include  only  N02  because  very  little  of  the  fundamental  band  of  ND  at 
5.3  urn  extends  into  the  LWIR  region  at  ambient  temperatures. 
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The  transmission  of  NO2  in  the  v2  and  v3  fundamental  bands  and 
the  vi  +  u3  combination  band,  spanning  the  wavenumber  interval  655  to  2950 
cnr1,  has  been  modeled  by  Pierluissi  and  Tomiyama  (Reference  2-9)  and  can 
be  written  in  the  form 


where 


,  3[+a?X . 

T  =  exp(-10 


X  =  c‘  +  loqinW 


(2-1) 


(2-2) 


n  j  m 

H  -  (p~)  (-2)  u  . 


(2-3) 


Here,  P  is  the  total  air  pressure,  T  is  the  temperature  (°K),  P0  and 
T o  (273  °K)  are  the  corresponding  STP  values,  U  is  the  absorber  thick¬ 
ness  (atm  cm),  ai ,  a2,  n,  m  are  constants  given  in  Reference  2-1*,  and  c' 
is  tabulated  in  Reference  2-1  as  a  function  of  wavenumber  v. 

The  spectral  emissivity  at  wavenumber  v  is 


e(v)  =  1  .  T(v) 


(2-4) 


and  the  broad-band  emissivity,  e,  is 


/  e(v)R(v)dv 


—  _  0_ 


(2-5) 


/  R(v)dv 


*  The  value  of  a^  is  erroneously  given  in  Reference  2-1  as  -0.25653.  The 
correct  value  is  1.19711. 


t  V 

,  -\  -v *,-*v  *,*  v  v*  v  >  v  v\  WvSV.*\- 


Here,  B(v)  is  the  Planck  function 

B(v)  =  2hcv3(e1-435  v/T-l)-1  .  (2-6) 

We  have  evaluated  c  numerically  from  Equation  2-5  using  Equa¬ 
tions  2-1  through  2-4,  in  which,  for  the  nuclear  winter  application,  we 
assumed  an  average  N02  layer  height  of  5  km  with  P  =  540.5  mb  and 
T  =  256  °K.  The  calculations  were  performed  for  N02  layer  thickness  (U) 
ranging  from  10"4  to  102  atm  cm.  The  results  are  shown  in  Figure  2-1. 

For  the  maximum  value  of  U  for  N02  in  Table  2-3,  that  incor¬ 
porates  contributions  from  fires  and  nuclear  fireballs,  the  emissivity 
(Figure  2-1)  is  seen  to  be  less  than  about  6  x  10"3.  This  corresponds  to 
a  transmission  of  thermal  radiation  from  the  earth's  surface  of  at  least 
99.4  percent. 

H20 

The  broad-band  IR  emissivity  for  water  vapor,  as  reported  by 
Staley  and  Jurica  (Reference  2-10),  is  shown  in  Figure  2-2  for  two  differ¬ 
ent  temperatures.  In  this  figure  the  path  length  is  in  units  of  gm  cm-2 
of  water  vapor  (1  gm  cm-2  =  1.2438  *  103  atm  cm).  For  an  atmospheric 
water-vapor  content  of  1.8  *  103  atm  cm  (Table  2-3),  corresponding  to  1.45 
gm  cm-2,  the  emissivity  is  seen  to  be  about  0.63,  corresponding  to  a 
transmission  of  thermal  radiation  of  37  percent.  These  values  are  for  the 
ambient  atmosphere.  However,  the  additional  amount  of  water  vapor  in¬ 
jected  by  the  fires  (3.6  atm  cm)  contributes  only  insignificant  increnents 
to  the  ambient  emissivity  and  transmissivity. 


WATER  VAPOR  PATH  LENGTH  u(?cm'£) 


Fioure  ?-?.  Rroad-band  IR  emissivity  for  water  vapor 
(from  Reference  2-10). 


Hydrocarbons 

Methane  (CH4).  For  methane  we  use  the  band  model  developed  bv 
Jarem,  et  al .  (Reference  2-11)  which  is  formally  similar  to  that  described 
above  for  N02 .  The  transmission  at  wavenumber  v  is  qiven  by 

T(v)  =  exp[-  ioa(c,(w)+1°9ioH)]  (2_7) 

where  W  is  determined  from  Equation  2-3.  The  values  of  the  parameters 
a,  n,  m  are  qiven  in  Reference  2-11  as  are  also  the  parameters  c'(v)  for 
v  ranqinq  from  10R5  to  3215  cm'1.  This  interval  spans  the  reqion  of  the 
two  IP-active  fundamentals  v3(1306  cm-1)  and  v4(3020  cm-1).  Aqain,  the 

spectral  emissivity  is  qiven  by  Equation  2-4  and  the  broad-band  emissivity 
by  Equation  2-5. 


The  broad-band  emissivity,  at  a  temperature  of  256  °K  and  total 
pressure  of  0.54  atm,  is  shown  in  Figure  2-3  (along  with  results  for  the 
other  hydrocarbons  described  below)  as  a  function  of  path  lenqth  (atm  cm) 
through  the  methane  layer.  Arrows  near  the  bottom  of  Figure  2-3  show  the 
values  of  U  from  Table  2-3  corresponding  to  the  fire  contribution  in  the 
baseline  case.  For  methane,  the  ambient  value  for  U  is  also  shown.  It 
corresponds  to  an  emissivity  of  about  10-2,  which  is  eguivalent  to  a 
transmissivity  of  99  percent.  The  incremental  amounts  from  the  products 
of  combustion  are  only  about  0.1  percent. 

Ethane  (C2H6).  For  ethane,  we  have  been  able  to  find  absorp¬ 
tion-coefficient  data  for  only  the  v9  band  centered  at  821.52  cm*1 
(12.2  urn).  Although  this  is  believed  to  be  the  major  band  of  interest  to 
us,  the  results  obtained  below  should  be  considered  as  lower  limits  to  the 
broad-band  emissivity. 


The  spectral  absorption-coefficient  data  for  the  vg  band, 
obtained  by  Varanasi,  et  al .  (Reference  2-12),  can  be  represented  as 

k_  =  s  exp[-a(v-v0)2]  (atm*1  cm*1)  ,  (2-8) 

V  *  n 

where 

oo 

S  =  J  k.  d(v-vQ)  =  31  (ML)  (atm*1  cnr1) 

oo  ^ 

a  =  5.79  x  10-4  cm2 
and 

'•>  =  825  cm-1 

o 
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Assuminq  the  statistical  (Lorentz  line)  band  model,  the  trans¬ 
mission  is  given  by  (Reference  2-13) 


T(  v)  =  exp^  -k  -U/^l+k  ll/(4a  )]  (2-R) 

v  VC 

where,  aqain,  H  (atm  cm)  is  the  ahsorber  amount  and 

ac  =  Yc/d  .  (2-10) 

Here,  Yc  and  d  are,  respectively,  the  mean  (collision  broadened)  line 
width  and  mean  line  soacinq. 

For  ethane,  we  take  (Reference  2-12) 


1/2 

Yc  -  0.102  (M)  Pgtm  ;  d  -  2.0  cm"1  (2-11) 

where  Daim  is  the  total  pressure  in  atmospheres.  Althouqh  ac  is 
qenerallv  a  function  of  v,  we  assume  it  to  be  a  constant.  With  the  above 
parameter  values,  alonq  with  values  for  T  and  Patm  °f  25^  °K  anc^ 

0.54  atm,  respect i vely,  Fauations  2-4,  2-5,  and  2-Q  lead  to  the  broad-band 
emissivitv  for  ethane  shown  in  Fiqure  2-3.  For  the  amount  of  ethane  in¬ 
jected  into  the  atmosphere  by  massive  burning  (2.1  *  10-3  atm  cm),  the 
emissivitv  is  seen  to  be  very  small  (6  *  10~5),  corresponding  to  almost 
complete  transparency  of  the  thermal  radiation.  It  is  unlikely  that  this 
result  would  he  siqnif icantl v  different  had  the  other  bands  of  ethane  been 
included. 


Propane  (r3Hfl).  For  propane  we  use  the  band-model  parameters  of 
diver,  et  al .  (Reference  2-14).  The  snectral  transmission  can  be  repre¬ 
sented  hv  the  expression 


where 


T(v)  =  exD[-kvU/*'l+kyU/(Tfac)  ] 


(2-12) 


ac  = 


atm 


(2-13) 


From  plots  of  a°(atm"1)  and  kv(atm_1  cm"1)  as  functions  of  v  given  in 
Reference  2-14,  we  have  computed  the  spectral  transmissivity  and  emis¬ 
sivity  from  Fquations  2-12  and  2-4,  respecti vely,  and  the  broad-band  IR 
emissivitv  from  Fquation  2-5,  as  functions  of  U.  The  results  are  shown  in 
Fiqure  2-3. 


From  the  amount  of  propane  injected  into  the  atmosphere  in  the 
baseline  case  (1.4  *  in-3  atm  cm),  we  see  that  the  emissivity  is  nearly 
the  same  as  that  for  ethane  (6  *  10" 5 )  and  that  any  decrease  in  the  trans¬ 
mission  of  thermal  radiation  due  to  injected  propane  is  quite  neqligible. 

Acetylene  (C2H2)-  Integrated  band  intensities,  S,  at  300  °K 
have  been  measured  for  several  bands  of  acetylene  by  Varanasi  and  Ranqaru 
(Reference  2-15).  These  are  shown  in  Table  2-4.  Of  the  four  shown,  two 
are  fundamentals  and  two  are  intercombination  bands.  The  v3  fundamental 
and  the  v2  +  v4  +  v5  intercombination  band  severely  overlap  and  are 
treated  as  a  single  band.  We  call  this  "band  3."  The  vs  fundamental  and 
the  v4  +  v5  intercombination  band  will  be  referred  to  as  bands  "1"  and 
"2,"  respectively.  Only  for  band  1  are  spectral  absorption  coefficients 
given  in  Reference  2-15.  To  obtain  the  broad-band  emissivity,  including 
contributions  from  bands  2  and  3,  we  use  an  approximation  developed  in 
Appendix  A. 


Table  2-4.  Rand  intensities  for  acetylene  at  300  *K  (Reference  2-15). 


Band 

Band  Center  (cm-1) 

Wavelenqth  (pm) 

S  (atm-1 cm"2) 

(1)  v5 

729 

13.7 

729 

(2)  +  v5 

1328 

7.53 

87 

(3)1  "3  1 

'jv2  +  v4  +  v5j 

(32951 
(3282(  3?qn 

3.04 

294 

Consider  first  "band  1."  As  for  the  case  of  ethane,  we  write 
the  soectral  transmissivity  in  the  form 


T(v)  =  exo[-k~H/*'l+k~U/(4ac)  ] 


(2-14) 


where 


a  =  Y  /d  =  a  P  . 
c  c  c  atm 


(2-15) 


We  have  been  unable  to  find  values  of  a°  for  acetylene.  However,  for 

Drooane  near  750  cm-1,  a°  “  1  atm-1  (Reference  2-14),  and  we  arbitrarily 

adopt  this  value  here.  Values  for  k-  (atm-1  cm-1)  from  550  to  810  cm-1 

v 

are  available  from  Reference  2-15.  These,  together  with  equations  2-14, 
2-15,  2-4,  and  2-5,  were  used  to  compute  the  broad-band  emissivity  for  the 
v5  band. 

In  Appendix  A  we  show  that  if  an  absorbinq  layer  is  optically 

thin  to  a  molecular  band  (k-l)  «  1),  then  the  broad-band  emissivity  for 

v 

that  hand  is  given  bv 
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•A  v*i 


(2-16) 


e  -  1.53  B*  (  v0)SU/T1* 

where 

B'(v0)  =  v03(e1-435  V0'/T-1  )-1  ,  (2-17) 

v0  is  the  hand-center  wavenumber,  S  is  the  integrated  band  intensity 
(atm~  cm*  ),  and  H(atm  cm)  is  the  layer  thickness.  It  turns  out  that 
under  nuclear-winter  conditions,  values  of  l)  for  acetylene  are  suffi¬ 
ciently  small  (Table  2-3)  that  the  layer  is  optically  very  thin  and  so 
the  condition  for  applicability  of  Eauation  2-16  is  easily  satisfied. 
Consequently,  for  bands  2  and  3,  at  a  temperature  of  256  ‘K,  and  with 
values  for  v0  and  S  from  Table  2-4,  Equation  2-16  leads  to  the  result 

^2  +  73  =  1.8  *  lO”2  IJ  .  (2-18) 

Equation  2-12  for  the  spectral  transmissivity  is  called  the 

"curve  of  qrowth."  For  small  values  of  k-tl,  the  exDonent  is  linear  in 

v 

k-(J.  Rut  for  larqe  enouqh  k -II  the  exponent  varies  as  the  square  root  of 
v  v 

k-tl.  In  the  linear  reqion  it  turns  out,  from  our  detailed  calculations, 
or  from  Fquation  2-16,  that  for  "band  1"  we  can  write 

e\(k-U  «  1)  =  7.3  x  in-1  u  .  (2-10) 

v 

TomDarinq  Fquations  2-18  and  2-10  we  see  that  the  contribution  to  the 
broad-band  emissivity  from  "band  1"  is  about  42  times  qreater  than  that 
from  the  sum  of  bands  2  and  2.  Consequently,  bands  2  and  3  can  be  iqnored 
with  very  little  error.  This  permits  us  to  plot  the  broad-band  emissivity 
for  acet.vlene  as  a  function  of  M,  even  to  larqe  values  of  II. 


if 


m 
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The  results  are  shown  in  Fiqure  2-3.  For  the  value  of  U  in  the 
baseline  case  (7.4  *  10"  ),  we  see  that  the  emissivity  is  about  5  *  10"  , 
which  corresDonds  to  a  decrease  in  the  transmission  of  only  0.5  percent. 

Ethylene  (C 2^4).  We  have  found  no  band  model  nor  spectral 
absorption  coefficients  in  the  literature  for  ethylene.  However, 

Nakanaqa,  et  al.  (Reference  2-16)  have  measured  inteqrated  band  inten¬ 
sities,  5,  for  the  v7,  v9,  vio»  vn,  and  v12  bands.  The  results  are  shown 
in  Table  2-6.  For  sufficiently  small  values  of  1),  we  can  determine  the 
broad-band  emissivity  from  the  equation  (Appendix  A) 

7  -  1.53  l  R’ (v  )S. U/T4  (2-20) 

i 

Table  2-5.  Rand  intensities  for  ethylene  (Reference  2-16). 


Rand 

v^cm"1) 

Wavelenqth  (urn) 

S^atm"1  cm"2) 

v; 

94Q 

10.54 

376.9 

v9 

3105 

3.22 

116.1 

v10 

826 

12.11 

0.134 

vn 

2989 

3.35 

63.9 

v12 

1443 

6.93 

4614 

where  the  sum  in  Fouation  2-20  is  taken  over  all  molecular  bands  and 
R‘ (vf)  is  qiven  hv  Equation  2-17.  When  the  parameters  from  Table  2-5 
are  used  in  Equation  2-20,  we  obtain,  for  T  =  256  °K, 


—  optical  1  v  thin  s  n  ?t.  n 

ethylene  '  '  atm  cm 


(2-21) 


which  is  plotted  in  Fiqure  2-3. 


Equation  2-21  should  be  valid  for  the  nuclear-winter  application 
in  which,  for  the  baseline  case,  U  -  8.2  *  10'3  and  so  e  ■  2  *  10"3. 

This  corresponds  to  a  decrease  in  the  transmission  of  thermal  radiation 
due  to  the  ethylene  layer  of  only  0.2  percent. 

If  we  add  the  effects  of  all  the  hydrocarbons  resulting  from 
massive  burninq  (baseline  case),  we  find  a  total  broad-band  emissivity  of 
about  8.1  *  10-  #  However,  due  to  band  overlap  among  the  different 
soecies,  the  resultant  emissivity  will  be  somewhat  less  than  this  (see 
Appendix  R).  We  therefore  expect  that  the  injected  hydrocarbons  will 
decrease  the  transmission  of  thermal  radiation  through  the  layer  by  less 
than  1  percent. 

BROAD-BAND  IR  PROPERTIES  OF  PARTICULATE  LAYFRS 

In  this  subsection  we  describe  calculations  performed  to  esti¬ 
mate  the  broad-band  IP  properties  of  particulate  layers  of  interest  in  the 
nuclear-winter  problem.  Of  special  relevance  are  the  emissivity,  trans¬ 
missivity,  and  reflectivity  of  water  clouds,  dust  layers,  soot  (carbon) 
layers,  and  layers  containing  a  mixture  of  soot  particles  and  water 
droplets.  We  first  describe  the  method  employed  and  then  present  the 
results  for  specific  particle  types. 

Method 

We  are  interested  in  determining  the  broad-band  IR  emissivity, 

7,  transmissivity,  T,  and  reflectivity,  p,  for  particle  layers  of  thick¬ 
ness  U  (gm  cm-2).  These  quantities  are  related  by  the  expression 

7  +  7  +  p”  =  1  .  (?-??) 

Here,  e  is  defined  bv  the  relation  2-8,  and  T  and  R  are  defined  by 
corresponding  expressions  in  which  the  spectral  emissivity,  e(v),  is 
replaced  by  the  spectral  coefficients  t(v)  and  R(v),  respectively. 


Initially,  we  attempted  to  calculate  the  spectral  coefficients 
e(v),  T(v),  and  R(v)  usinq  a  Mie-scatterinq  code,  and  a  radiation  trans¬ 
port  alqorithm  similar  to  one  in  the  ROSCOE  and  NORSE  codes.  However,  we 
ran  into  numerical  difficulties  with  the  method,  especially  when  consider- 
inq  mixtures  of  particles  where,  for  certain  values  of  l)  and  density 
ratios  of  particle  types,  the  sum  2 -22  exceeded  unity.  Subsequently,  we 
adopted  an  approximate  procedure,  summarized  by  Bartky  and  Bauer  (Refer¬ 
ence  2-17),  for  calculatinq  the  spectral  coefficients  in  terms  of  the  Mie- 
scatterinq  cross  sections  that  we  obtained  from  our  Mie  code. 

The  appropriate  formulas  for  the  spectral  coefficients,  using 
the  parallel  beam  or  one-dimensional  approximation,  are  as  follows: 


e  =  2<[ekT(k+<)  +  e'kT(k-<)  -  2k] /D 
R  =  (k2  -  <2)(ekT  -  e’kT)/n 


where 


T  =  4<k/n 


n  =  (k  +  <)2  ekT  -  (k  -  <)2  e'kT 
k  =  {<[2(1  -  <16+  <] } */2 
*  =  +  a  J  =  d,/(o,  +  a  ) 

d  d  b  d  d  b 


a  +  B  =  1 


(2-23) 

(2-24) 


(2-25) 


(2-26) 

(2-27) 

(2-28) 


(2-29) 


B  =  0.5  J  P(9)sin0d0 
*/2 


0.5  /  P(u)du 
-1 


(2-30) 


Here,  o  and  a  a re,  respectively,  the  absorption  and  scattering  cross 
^  2  ^ 

sections  (cm  ),  a  and  a  are  the  related  mass  absorption  and  mass- 
a  s  2  i 

scattering  cross  sections  (cm  qm-  ),  a  and  e  are  the  fractions  of  the 
radiation  scattered  forward  and  backward,  respectively,  at  a  single  scat- 
terinq,  and  P(9)  is  the  scattering  pattern  normalized  such  that 

1 

n.R  /  P( u)du  =  1  ,  (2-31) 

-1 

where  u  =  cos9  and  0  is  the  scattering  angle.  The  quantity  t  is  the 
optical  depth  which,  for  a  distance  z  into  a  medium  of  particle  number 
densitv  N,  is  given  by 

z 

*  =  /  (o  +  ojNdz1  .  (2-32) 

a  S 

0 

Alternatively, 

t  =  (q  +  o  )••  =  q  II  ,  (2-33) 

a  b  tr 

where  a  ,  the  sum  of  q  and  q  ,  is  the  mass-extinction  coefficient, 
e’  a  s 

In  comoutinq  the  SDectral  coefficients  e,  R,  and  T  from  Equa¬ 
tions  ’-23  to  2-25,  the  quantities  required  as  input  are  the  absorption 
and  scattering  cross  sections  (for  determining  <)  and  the  scattering  pat¬ 
tern  P(9)  (for  determininq  6).  ue  have  evaluated  these  quantities, 
averaged  over  specified  oarticle-size  distributions,  from  Mie  codes 
developed  oriqinally  bv  Thompson  at  General  Flectric  Co.,  TFMPO,  and 
described  in  Reference  2-lfl.  Specifically,  we  have  performed  calculations 
for  two  different,  size  distributions  of  water  droplets,  for  a  distribution 
of  dust  (basaltic  glass)  particles,  for  a  distribution  of  soot  (carbon) 
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particles,  and  for  a  combined  distribution  of  water  droplets  and  soot  par¬ 
ticles.  The  basic  inqredient  in  these  calculations  is  the  complex  refrac¬ 
tive  index  of  the  particles  as  a  function  of  wavelenqth.  Values  used  for 
the  real  and  imaginary  components  were  obtained  from  References  2-19  to 
2-21  and  are  tabulated  in  Appendix  C. 

For  a  mixture  of  soot  particles  and  water  droplets,  we  used  the 
followinq  expressions  for  the  parameters  <  and  B  that  appear  in  the 
foreqoinq  relations: 


(Pw/P5)q«  +  p* 

("pW/pS)q”  ♦  q* 


(2-34) 


(pw/pS)qV  +  qV 
"  (PW/PS)q*  +  q* 


(2-35) 


„  w,,w  ,  s,.s 
T  =  +  a 

e  e 


(2-36) 


Here,  pw  and  os  are  the  mass  density  of  water  drops  and  soot  particles, 
respect i vel v.  The  other  parameters  have  their  usual  meanina,  with  super¬ 
scripts  "s"  and  "w"  referring,  respectively,  to  soot  particles  and  water 
droplets . 


Particle  Size  nistributions .  We  assume  the  loa-normal  distri¬ 


bution  of  particle  diameters,  0: 


f  (n)  = 


-0.5'  *n(n/nm)/tnS.]; 


/?n  n  tn  5 


(2-37) 
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where  D  and  S  are  the  mean  diameter  and  standard-deviation  param- 
m 

eters,  respectively.  The  values  used  for  0m  and  S  for  each  particle 
type  are  shown  in  Table  2-6.  In  the  case  of  water  droplets  (clouds), 
calculations  were  made  for  two  different  distributions.  The  first  distri¬ 
bution  approximates  measurements  made  on  altostratus  clouds  (see  Reference 
2-22);  the  second  contains  a  larger  fraction  of  larger  droplets.  The  dust 
and  soot  distributions  are  those  used  for  the  baseline  calculations 
(Reference  2-1).  For  calculations  involvina  a  combination  of  soot  and 
water  droplets,  distribution  1  for  water  was  used.  These  distributions 
are  shown  in  Figure  2-4. 

Table  2-6.  Log-normal  parameters  used  in  the  Mie  calculations. 


Particle 

nm(pm) 

S 

water  (1) 

10.477 

1.6 

water  (?) 

17.918 

1.4 

dust  (basaltic  glass) 

0.6 

2.0 

soot  (carbon) 

0.2 

2.0 

Results 


Fiaures  2-6  to  2-8  show  the  calculated  broad-band  IR  properties 
for  water  clouds,  dust  layers,  and  soot  layers  as  functions  of  the  par¬ 
ticle  content  (am  cm-2)  along  the  path.  Figure  2-9  shows  the  broad-band 

emissivitv  of  a  layer  containing  a  mixture  of  soot  particles  and  water 

w  s 

drops  (distribution  1)  for  selected  values  of  the  p  /p  ratio.  More  com¬ 
plete  results  for  this  case,  including  the  transmissivities  and  reflec¬ 
tivities,  are  given  in  Table  C-4  of  Appendix  C. 
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OPTICAL  PROPERTIES  OF  PARTICULATE  LAYERS  IN  THE  VISIBLE 


The  formulas  2-2 3  to  2-36  have  been  applied  to  the  calculation 
of  the  optical  properties  of  Darticulate  layers  at  a  wavelength  of  0.55  urn 
in  the  visible.  The  calculations  were  done  for  a  water  cloud  (size  dis¬ 
tribution  1)  a  dust  (basaltic  qlass)  layer,  a  soot  (carbon)-particle 
layer,  and  a  layer  consisting  of  a  mixture  of  water  droplets  (distribution 
1)  and  soot  particles.  For  the  size-distribution  parameters  shown  in 
Table  ?-6,  the  values  obtained  from  our  Mie-code  calculations  for  the 
parameters  a,,  cr,  q  .  *,  and  B,  required  for  input  to  the  formulas,  are 
shown  in  Table  2-7.  Also  shown  in  Table  2-7  are  the  values  for  the  real 
and  imaginary  components  of  the  refractive  index  (n  =  nR  -  in^)  used  in 
the  Mie  calculations.  For  water,  dust  (basaltic  glass),  and  soot,  the 
index  values  were  obtained  from  References  2-20,  2-21,  and  2-1,  respec¬ 
tively. 

Table  2-7.  Input  parameters  to  optical  formulas  for  A  =  0.55  Mm. 


Part ic  le 

nR 

"1 

qa(cm2  qm'1  ) 

qs(<V  qm’1  ) 

qpUm2  qm'1  ) 

r 

6 

watPr  (1) 

1.33b 

1.0  .  Ur’ 

2. 568  -  101 

1.715  -  103 

1.741  «  103 

1.475  .  1U-7 

4.603  «  !()-? 

dust 

1.57 

4.6  .  10-“ 

8.010  «  10* 

9.22R  ■  103 

9.308  *  103 

8.605  «  10-’ 

1.090  «  10-1 

soot 

1.55 

0.1 

1.770  -  10“ 

2.652  «  lit 

4.42?  *  1ft 

4.1X13  «  10-' 

6.464  «  10-7 

Fiqures  2-10,  2-11,  and  2-12  show  the  resultinq  values  of  the 
emissivity,  transmissivity,  and  reflectivity  as  functions  of  the  layer 

o 

thickness  (qm  cm"  )  throuqh  the  water,  dust,  and  soot  layers,  respective¬ 
ly.  For  the  layer  consistinq  of  a  mixture  of  soot  particles  and  water 

w  s 

droplets,  results  have  been  obtained  for  parametric  values  of  the  P  /p 

ratio  and  the  optical  depth.  These  are  shown  in  Table  F-5  of  Appendix  C. 

w  s 

The  emissivity  results  in  this  case,  for  selected  values  of  the  P  /p 
ratio,  are  shown  in  Figure  2-13. 


=0.55  ym)  of  dust  (bas 


SECTION  3 

MODEL  CALCULATIONS 


THE  BASIC  MODEL 


Most  of  the  models  we  shall  discuss  in  this  report  will  consist 
of  two  or  three  atmospheric  layers  over  an  earth  which  is  assumed  either 
to  be  black  or  to  have  some  significant  albedo  (the  last  case  is  for  study 
of  the  effect  over  snowpack).  It  is  important  to  understand  that  within  a 
given  layer  only  the  optical  properties  of  the  medium  are  taken  as 
constant;  the  radiative  fluxes  and  the  temperature  vary  continuously 
throughout  a  given  layer.  Before  discussing  these  models  let  us  consider 
the  case  of  a  single  layer  which  we  shall  take  to  be  infinitely  deep.  We 
shall  assume  that  a  flux  of  solar  energy  of  intensity  z  is  incident  on  the 
surface  of  the  layer  which  we  take  to  be  z=0.  If  such  a  flux  is  incident 
on  a  blackbody  the  body  will  reach  an  equilibrium  temperature,  T^g  ,  given 
by 


(3-1) 


where  o  is  the  Stefan-Bol tzmann  constant.  We  shall  assume  that  the 
medium  in  the  layer  is  character i zed  by  a  visible  opacity  but  that  the 
medium  does  not  emit  in  the  visible.  We  shall  also  assume  that  the  medium 
possesses  an  opacity  in  the  infrared  and  that  the  infrared  radiation 
consists  of  an  upward  flux,  U,  and  a  downward  flux,  0,  since  the  medium 
does  not  emit  in  the  visible  the  solar  flux  is  downward  and  we  shall  call 
it  S.  The  equations  for  the  fluxes  are 


dU  s  e  U  -  eT  cT4 
dz  I  I 

dO  =  _  e  o  +  et  oT4 
dz  I  1 

2e.  oT4  =  e  T  (U+D )  +  e  S 
I  I  s 

In  order  to  produce  transport  equations  of  this  form,  it  is 
necessary  to  assume  local  thermodynamic  equilibrium,  ignore  scattering  and 
provide  an  ansatz  for  the  angular  distribution  of  the  radiation.  If  the 
radiation  fields  are  taken  to  be  uniform  in  angle  over  each  hemisphere  (up 
and  down)  but  possibly  of  a  different  amplitude  (Eddington's  approxima¬ 
tion),  the  above  equations  result  provided  that  e  is  interpreted  as  twice 
the  absorption  opacity.  The  general  solution  to  these  equations  (which  we 
shall  need  later)  is 


S  =  Le 


e  T ,  -e  z 


U  =  L  |  A(1  -  JL)  e  S  +  B(e.z+1)  +  C} 
2  £  1 


D  *  1  { -A(l  +  !l}  e"^2  +  B(e.z-l)  +  C} 
2  e  1 

s 


(3-3) 


T  5  of 


\  iA(-i 

2  £j 


-1)  e  £s  +8  e jZ  +  C} 


where  we  have  defined  x  to  be  oT4  in  order  to  simplify  the  notation, 
boundary  conditions  for  the  problem  at  hand  are 


S  z=0  =  1 


n  z=0  =  0 


(3-4) 


U  z=0  =  S  z=Q 
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For  the  temperature  we  get 


t  =  L  [1  +  JL  +  (l!  -  11) e  s‘j 

2  e  s  ei  es 


(3-5) 


Many  of  the  1R  effects  on  the  nuclear  winter  phenomenon  can  be 
understood  in  terms  of  this  equation.  If  we  find  that 


T  =  oT4bb 


(3-6) 


independent  of  z.  If  e^>  we  find  that  at  the  top  of  the  layer  the 
temperature  is  hotter  than  T0R: 


TU  *  ”4  \  b  *  Ef 


(3-7) 


The  tenperature  decreases  from  this  value  as  one  proceeds  into  the  layer 
asymptotically  approaching  a  value  less  than  T  • 


1 1  =  oj4  1  (l  +  1) 

lz=«  BB  ?  u  e 


(3-3) 


If  <  ej  the  situation  is  reversed:  Equations  (3-7)  and  (3-8)  are  still 
valid  but  they  now  predict  that  the  temperature  at  the  top  of  the  layer  is 
less  than  TgB  while  the  temperature  deep  within  the  layer  is  greater  than 


In  the  normal  atmosphere  the  solar  absorption  is  small  but  there 
is  apDreciable  IR  absorption  due  mostly  to  water  vapor;  thus  ej  is  larger 
than  and  the  stratosphere  is  colder  than  T while  the  earth's  surface 


is  warmer  than  TDD  .  The  essence  of  the  nuclear  winter  phenomenon  as  it 

DD 

has  been  previously  reported  is  that  a  layer  of  soot  is  injected  into  the 
atmosphere  whereupon  becomes  greater  than  e  j ;  therefore  the  top  of  the 
layer,  which  is  at  several  kilometers  altitude,  becomes  hotter  than  T 

DD 

while  the  surface  temperature  may  cool  to  near  TDD  or  below  (the  reason  it 

DD 

is  usually  not  below  TDD  is  that  the  layer  is  not  infinitely  thick). 

DO 


Before  considering  numerical  examples  we  need  a  slightly  more 

complicated  model  so  that  the  effects  of  the  finite  thickness  of  the 

1 ayers  may  be  considered.  We  shall  study  the  model  shown  in  Figure  3-1. 

Layer  1  is  composed  primarily  of  dust  and  water  vapor  and  we  shall  assume 

that  it  has  a  solar  albedo,  ,  an  IR  opacity  of  eT  with  a  total  optical 

depth  in  the  IR  of  Uj  ,  and  no  solar  absorption;  such  a  layer  is  probably 

there  in  realistic  cases  but  its  optional  depth  is  probably  not  very 

great.  Layer  2  is  the  soot  layer  -  the  main  component  of  the  nuclear 

winter  phenomenon;  we  shall  assume  that  it  has  a  solar  albedo  of  A2,  a 

solar  opacity  of  ec  with  a  total  optical  depth  in  the  solar  of  uc  and 
s2  s2 

a  ratio  of  e  j  /e  a  •  Layer  3  is  included  in  the  model  to  represent  a 
possible  rain-out  layer  which  may  or  may  not  be  present  in  realistic 
cases.  We  assume  it  has  no  solar  absorption  but  has  an  IR  opacity  of  eT 
with  a  total  depth  pj  due  mostly  to  water  vapor.  Linder  the  third  layer 
is  the  surface  of  the  earth  which  is  taken  to  be  a  black  body. 


Within  any  layer  the  tenperature  and  the  fluxes  satisfy  Equa¬ 
tions  (3-2)  and  so  have  the  general  solution  (3-3).  A  complete  solution 
is  easily  written  down  but  the  collection  of  formulae  is  somewhat 
tedious.  Here  we  shall  study  the  behavior  of  the  temperature  at  two 
points:  ,  the  temperature  at  the  top  of  the  soot  layer  and  T  ,  the  ten¬ 

perature  of  the  earth;  the  complete  solution  can  be  found  in  Appendix  D. 
We  have: 
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LAYER  1 


:  DUE  TO  H20,  C02,  N02 

Ai  (Albedo):  DUE  TO  DUST 


LAYER  2 


£  :  DUE  TO  SOOT 

$2 

e  :  DUE  TO  SOOT,  H20,  HYDROCARBONS,  WATER  CLOUDS 
A2:  DUE  TO  SOOT-CLOUD  MIXTURE 


LAYER  3  (RAINOUT  LAYER)  £I3: 


DUE  TO  HaO 


EARTH  (BLACK  BODY) 

m//////////////////////////////m 


Figure  3-1.  Schematic  of  the  simple  model  nuclear  winter 
which  ignores  convective  instabilities. 


T  =ii  (1-A1)(1-A2)  ll  +  I  +  u  j 

s  2  a  n 


(3-9) 


-u  -u 

T  =  I  2.  (1-Ai)(1-A2)  1.1  +  ur  +  a(l  -  e  S2)  +  (1  +  u.  )  e  2j  (3-10) 
t  2  *i 


Let  us  now  look  at  some  examples  of  using  equations  (3-9)  and 

(3-10).  We  shall  assume  that  L  is  equal  to  the  solar  constant  divided  by 

4:  l  -  3.5  x  105  ergs/cm2 sec.  If  we  assume  that  layers  1  and  2  are  miss- 

inq,  that  ur  =  1.3  and  that  the  total  albedo  is  .2  we  obtain  a  temperature 
I3 

for  of  300K.  This  temperature  represents  approximately  ambient  condi¬ 
tions  in  our  model  and  should  be  used  as  a  standard  of  comparison  for 
other  temperatures  we  shall  calculate  later.  Let  us  now  consider  a  case 
which  closely  resembles  the  NRC  baseline  scenario:  we  shall  choose  a  to 

be  .25,  the  NRC  baseline  value  for  pure  soot1,  u  to  be  2,  the  value 

$2 

attained  in  the  NRC  baseline  scenario  for  the  time  of  the  minimum  surface 
temperature;  u.  we  shall  choose  to  be  .2  based  on  the  assumption  that 

U 

there  will  be  a  few  times  10_;5g/cnr  of  water  vapor  above  the  soot  layer 
and  shall  choose  the  albedos  to  be:  Ai  =  .1  based  on  as10"1+  g/cm2  of  dust 
and  Fig.  2-11;  and  A2  =  .04  based  on  Fig.  2-12.  Equations  (3-9)  and 
(3-10)  give  Ts  =  340K  and  =  255K.  We  thus  get  a  temperature  for  the 
earth  which  is  45K  below  the  ambient  value  and  a  temperature  for  the  top 
of  the  soot  layer  (which  would  have  been  at  237K  under  ambient  conditions) 
103K  hotter  than  ambient;  these  results  are  in  good  agreement  with  the  NRC 
study. 


Of  principal  interest  to  the  present  study  is  the  variation  of 
tnese  results  with  changes  in  the  parameter  a.  In  Figure  3-2  we  show  a 
plot  of  and  as  functions  of  a  with  all  other  parameters  fixed  at  the 
values  we  have  just  used.  The  values  of  a  shown  in  the  graph  range  from 
.1  to  5;  according  to  the  NRC  the  likely  range  of  values  for  a  due  to  soot 
alone  is  .03  to  l2  .  We  believe  that  there  will  always  be  enough  water 
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vapor  present  such  that  the  value  of  a  for  the  mixture  will  not  fall  below 
.1;  to  raise  the  value  of  a  to  something  greater  than  one  would  require 
the  presence  of  absorbers  which  are  more  effective  in  the  IR  than  they  are 
in  the  solar.  The  likely  candidates  for  such  absorbers  are  hydrocarbon 
vapors,  water  vapor  and  condensed  water  (clouds). 

By  looking  at  Fig.  3-2  we  can  see  that  to  make  an  appreciable 
change  in  the  temperature  structure  we  must  change  a  by  at  least  .1  or 
more;  at  the  NRC  value  of  .25  a  change  of  .1  in  a  would  lead  to  a  3K 
change  in  and  a  17K  change  in  T  .  For  the  NRC  range  of  values  for  pure 
soot  (.l<a<l)  we  see  that  the  temperature  of  the  earth  ranges  from  251K  to 
278K  while  the  temperature  at  the  top  of  the  layer  ranges  from  41 3K  to 
272K.  While  convective  instabilities,  discussed  below,  will  alter  these 
ranges,  particularly  the  range  for  the  top  of  the  layer,  somewhat,  the 
variation  implied  by  the  uncertainties  in  the  properties  of  the  soot  alone 
is  substantial  and  studies  to  determine  these  properties  with  less  uncer¬ 
tainty  are  clearly  warranted. 

The  information  on  hydrocarbon  emissions  given  in  section  2 
along  with  the  information  on  the  IR  properties  of  the  hydrocarbons, 
especially  Fig.  2-3  makes  it  clear  that  hydrocarbon  emission  cannot  change 
a  by  an  amount  nearly  as  large  as  .1.  The  most  important  hydrocarbon  is 
probably  the  ethylene  although  the  acetylene  appears  to  be  stronger  in 
Figure  2-3.  The  reason  for  the  enhanced  importance  of  ethylene  is  that 
most  of  its  emission  lies  in  the  window  where  water  vapor  is  transparent 
while  the  overlap  between  water  vapor  and  acetylene  will  partly  cancel  the 
importance  of  the  acetylene.  Even  if  we  ignore  any  questions  of  overlap 
we  see  that  the  total  of  all  the  additional  hydrocarbons  is  only  about 
equal  in  IR  emissivity  to  the  ambient  methane.  To  estimate  the  effects  of 
these  gases  on  the  temperature  structure  we  recall  that  the  optical  depth 
of  the  visible  is  taken  to  be  2.  Since  the  soot  and  the  additional  hydro¬ 
carbons  came  from  the  same  source  (the  fires)  they  must  be  approximately 
uniform  y  mixed.  Thus  we  find  that  a  is  given  by 
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a 


optical 


dep 


h  in  the  IR 


The  IR  optical  depth  of  the  soot  is  .5;  the  change  in  a  due  to  the  hydro¬ 
carbons  is 


Aa  =  IR  optical  depth  of  the  hydrocarbons 

2 

The  IR  optical  depth  of  the  hydrocarbons  assuming  no  overlap  is  about 
2*lCr2  which  gives 

Aa,  ,  ,  =  .01 

hydrocarbons 

The  specific  optical  depths  we  have  been  using  assume  a  particular  distri¬ 
bution  of  the  soot  and  hydrocarbons  over  the  earth's  surface,  a  and  Aa 
will,  however,  be  approximately  independent  of  this  assumption  since  if 
the  hydrocarbons  are  concentrated  due  to  patchiness  the  soot  will  be 
also.  Even  if  there  were  to  exist  some  mechanism  whereby  the  hydrocarbons 
were  ten  times  as  concentrated  as  the  soot  the  effect  of  the  hydrocarbons 
would  have  only  a  marginal  effect  on  the  temperature  structure.  We  con¬ 
clude  that  the  "greenhouse"  effect  of  the  hydrocarbons  cannot  significant¬ 
ly  counteract  the  "anti-greenhouse"  effect  of  the  soot. 


Much  the  same  situation  is  encountered  in  the  case  of  C02  as  was 
found  in  the  case  of  the  hydrocarbons.  According  to  the  discussion  which 
can  be  found  in  Chapter  2  the  amount  of  C02  added  to  the  atmosphere  by  the 
fires  is  small  compared  with  the  amount  in  the  ambient  atmosphere.  Fur¬ 
thermore  the  ambient  C02  is  mixed  nearly  uniformly  over  the  portion  of  the 
atmosphere  which  is  affected  by  the  soot  so  the  blasts  and  fires  will  not 
cause  an  appreciable  redistribution  of  the  C02 .  Thus  the  IR  optical  depth 
of  the  C02  will  be  the  ambient  value  of  about  .004  which  is  not  enough  to 
change  a  by  an  appreciable  amount. 
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In  the  case  of  NO2  it  is  the  value  of  Uj  ,  rather  than  a,  which 

is  most  likely  to  be  affected.  Much  of  the  NO2  created  by  the  blasts  will 

be  carried  to  altitudes  where  it  will  remain  above  the  soot.  Figure  2-1 

shows  that  if  the  high  altitude  NO2  is  distributed  more  or  less  uniformly 

its  contribution  to  m r  will  be  less  than  the  contribution  we  have  assum- 

1 1 

ed  for  water  vapor  by  a  factor  of  twenty;  thus,  even  for  relatively  con¬ 
centrated  patches  the  NO2  will  not  have  a  major  effect  on  the  temperature. 

The  question  of  the  effects  of  water  vapor  is  much  more  compli¬ 
cated.  One  reason  is  that  the  water  vapor  has  a  much  greater  chance  of 
being  important  so  we  cannot  perform  a  quick  calculation  and  dismiss  it. 
Another  reason  is  that  the  ambient  water  vapor  is  not  uniformly  distrib¬ 
uted  over  the  part  of  the  atmosphere  where  the  soot  is  found  (it  is  con¬ 
centrated  at  the  lower  altitudes)  so  the  redistribution  of  the  ambient 
water  vapor  in  addition  to  the  addition  of  new  water  vapor  must  be  con¬ 
sidered.  A  final  complication  is  that  the  absorption  coefficient  of 
water  vapor  has  a  complicated  behavior  as  a  function  of  frequency  across 
the  IR  band  whereas  our  simple  calculations  treat  it  as  constant.  To 

correct  this  last  problem  we  should  keep  more  IR  groups  (the  absorption 
coefficient  of  water  vapor  can  be  well  approximated  by  four  IR  groups) 

which  would  considerably  complicate  our  formalism.  We  shall  instead 
observe  that  optical  depths  on  the  order  of  one  are  the  ones  important  to 
us  and  water  vapor  has  an  effective  absorption  coefficient  of  about  1 

cm2/g  in  this  range.  More  complicated  calculations  using  more  than  one  IR 
group  have  confirmed  that  this  procedure  is  sufficient  to  our  needs. 

Estimating  the  amount  of  water  vapor  which  will  be  mixed  with 

the  smoke  is  a  difficult  problem.  We  first  consider  the  water  released  by 
the  combustion  process.  According  to  the  NRC  most  combustable  materials 
generate  about  1  g-H20/g-burned ;  the  emission  factor  for  the  soot  is  about 
.02g/g.  However,  the  absorption  coefficient  for  the  soot  in  the  visible 
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(2x10‘4cm2/g)  is  2X101*  larger  than  the  one  we  are  assuming  for  water  vapor 
in  the  IR.  Thus  the  increment  in  a  due  to  the  water  of  combustion  is  only 
.0025. 

While  the  water  of  combustion  is  not  sufficient  to  significantly 
change  the  value  of  a  the  atmospheric  water  vapor  redistributed  by  the 
fire  plumes  is  much  greater.  The  NRC  estimates  that  a  total  of  40,000Tg 
of  water  would  probably  be  drawn  into  the  fire  plumes  with  an  upper-limit 
estimate  of  500,000  Tg.  Since  the  estimate  for  the  total  soot  content  is 
180  Tg  these  numbers  correspond  to  changes  in  a  of  .011  and  .14  respec¬ 
tively.  The  latter  number  is  just  at  the  borderline  where  significant 
changes  in  the  phenomenon  might  be  observed.  Even  if  the  smaller  number 
for  the  total  entrained  water  is  correct,  which  is  likely  the  case,  the 
water  vapor  will  undoubtedly  be  distributed  very  unevenly  depending  on  the 
particular  location  and  also  depending  on  the  time  of  year  of  the  event. 
In  spite  of  these  reservations  we  believe  that  water  vapor  will  not  make  a 
significant  impact  on  the  nuclear  winter  phenomenon  by  causing  an  appre¬ 
ciable  change  in  a.  Even  the  extreme  amounts  of  water  vapor  seem  to  fall 
just  short  of  being  enough. 

We  thus  come  finally  to  the  question  of  the  possible  effect  on  u 
of  condensed  water.  For  condensed  water  one  problem  we  had  for  water- 
vapor  is  removed:  the  strong  dependence  of  the  IR  absorption  coefficient 
on  frequency;  the  absorption  coefficient  of  water  clouds  is  quite  constant 
as  a  function  of  frequency.  It  does  depend  on  the  drop  size  distribu¬ 
tion.  For  the  size  distribution  of  Figure  2-5  we  find  an  absorption  coef¬ 
ficient  of  890  cm2 /g;  for  that  of  Figure  2-6  we  find  730  cm2/g;  we  have 
seen  values  as  high  as  2.2*103  cm2/g  for  naturally  occurring  clouds  quoted 
in  the  literature.  We  shall  adopt  a  value  of  800  cm2/g  for  the  moment. 
If  we  take  the  NRC's  estimate  for  the  most  likely  amount  of  entrained 
water,  40,000  Tg,  we  find  a  ratio  of  water  mass  to  soot  mass  of  220.  If 
all  the  water  condensed  out  we  would  find  a  value  of  a  of  9. 


We  thus  see  that  if  any  appreciable  portion  of  the  water 
condenses  out  it  will  have  a  dominant  effect  on  a.  It  will  also  have  a 
significant  effect  on  A2  the  albedo  of  the  cloud-soot  layer.  From  Table 
C-5  we  see  that  if  the  ratio  of  water  to  soot  is  200  to  1  the  albedo  will 
be  about  .26  while  if  it  is  100  to  1  the  albedo  will  be  about  .18;  if  the 
ratio  is  50  to  1  which  would  correspond  to  a  value  of  a  of  about  2.5  the 
albedo  is  about  .12.  (Note  that  the  optical  depth  shown  in  the  table  is 
for  scattering  and  absorption;  we  are  interested  in  an  absorption  optical 
depth  of  2  which  requires  about  10-4g/Cm2  of  soot.)  In  Figure  3-3  we  plot 
and  as  given  by  Equations  (3-9)  and  (3-10)  using  the  value  of  A2 
given  in  Table  C-5  (recall  that  convective  instabilities  were  ignored  in 
deriving  these  equations). 

Water  clouds,  if  they  are  present,  have  the  potential  to  play  a 
dominant  role  in  the  nuclear  winter  phenomenon.  Furthermore  clouds  typic¬ 
ally  do  form  in  fire  plumes  and  once  the  cloud  forms  its  presence  keeps 
the  top  of  the  cloud  cold  which  keeps  the  cloud  there.  We  believe  that 
the  calculations  just  presented  indicate  that  water  clouds  cannot  be 
iqnored  and  warrant  further  study;  their  temperatures  given  in  Figure  3-3 
cannot  be  considered  realistic  due  to  the  fact  we  have  ignored  the  effects 
of  convective  instabilities.  Convective  instabilities  will  be  triggered 
once  a  becomes  greater  than  I  by  any  appreciable  amount  (something  less 
than  1.2  for  typical  conditions).  Having  convinced  ourselves  that  a 
might  be  raised  to  such  values  by  the  presence  of  clouds  we  shall  now 
proceed  to  incorporate  convection  instabilities  into  our  calculations. 

THE  EFFECTS  OF  CONVECTIVE  INSTABILITIES 

Up  to  this  point  we  have  been  ignoring  the  fact  that  if  the 
derivative  of  the  atmospheric  temperature  with  respect  to  z  (recall  that 
we  are  taking  the  direction  of  positive  z  to  be  downward)  exceeds  some 
value  the  lower  hotter  air  will  rise  spontaneously  and  convection  will  set 
in.  We  shall  take  the  critical  value  for  the  derivative  to  be  7K/Km,  the 


dry  adiabatic  lapse  rate,  and  shall  not  in  this  report  attempt  to  correct 
for  condensation. 

From  the  formula  which  gives  the  temperature  in  the  cloud-soot 
layer  (Equation  D-2)  we  calculate  the  derivative  of  the  temperature  in 
that  layer  to  be: 


-e.  2 

.  tu  .25  e  [a  -  lje  S2 
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where  z  is  zero  at  the  top  of  the  cloud-soot  layer  and  increases  as  we 
proceed  deeper  into  the  layer.  The  maximum  value  of  Al  is  at  z=0  where 
the  value  is  dZ 


dz| z=0 
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(3-12) 


If  we  assume  that  the  total  optical  depth  in  the  solar  of  2  is  spread  over 

4  Km  of  altitude  c  will  be  .5  Km-1.  For  typical  values  of  the  other 
s2 

parameters  we  find  that  the  temperature  profile  will  be  unstable  for 
values  of  a  greater  than  about  1.2.  For  a  value  of  a  of  4  the  region 
z<3.2  Km  has  a  qradient  larger  than  7K/Km. 

The  details  of  the  way  the  unstable  temperature  profile  is 
relieved  by  convection  to  become  stable  will  depend  on  the  way  in  which 
the  instability  develops  and  on  the  character  of  the  atmosphere  above  the 
cloud-soot  layer.  Under  the  assumptions  we  have  been  making  the  tempera¬ 
ture  in  layer  1  is  nearly  constant  so  that  layer  is  quite  stable. 

If  the  cloud-soot  layer  mixes  with  the  air  above  it  the  value  of 
u  would  remain  unchanged  but  the  individual  opacities,  £c  and  er  would 


be  reduced;  that  would  lead  to  a  smaller  gradient  in  the  temperature  and  a 
more  stable  profile  as  can  be  seen  by  examination  of  Equation  (3-12).  If 
the  convection  is  confined  to  occur  entirely  within  the  cloud-soot  layer 
the  values  of  a  and  the  opacities  will  be  unchanged.  Since  our  goal  here 
is  to  estimate  the  potential  effects  of  the  convection  on  the  temperature 
profile,  and  the  assumption  that  the  convection  takes  place  entirely  with¬ 
in  the  layer  will  provide  the  maximum  effect  on  that  profile,  we  shall 
assume  that  the  convection  takes  place  within  the  layer. 

We  seek  a  solution  of  the  form  shown  in  Figure  3-4.  At  the  top 
of  the  cloud-soot  layer  there  is  a  convection  zone  in  which  the  tempera¬ 
ture  gradient  is  7K/Km.  Below  the  convection  region  are  stable  layers  of 
the  type  we  have  been  studying.  The  solution  to  the  problem  will  be 
determined  such  that  the  depth  of  the  mixing  layer  is  the  smallest  it  can 
be  and  still  allow  us  to  satisfy  the  following  requirements:  the  radia¬ 
tion  transport  equations  are  satisfied  in  all  regions;  energy  conservation 
is  produced  by  a  physically  realizable  convective  flux  in  the  convection 
zone.  This  last  requirement  amounts  to  insisting  that  the  convective  flux 
necessary  to  conserve  energy  in  the  convection  zone  must  vanish  at  the  top 
and  bottom  of  the  zone. 


In  addition  to  the  solutions  already  obtained  we  shall  need  the 
solution  to  Equations  (3-2)  except  with  the  last  of  these  replaced  by  the 
specification  of  the  temperature  gradient: 


dS  =  _  e  $ 
dz  s2 
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T  =  To  +  a  z 
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SAME  OPTICAL  PROPERTIES  AS  IN  FIGURE  1.  j 
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CONVECTION-FREE  REGION  SAME  OPTICAL  PROPERTIES 

OF  LAYER  2  AS  IN  FIGURE  1. 


LAYER  3  (RAINOUT  LAYER) 


SAME  OPTICAL  PROPERTIES 
AS  IN  FIGURE  1. 


EARTH  (BLACK  BODY) 


Figure  3-4.  Schematic  of  the  mode]  nuclear  winter  with  a 

convection  region  at  the  top  of  the  cloud-soot  layer 
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where  To  is  the  temperature  at  z=0  (the  top  of  the  cloud-soot  layer)  and 
we  shall  take  6  to  be  7K/Km.  We  define  a  new  independent  variable  to  be 


y  =  To  +  3z 


(3-14) 


which  is  just  the  temperature.  The  equations  become 
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The  general  solution  to  these  equations  is 
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Let  F  be  the  flux  of  convective  energy  in  the  convection  zone.  F  must  be 
zero  at  the  top  of  the  cloud-soot  layer  and  at  the  bottom  of  the 
convection  zone.  The  conservation  of  energy  requires 


£S2S  +  Zh  (U+D)  =  2CI2  T  +  ^ 


(3-17) 
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F  is  negative  throughout  the  convection  zone  and  is  approaching  zero  at 

the  bottom  of  it.  Therefore  df.  must  be  positive  or  zero  at  the  bn* ‘•ah  0f 

dz 

the  zone.  The  radiative  fluxes  are  continuous  at  the  bottom  of  the  zone 

so  if  df.  not  zer0  x  must  be  discontinuous  there  and  must  be  colder  at 
dz 

the  bottom  of  the  convection  zone  than  at  the  top  of  the  convection-free 

region  beneath  it;  a  situation  which  is  unstable.  To  achieve  a  stable 

atmosphere  we  must  therefore  have  df  equal  to  zero  at  the  bottom  of  the 

dz 

convection  zone  which  is  equivalent  to  having  the  temperature  be 

continuous  there.  We  now  have  five  conditions  which  must  be  satisfied  in 

the  convection  zone:  the  continuity  of  the  fluxes  U  and  D  at  the  top  of 

the  zone;  the  vanishing  of  F  at  the  top  and  bottom  of  the  zone;  and  the 

vanishing  of  ^f  at  the  bottom  of  the  zone.  The  unknown  parameters  asso- 
dz 

dated  with  the  zone  are  the  temperature  at  the  top  of  the  zone,  T0 ,  the 
constants  and  C2  in  Equation  (3-16),  the  constant  of  integration  which 
appears  when  we  integrate  Equation  (3-17)  to  solve  for  F,  and  the  depth  of 
the  layer.  The  solution  is  therefore  determined;  we  do  not  have  the 
freedom  to  impose  stability  criteria  on  the  layer  under  the  convection 
zone  and  we  must  concern  ourselves  with  whether  or  not  our  solution  will 
have  a  stable  temperature  gradient  in  that  region.  We  have  shown  above 
that  the  largest  temperature  qradient  in  the  convection-free  region  will 
occur  at  the  interface  with  the  convection  zone.  Differentiating  Equation 
(3-17)  with  respect  to  z  we  find: 


e  dS  +  e  (dU  +  dDj  =  ?e  dh  +  d2  F 
s2  dz  h  dz  dz  Iz  dz  dz2- 
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which  we  can  use  Equation  (3-13)  to  rewrite  (3-18)  as: 
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The  left  hand  side  of  this  equation  is  continuous  across  the  interface.  F 
is  negative  throughout  the  convection  zone  and  both  it  and  its  first 
derivative  vanish  at  the  bottom  of  the  convection  zone;  therefore  the 
second  derivative  is  either  negative  or  zero  at  the  interface.  We  thus 
conclude  that  the  gradient  of  the  temperature  at  the  top  of  the 
convection-free  region  is  less  than  or  equal  to  that  in  the  convection 
zone  which  is  stable  by  construction  so  the  entire  cloud-soot  layer  is 
stable. 


The  radiative  fluxes  in  the  convection  zone  are  given  by  Equa¬ 
tion  (3-16).  We  can  use  Equation  (3-13)  to  rewrite  Equation  (3-17)  as 


iL  (-  S  +  U  -  D) 
dz 


(3-20) 


which  gives 


S  +  U  -  0  +  C3 


(3-21) 


where  C3  is  an  integration  constant.  To  simplify  the  notation  we  shall 
assume  that  is  zero  and  shall  include  the  albedo  of  layer  1  into  that 
of  layer  2  by  defining  A  =  (1-A1)(1-A2).  The  boundary  conditions  at  the 
top  of  the  cloud-soot  layer  are  then: 
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solution  is  given  in  Appendix  D.  Here  we  shall  focus  our  attention  on  the 
temperature  at  three  points:  the  top  of  cloud-soot  layer  as  given  by 
Equation  (3-24);  the  temperature  at  the  bottom  of  the  cloud-soot  layer, 
which  we  call  and  which  is  given  by 
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T.  =  i  L 20 (To  +  Bd)4  +  AE  -  -l2J(e 

L  2  e  e 
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and  the  temperature  of  the  earth. 


T£,  given  by 


i  =  1  [ 2o (T0  +  S  d)4  +  AIT  U-i2  -  -iij  e  S2  +e  **S2  U  +  u  - 
L  2  e  e .  I } 

S2  I2  3 


These  three  temperatures  are  shown  as  functions  of  a  in  Figure  3-5  with 
the  assumption  that  the  cloud-soot  layer  is  distributed  over  8  Km  of  alti¬ 
tude. 


Comparing  the  results  shown  in  Figures  3-3  and  3-5  we  see  that 
the  effect  of  the  convection  has  been  to  raise  the  temperature  and 
lower  the  temperature  .  The  reason  we  do  not  show  results  for  values  of 
a  greater  than  about  5  is  that  for  higher  values  the  depth  of  the  convec¬ 
tion  zone,  d,  becomes  larger  than  the  8  Km  we  have  assumed  for  the  total 
depth  of  the  cloud-soot  layer.  Calculations  for  these  cases  would  require 
an  extension  of  our  methodology.  The  temperature  profile  from  the  top  of 
the  cloud-soot  layer  to  the  earth's  surface  is  a  stable  one.  The  inter¬ 
face  between  the  top  of  the  cloud-soot  layer  and  the  air  above  it  is  prob¬ 
ably  unstable.  The  temperatures  at  the  layer  top  are  around  250K;  the 
temperature  just  above  it  is  probably  somewhat  less  since  the  top  of  the 
layer  is  at  an  altitude  of  =9Km.  If  this  instability  is  present  it  will 
cause  the  top  of  the  layer  to  rise  which  in  turn  will  cause  the  top  of  the 
layer  to  become  colder  and  will  increase  the  surface  temperatures.  Thus  a 


temperature  profile  between  that  predicted  by  ignoring  convective 
instabilities  and  that  we  have  calculated  by  allowing  convection  only 
within  the  cloud-soot  layer  may  be  expected. 

Another  somewhat  unrealistic  feature  of  the  present  calculations 
is  seen  in  the  results  for  T^,  the  temperature  at  the  bottom  of  the  cloud- 
soot  layer.  This  quantity  is  seen,  in  Figure  3-5,  to  range  from  290K  to 
310K;  it  is  unlikely  that  there  will  be  sufficient  water  to  maintain  a 
condensation  cloud  to  such  high  temperatures  and  we  should  assume  that  the 
cloud  disappears  at  some  temperature  which  depends  on  the  amount  of  water 
present.  When  the  cloud  disappears  the  value  of  a  will  change  dramatical¬ 
ly,  probably  falling  below  1.  Our  present  methodology  can  allow  for  such 
an  effect  but  we  shall  postpone  such  calculations  until  a  more  realistic 
treatment  of  the  interface  at  the  top  of  cloud-soot  layer  has  been  devel¬ 
oped. 

We  should  perhaps  remark  on  the  comparison  between  the  effects 
of  the  cloud-soot  layers  we  have  been  studying  and  those  of  naturally 
occurring  clouds.  It  may  seem  strange  that  either  a  soot  cloud  or  a  water 
cloud  would  cool  the  surface  while  a  mixture  »f  the  two  would  warm  the 
surface.  The  source  of  the  effect  is  the  difference  in  the  albedos  of 
water  clouds  and  cloud-soot  mixtures.  Water  clouds  as  heavy  as  those  we 
have  been  studying  would  have  albedos  in  the  solar  of  .8  or  more;  the 
presence  of  the  soot  lowers  this  to  about  .2.  Thus  while  the  presence  of 
naturally  occurring  water  clouds  leads  to  an  increased  IR  opacity  for  the 
atmosphere,  which  tends  to  raise  surface  temperatures,  the  high  albedo  of 
the  clouds  tends  to  lower  temperatures  throughout  the  profile;  the  net 
effect  is  slightly  lower  surface  temperatures.  The  presence  of  the  soot 
mixed  with  the  cloud  allows  for  the  presence  of  high  IR  opacity  without  a 
high  solar  albedo. 


NUCLEAR  WINTER  OVER  SNOWPACK 


The  nuclear  winter  phenomenon  is  always  warming  for  the  surface 
if  it  occurs  over  snowpack.  By  this  we  mean  that  the  snow  will  ultimately 
absorb  more  of  the  solar  energy  directed  at  it  if  there  exists  a  soot 
layer,  with  or  without  an  accompanying  water  cloud,  than  it  would  under 
clear  sky  conditions.  The  reason  for  this  is  as  follows:  the  absorption 
of  solar  energy  by  the  soot  ultimately  results  in  a  downwardly  propagating 
IR  flux.  The  magnitude  of  the  IR  flux  at  the  surface  depends  on  the  opac¬ 
ity  structure  of  the  atmosphere  but,  as  we  shall  see  below,  it  is  at  least 
one  half  the  value  of  the  absorbed  solar  flux.  The  albedo  of  the  snow  in 
the  visible  is  .8  or  more  while  the  snow  is  very  nearly  black  in  the  IR. 
Thus  more  of  the  solar  energy  will  be  absorbed  by  the  snow  if  it  is  first 
absorbed  by  the  soot  than  if  it  propagates  unimpeded  to  the  surface. 


To  study  the  effect  over  snowpack  we  must  allow  for  an  upward  as 
well  as  a  downward  solar  flux.  We  call  these  respectively  S(j  and  Sg. 
They  satisfy  the  equations 


dSn 
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All  of  Equations  (3-2)  remain  valid  except  the  last  which  becomes: 


£slSU+SOj  +  £ilU+Dj  '  2tI 1  =  0 
the  general  solution  to  the  system  is 


Sy  «  B  ef  sz 


(3-28) 
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U  =  I  (l-a)A  e’£sZ  -  1  ( 1  +a  )  B  e^2  +  C(z+1)  +  E 
2  2 


D  =  -  1  (l+a)A  e’£sZ  +  1  (l-a)B  e£sZ  +  C(z-l)  +  E 
2  2 


t  =  I  U  -  a J i A  e  £$Z  +  B  e^J  +  Cz  +  E 
2  a 

where  A,  B,  C  and  E  are  constants.  We  shall  consider  the  same  model  shown 
in  Figure  1,  with  the  same  values  for  the  parameters  we  used  earlier, 
except  that  we  shall  assume  that  the  earth  has  an  albedo  to  solar  radia¬ 
tion  of  .8.  The  solar  fluxes  can  be  solved  for  independently  of  the  IR 
fluxes;  the  solution  is  given  in  Appendix  0.  Here  we  shall  concentrate  on 
the  temperature  at  the  surface.  The  temperature  with  all  the  layers 
present  is  given  by: 


te  -  i  Kin.,  )is„  e  52 
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where 
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where  A  is  the  albedo  of  the  earth, 
e 


To  render  these  expressions  somewhat  more  transparent  we  shall  assume  that 
layer  1  is  missing  (it  does  little  anyway)  and  that  A2  is  zero  (it  is 
always  small  and  is  very  small  for  a  soot  layer  with  no  water  cloud).  We 
get 


-2u  -y  -u  -p 

t  =Ie[1-A  e  2  +  (1-A  )(l+u  )e  2  +  a  (1-e  2 )  (1+A  e  2)J 

E  7  e  e  13  e  (3- 


(3-32) 


For  clear  sky  condition  we  get 


tEC  = 


(3-33) 


The  change  is  thus 


-2y  -u  -2y 

At E  =  1  £  ( ( 1-e  2)  i(2Ae-l)  +  a(l+Aee  2 ) |  +  y ^ ( 1 -Ae)e  2J  (3-34) 


As  long  as  A  >.5  At  is  always  positive,  even  if  a  and  u.  are  zero. 

6  u  1 3 

Perhaps  a  better  comparison  is  between  the  conditions  with  layers  2  and  3 
present  (the  nuclear  winter  condition)  with  those  for  layer  3  alone  (the 
normal  troposphere).  For  that  case  it  is  not  quite  guaranteed  that  the 
soot  cloud  will  lead  to  warming  although  it  will  be  for  most  conditions 
expected  to  be  encountered.  If  only  layer  3  is  present  we  get 


TEa‘?£  (1-Ae>  ‘2  *“I3J 

Comparing  this  formula  with  Equation  (3-32)  we  see  that 


(3-35) 


^  *  U  i(l-e'2M^)  (2Ae  -  1  -  ul3(l-Ae)  +  o(l  +  Ae  e'2pS2 )  j  j  (3-36) 
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This  quantity  is  positive  unless  a  is  very  small  and  p.  is  very  large. 

13 

The  reason  that  at  can  be  negative  is  that  the  soot  layer  reduces  the 
effectiveness  of  the  tropospheric  layer.  In  spite  of  this  At  will  be 
positive  unless  pj  is  very  large  (greater  than  3)  and  a  is  near  zero.  If 
a  water  cloud  is  present  along  with  the  soot  it  will  probably  melt  the 
snow;  if  we  assume  that  the  optical  depth  in  the  solar  is  2,  Equation 
(3-32)  predicts  melting  of  the  snow  if  a  is  greater  than  .77.  Thus  even 
if  the  effective  solar  flux  is  considerably  below  our  assumed  value,  which 
will  probably  be  the  case  since  we  are  presumably  dealing  with  winter  con¬ 
ditions  or  quite  high  latitudes,  the  values  of  a  associated  with  soot- 
water  cloud  mixtures  would  in  most  cases  be  high  enough  to  melt  the  snow. 

The  importance  of  the  fact  that  the  nuclear  winter  phenomenon  is 
warming  over  snow  pack  is  difficult  to  access.  The  fact  that  the  regions 
of  snowpack  will  be  warmer  than  normal  might  alter  normal  circulation  pat¬ 
terns.  Whether  that  would  occur  and,  if  it  did,  what  the  consequences 
would  be  are  questions  beyond  the  scope  of  the  present  study. 


SECTION  4 

CONCLUSIONS 


The  principal  conclusions  of  the  present  study  pertain  to  the 
possible  importance  of  various  IR  radiators  on  the  nuclear  winter  phenom¬ 
enon.  We  find  that  hydrocarbons,  CO2  and  NO2  will  not  have  a  major  impact 

on  the  phenomenon.  Water  vapor  will  have  at  most  a  marginal  effect  on  the 

temperature  profile  if  convective  instabilities  are  ignored;  water  vapor 
may  be  important  to  the  development  of  the  instability  which  will  exist  at 
the  top  of  the  soot  layer  if  no  condensation  cloud  is  present.  Condensed 
water  has  the  potential  to  play  a  major,  perhaps  dominant  role  in  the  nuc¬ 
lear  winter  phenomenon. 

If  water  clouds  are  present  at  late  times  they  will  probably 
alter  the  entire  character  of  the  temperature  profile,  changing  it  from 
one  which  is  hot  at  the  cloud  top  and  cold  at  the  earth's  surface  to  one 
which  is  cold  at  the  top  and  hot  at  the  earth's  surface.  The  amount  of 

water  needed  to  produce  such  clouds  is  well  within  the  bounds  of  what 

might  reasonably  be  expected  to  be  entrained  into  the  fire  plumes.  If  a 
cloud  exists  at  the  top  of  the  soot  layer  the  temperature  is  cold  there 
and  the  cloud  tends  to  persist.  If  the  top  of  the  soot  layer  is  cloud- 
free  the  temperature  is  hot  there  and  clouds  will  form,  if  at  all,  only  at 
lower  altitudes  where  they  have  a  much  smaller  effect  on  the  temperature 
profile.  This  hysteresis  effect  greatly  complicates  the  analysis  as  to 
whether  or  not  a  condensation  cloud  will  be  present.  To  analyze  the  prob- 
lem  one  must  calculate  the  development  of  the  plumes,  many  of  which  will 
be  cloud  topped,  some  of  which  may  not  be,  as  they  spread  and  merge  to 
form  the  heni sphere-wide  layer.  The  condensation  clouds  will  probably 
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affect  this  merging  process.  If  they  survive  to  later  times  they  may  have 
a  profound  effect  on  the  phenomenon.  The  models  we  have  described  in  this 
report  form  a  framework  in  which  one  may  think  about  the  problem  and  they 
provide  a  starting  point  from  which  to  develop  the  much  more  detailed 
calculations  which  will  be  necessary.  Performing  sufficiently  complete 
calculations  in  the  necessary  detail  presents  a  challenge. 
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APPENDIX  A 

BROAD-BAND  EMISSIVITY  FOR  OPTICALLY-THIN  BANDS 


If  a  band  model,  or  the  spectral  absorption  coefficient  for  an 
infrared  molecular  band,  is  unavailable,  but  the  integrated  band 
intensity.  Si,  is  known,  then  one  can  still  obtain  the  broad-band 
infrared  emissivity  for  the  band  in  the  optical ly-thin  case.  This  is 
easily  demonstrated  as  follows. 


If  e.(v)  is  the  spectral  emissivity  of  a  layer  at  wavenumber  v 
within  band  v.,  and  B(^)  is  the  Planck  function,  then  the  broad-band 
emissivity  for  the  band  is  defined  as 


/  ei(v)B(v)dv 

band 

ao 

/  B(v)dv 
o 


(A—  1 ) 


But  e.j(v)  is  related  to  the  spectral  transmissity,  T-(v),  by  the  relation 


e.(v)  =  1  -  T.(v) 


(A-2) 


However,  in  the  optical ly-thin  case, 

T.(v)  -  e"kvU  -  1  -  k~U 


( A- 3 ) 


where  kO  is  the  linear  absorption  coefficient  (atm-1  cm-1)  for  band 
and  U  is  the  path  length  (atm  cm)  through  the  absorber.  Thus,  in  this 
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e^v)  -  k^U  .  (A-4) 

Consequently,  because  the  Planck  function  generally  varies  slowly  over  a 
molecular  band,  we  can  write 

I  e  (v)B(v)dv  «  R(v.)l)  /  k~dv  (A-5) 

band  band 

where  v.  is  the  band-center  wavenumber.  Furthermore,  by  definition  of  the 
integrated  band  intensity, 

s,  =  /  ,  (A-6) 

band 


and  we  thus  obtain,  from  Equations  A-l,  A-5,  and  A-fi, 


where 


P'tv^S.jU 

7  B'(v)dv 
o 


R'(v)  =  - 

1.435V/T 

e 


( A- 7 ) 


(A-8) 


Rut,  integration  of  A-R  yields 


/  R'(v)dv  =  1.53  T4  ,  (A-9) 

o 


and  so 

e  "  1.53  B'(v.)S.U/T4 


(A- 10) 
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APPENDIX  B 

BROAD-BAND  EHISSIVITY  FOR  TWO  SPECIES  WITH  OVERLAPPING  BANDS 


If  we  have  a  layer  composed  of  two  different  species  whose 
absorption  bands  partially  overlap,  then  the  broad-band  emissivity  of  the 
layer  is  less  than  the  sum  of  the  broad-band  emissivities  for  the  indi¬ 
vidual  species.  This  result  is  proved  as  follows. 

Let  Ti(v)  and  T2(v)  be  the  spectral  transmissivities  for 
species  1  and  2,  respectively.  Suppose  that  an  absorption  band  for 
species  1  lies  between  wavenumbers  0  and  v2  and,  for  species  2,  an 
absorption  band  lies  beyond  wavenumber  v1#  where  ^  <  v2.  Then,  for 
0  <  v  <  vls  only  species  1  absorbs;  for  vi  £  v  £  v2.  both  species  absorb; 
for  v  >  \>z,  only  species  2  absorbs. 

For  simplicity,  define  a  normalized  Planck  function,  B(v),  such 

that 

00 

/  8(v)dv  =  1  .  (B-l) 

o 

Then  the  broad-band  transmissivity  of  the  layer,  composed  of  the  two 
species,  is  defined  as 

co 

T  =  J  T!(v)T2(^)B(v)dv  .  (B-2) 

o 
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But,  considering  the  absorption  regions  described  above.  Equation  B-2  can 
be  written  as 


x>2  00 

T  =  /  Ti(0)B( v)dv  +  /  T1(v)T2(v)B(G)dv  +  /  T2(v)B(v)dv  .  (B-3) 

0  V1  v2 

If  Ti  is  the  broad-band  transmissivity  of  the  layer  if  only  species  1 
were  present,  and  T2  is  the  transmissivity  if  only  species  2  were 
present,  then 


T 1  =  /  T ! ( v)B( v)dv 
0 

\>1  \>2  00 

=  /  T1(v)B(v)dv  +  /  T^vjBfvJdv  +  {  B(v)dv 

0  ”1  ^2 


(B-4) 


and 


I2  =  J  T2(v)B(v)dv 
0 

\>2  00  v  ^ 

=  /  T2(v)B(v)dv  +  J  T2(v)B(v)dv  +  J  B(v)dv  (B-5) 

^  v2  0 

Equations  B-3  to  B-5  can  now  be  combined  to  yield 

\>2  00 
T  =  ?!  +  T2  +  /  :Ti(v)T2(v)_T1(0)-T2(v) ;  B(v)dv  -  j  B(v)d^  -  J  B(v)dv 


But,  since 


TlT2  -  Tl  -  T2  =  (1-Ti )(1-T2)  -  1  ,  (B-7) 

and,  in  view  of  the  normalization  R-l,  Eauation  B-6  can  be  written  as 

v2 

T  =  Tl  +  T2  -  1  +  /  [l-T1(v)][l-r2(v)]B(v)dv  .  (B-8) 

Eauation  B-8  can  be  expressed,  alternatively,  in  terms  of  the  broad-band 
emissi vities,  by  notinq  that  both  the  spectral  and  broad-band  emissivi- 
ties,  e,  are  related  to  the  spectral  and  broad-band  transmissivities,  T, 
through  the  relation 

e  «  1  -  T  .  (B-9) 

We  thus  obtain  the  result 

v2  _  _ 

T  =  T L  +  r2  _  /  e1(v)e2(v)R(v)dv  (B-10) 

where  and  T2  are  the  broad-band  emissivities  of  the  layer  with  only 
species  1  and  2,  respectively,  present. 

This  proves  that  the  broad-band  emissivity  for  the  mixture  is 
less  than  the  sum  of  the  individual  emissivities  by  the  averaqe  cross 
product  of  the  two  emissivities  in  the  region  of  band  overlap.  Although 
we  have  made  no  attempt  to  prove  that  this  result  is  valid  for  an  arbi¬ 
trary  number  of  different  species,  we  see  no  reason  to  doubt  that  when 
band  overlap  occurs,  the  broad-band  emissivity  is  always  less  than  the  sum 
of  the  individual  emissivities. 
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APPENDIX  C 

TABULATIONS  OF  OPTICAL  PROPERTIES 


Presented  here  are  tabulations  of  the  complex  index  of  refrac¬ 
tion,  especially  for  the  infrared  (IR)  region,  used  in  the  Mie-code  calcu¬ 
lations  reported  in  Section  2.  Detailed  tabular  results  are  also  pre¬ 
sented  of  the  emissivity,  transmissivity,  and  reflectivity  of  sooty  water 
clouds,  both  for  broad-band  IR  thermal  radiation  and  for  visible  (*  =  0.55 
urn)  radiation. 

Tables  C-l  to  C-3  give,  respectively,  the  real  and  imaginary 
components  of  the  complex  index  of  refraction  (n  =  nR  -  in j )  for  water, 
dust  (basaltic  glass),  and  soot  (carbon).  Table  C-4  gives  the  broad-band 
IR  properties  of  sooty  water  clouds  for  parametric  values  of  the  water-to- 
soot  density  ratio,  pw/ps,  and  the  cloud  thickness  U(gm  cm-2).  Table  C-5 
gives  corresponding  results  for  radiation  of  wavelength  0.55  pm,  and  in¬ 
cludes  as  well  the  optical  depth,  t,  of  the  cloud.  The  log-normal  distri¬ 
bution  parameters  used  for  Tables  C-4  and  C-5  are  those  for  water  (1)  and 
soot  shown  in  Table  2-6  of  Section  2.  In  all  cases  the  temperature  was 
taken  to  be  250  °K. 
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A 

(pm) 

nR 

nI 

A 

(pm) 

nR 

nI 

A 

(pm) 

Hr 

nI 

0.200 

1.396 

1.10'7 

2.85 

1.149 

0.185 

6.6 

1.334 

3.56'2 

0.225 

1.373 

4.90"8 

2.90 

1.201 

0.268 

7.0 

1.317 

3.20'2 

0.450 

1.337 

1.02’9 

2.95 

1.292 

0.298 

9.0 

1.262 

3.99"2 

0.475 

1.336 

9.35-10 

3.00 

1.371 

0.272 

10.0 

1.218 

5.08'2 

0.500 

1.335 

1.00"9 

3.10 

1.467 

0.192 

11.5 

1.126 

0.142 

0.725 

1.330 

9.15'3 

3.15 

1.483 

0.135 

12.0 

1.111 

0.199 

0.750 

1.330 

1.56"7 

3.2 

1.478 

-2 

9.24  c 

12.5 

1.123 

0.259 

0.775 

1.330 

1.48-7 

3.7 

1.374 

3.60"3 

13.5 

1.177 

0.343 

0.800 

1.329 

1.25*7 

3.8 

1.364 

3.40*3 

15.0 

1.270 

0.402 

0.925 

1.328 

1.06'6 

3.9 

1.357 

3.80'3 

17.5 

1.401 

0.429 

0.950 

1.327 

2.93"6 

4.6 

1.330 

1.47'2 

28 

1.549 

0.338 

0.975 

1.327 

3.48"6 

4.7 

1.330 

CSJ 

ID 

w—i 

30 

1.551 

0.328 

1.00 

1.327 

2.89'6 

4.8 

1.330 

1 . 50"2 

32 

1.546 

0.324 

1.20 

1.324 

9.89~6 

5.2 

1.317 

i.or2 

38 

1.522 

0.361 

1.40 

1.321 

1.38'4 

5.3 

1.312 

9.80'3 

40 

1.519 

0.385 

1.60 

1.317 

8 . 55*  5 

5.4 

1.305 

1.03'2 

42 

1.522 

0.409 

1.80 

1.312 

1.15'4 

5.8 

1.262 

3.30'2 

50 

1.587 

0.514 

2.00 

1.306 

1.10'3 

5.9 

1  248 

6.22'2 

60 

1.703 

0.587 

2.20 

1.296 

2.89'4 

6.0 

1.265 

0.107 

100 

1.957 

0.532 

2.40 

1.279 

9.56'4 

6.1 

1.319 

0.131 

150 

2.069 

0.495 

2.75 

1.157 

5.90-2 

6.2 

1.363 

3.80*2 

200 

2.130 

0.504 

2.80 

1.142 

0.115 

6.3 

1.357 

5.70  c 

VwVV 

rv 
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Table  C-d.  Broad-band  IR  properties  of  sooty  water  clouds 


°H2n/0soot 

Vo 

( gm  cm"z ) 

Usoot 
(am  cm*  ) 

a  a  r.  sr  —  a  7 

1 . 009E-06 

1 

a  * 

000E-03 

•  •?  0  0  £  “  •’?  3 

3 . 0 0 0 E -  94. 

tr 

000E-03 

^  A  a  a  r  -  a  “t 

1 . 000E-95 

i » 

OOOE-O 2 

p  oC'OE-03 

3 . 090E-95 

900E-02 

‘  .000E-03 

i . 090E-04 

i . 

000E-01 

.  0 0  0  E  -  0  3 

5 . 090E-04 

3. 

OOOE-O 1 

A  a  A  CT  -  A  7 

1 . 900 £-03 

1  . 

000E+09 

L  .  000E-03 

2, 000E-03 

•■> 

000E+90 

■  .  0  0  0  E  -  v  3 

4 . 000E-03 

4  . 

OOOE+OO 

: . 000E-O3 

i,  000E-03 

4  . 

OOOE+OO 

L  .  0  0  0 E  -  0  3 

1 . 0  0  0  E  -  0  2 

i 

OOOE+Ol 

L . OOOE-03 

5 . 000E-02 

e; 

OOOE+Ol 

1  .  OOOE-03 

1 . 000E-91 

l . 

0O0E+02 

1  •  000E-03 

3 . OOOE-O 1 

5. 

000E+02 

L  .  000 £-0  3 

1 . 000E+09 

1 

009E+03 

L . 000E-02 

1 . 009E-0  6 

1 . 

OOOE-O A 

1 . 0  0  0  E  -  0  2 

5.000E-06 

c- 
J  • 

900E-04 

l  .OOOE-02 

1 , 099E-95 

1 . 

000E-03 

1 .000E-02 

5.000E-0S 

c* 

,  009E-03 

1  ,  O'  0  0  E  -  0  2 

1 .009E-04 

1 . 

000E-02 

I  .090E-02 

5.  OOOE-O 4 

c: 

w  t 

OOOE-O 2 

1 . 0  0  0  E  -  0  2 

1 .000E-03 

t , 

.  000E-91 

1 . 000E-02 

2. 900E-03 

2  i 

,  OOOE-O i 

1 . 9 0  0 £-02 

4 , 0  0  9  E  -  9  3 

4  , 

.  OOOE-O 1 

;  , o  o  E  -  0  2 

4 . 0  0  0  E  -  0  3 

c 

.  OOOE-O  1 

1  -  000E-02 

1  .  OOOE-92 

1 

.OOOE+OO 

1 . OCOE-02 

5  .  OOOE-O 2 

5  . 

.  9  0  0  E  +  0  0 

1 .900E-02 

1 . OOOE-O 1 

1  , 

.OOOE+Ol 

i  .  0  9  0  E  -  0  2 

3  .  OOOE-O ). 

5 . 

.OOOE+Ol 

i . 000E-02 

1 . OOOE+OO 

.  000E  +  02 

i  .  0  0  0  E  -  0  1 

i  .000E-06 

« 

OOOE - 05 

1 .0C0E-91 

3 . 009E-36 

j 

,  OOOE -05 

1 . 000E-01 

1 . 000E-95 

i 

.  OOOE-O  1 

L - OOOE-O 1 

5. 000E-95 

j 

.  OCoE-04 

t  ,  ft.  ft  ft,  c  _  n  i 

i  .  OOOE-O 

i 

,  o  0  0  E  -  0  3 

L  •  0  0  0  E  -  0  1 

3 . 009E-04 

'S 

.  OOOE -93 

_  0  0  0  £  -  0 1 

1  ,  ft  o  0  E  -  0  3 

•OOOE-92 

• . OOOE-O 1 

2 . OOOE-O 2 

2 

.OOOE- >2 

; .  oo ■:•£ -oi 

4 . 000E-03 

i 

.  OOOE -0  2 

1 , 999E-91 

£  .000E-03 

u 

.  OOOE-92 

•  .  v •'•Or  -01 

1 . 0  0  0  £  -  0  2 

1 

,  ft  ft..  .  c  -  A  1 

1  0  0  0  E  -  0 1 

3 . 090E-92 

er 

J 

.  OOOE-O 1 

i  ♦  j  'J  22  V  i 

i  .  0  0  0  E  —  0  j 

[ 

.  OOOE  +  OO 

1  •  "j  v  z.  “  0  1 

3 . OOOE-O 1 

r-j 

.  OOOE  +  OO 

L  -  0  E  -  1 

|  X  ft ft 

t  r, ft,  r:  *  .  ^ 

Emissivity 

Transmissivity 

Re+lectivity 

—  OllC-.M 

1  .  3  °  1 E  -  0  1 

2 . 4  3  2  £  -  0  2 

4 , i 93E-0 1 

3 . 414E-02 

ft  ,  4  5  9  £  -  ft  2 

a  ,  ■»  2  °  E  -  0 1 

3 . 9  4  7  E  -  0  2 

2  .  "46E-02 

■4 . 5  a  4  E  -  0 1 

1  .  7  9  —  —  —  .3  3 

2  .  J  4  0  E  -  0  2 

° . 434E-01 

o  .  4  £•  4  E  -  0  3 

ft  _  ft  p  O  C  _  ft  ft 

9 . 698E-01 

9 . 1 3  5  E  -  0  5 

3 . 0 1 1 E  -  0  2 

°.499E-01 

4 . 47-’E-0  ■ 

3 . 0  1 1 E  -  0  2 

9.699E-91 

3 . 2 2  E  -  1 1 

3.011E-02 

9 . 4  9  9  £  -  0  1 

8.349E-20 

2 . 0  1 1  £  -  0  2 

9 . 499E-01 

2.131E-23 

3 . 0  1 1  £  -  0  2 

9 . 499E-0 1 

C  ,  OOOE+OO 

3 . 0 1 1 E  -  0  2 

9.699E-01 

0. OOOE+OO 

3.011E-02 

9 . 6  ?  9  £  -  0 1 

0  .  OOOE  +  OO 

3 . 01  IE -0  2 

9 . 4  9  9  E  -  9 1 

0 .  OOOE  +  OO 

3 . 0  1 1 E  -  0  2 

9 . 499E-01 

0. OOOE+OO 

3 . 0 1 1 E  -  0  2 

2 . 24 1 E -01 

7 , 3  4  o  £  -  i 

1  .  303E-02 

4.414E-01 

3.3S0E-01 

2 . 337E-02 

7 . 9  3  4  £  -  0 1 

1 . 7  7  3  E  -  0 1 

2 . 4B6E-02 

9. 214 E -01 

3.  1  9 3 E - 0  2 

2 .*  4  S  E - 0  2 

9 . 3  4  0  E  -  0 1 

3 . 4  4  7  E  -  0  2 

2.^55E-02 

? . 4 1 8E -0 1 

9. 123E-03 

2  .  °  0  ?  E  -  0  2 

9. 4S0E-01 

2.7  4  4E-03 

2 . °28E-02 

=■ .  7G3E-0  1 

3 . 4  ^  7  E  -  0  4 

2 . 930E-02 

9. 7  0  3  £ - 0 1 

3  ,  ^  8  4  E  -  0  4 

2 , 9  3  o  E  -  0  2 

9  *  -  *  •  7  E  -  0 1 

1  G  p  O  -  _  ft  •? 

2 . 930E-02 

9 . 70  7E-01 

1.12SE-10 

2 , 930E-02 

9 . 707E-01 

ft  ,  ft  ft  ft  e  +  0  0 

; . 930E-02 

3 . 70 7E-01 

0. OOOE+OO 

2  .  ? 20E-02 

9 . 707£-Ol 

.  .  ft  ft  C  X  ft  ft 

2  ■  9  3  0  E  -  0  2 

9.‘70',E-01 

0. OOOE+OO 

2 . 930E-02 

2 . 735E-02 

o  ,  7  r.  7  c  _  ft  i 

1 . 9E4E-03 

1 . 270E-01 

3 ,o49£-01 

3 . 1 31E-03 

-  3  2  3  E  -  0  1 

7.545E-01 

I.  .  321E-02 

6.339E-01 

3.203E-01 

2 . 3  6  9  E  -  0  2 

8.123E-01 

1 . 424E-01 

2 . 5  3  0  E  -  0  2 

5  7  O  X  C  _  .’l  1 

3 . 4 1 2  E  -  0  2 

ft  _  7  ft  9  g; ..  ft  - 

G  <7  i  *7  r  _  ft  1 

1 .373E-03 

2 . 7  3  2  E  -  0  7 

9  .  4  7  7  E  —  0 1 

4 . 2S1E-03 

2 . 900E-02 

9 . 714E-01 

9 . 0  2  3  E  -  0  4 

~i  t  3  ft  2  £  -  0  ' 

°  .  3 1 9  E  -  0  1 

4 . 219E-05 

2 . 3  0  2  E  -  0  2 

3  .  “  2  0  E  -  0  1 

j  ,  J  1  *2  ‘  .)  / 

2.  3  ■  1 2  E  -  0  2 

9 . 3  2  0  E  -  0 1 

2 , 302E-02 

-j  a  c  _  4 

r,  r  2  a  ■' 

"■  o r  _  ft  -I 

■'  *  ■>  - 

o  #  7  ^  0  £  -  0  i 

0. OOOE+OO 

ft  _  ft  ft  ft  C  _  ft  ft 

3  f  2  0  L  —  v  1 

0  ■  OOOE  +■•..; 

2  .£•:  2  E  -  0  2 
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Table  C-4.  Broad-band  IR  properties  of  sooty  water  clouds  (continued) 


PH20^PSOOt 

uh2o 

(qm  cm*  ) 

Usoot 
( qm  cm*  ) 

5 .OOOE-Ol 

1 . 000E-O6 

2.  OOOE -06 

5 . OOOE-O 1 

5 . 000E-C6 

1 . 000E-05 

3 . 000E-01 

1 .000E-05 

2. 000E-05 

5. OOOE-Ol 

5. OOOE -05 

1 .000E-0  4 

5 . OOOE-Ol 

1 . 000E-04 

2 . 000E-04 

5. OOOE-Ol 

5 . 000E-04 

1 . OOOE-03 

3 . OOOE-Ol 

1 . 0  0  0  E  -  0  3 

2. 000E-03 

5 . OOOE-Ol 

2.000E-03 

4. 000E-03 

2 . OOOE-Ol 

4 . OOOE -03 

3.  OOOE-03 

3 . OOOE-Ol 

£  .  000E-03 

1 . 200E-02 

3. OOOE-Ol 

1 . OOOE-02 

2 . OOOE-02 

2 . OOOE-Ol 

5 . 000E-02 

1  .  OOOE-Ol 

S. OOOE-Ol 

1  .OOOE-Ol 

2. OOOE-Ol 

5. OOOE-Ol 

5  .  OOOE-Ol 

1 .OOOE+OO 

3. OOOE-Ol 

1 .OOOE+OO 

2. OOOE+OO 

1 . OOOE+OO 

1 .000E-06 

I . 000E-06 

1 , OOOE+OO 

5 . 0  0  0  E  -  0  6 

5 . 000E-06 

1 .OOOE+OO 

1 . 000E-05 

1 .000E-05 

1 , OOOE+OO 

3 . OOOE-05 

5. 000E-05 

1 , OOOE+OO 

1 . 000E-04 

1. 000E-04 

1 .OOOE+OO 

5 . 000E-04 

5. 000E-04 

1 . OOOE+OO 

1 .000E-03 

1 . OOOE-03 

1 .OOOE+OO 

2  .  OOOE-03 

2.  OOOE-03 

1 .OOOE+OO 

4 . 000E-03 

4.  OOOE-03 

1 .OOOE+OO 

6 . 000E-03 

£  .  000E-03 

1 . OOOE+OO 

1 . 000E-02 

1 . 000E-02 

1  .  0  0  0  E  +  0  0 

5 .000E-02 

5. 000E-02 

1 , OOOE+OO 

1  .  OOOE-Ol 

1 . OOOE-Ol 

1 . OOOE+OO 

5  .OOOE-Ol 

3. OOOE-Ol 

1 . 0  0  0  E  +  0  0 

1 . OOOE+OO 

1 . OOOE+OO 

3 . OOOE+OO 

1 . OOOE -06 

2. OOOE-O? 

3 . OOOE+OO 

3 . OOOE -06 

1  .  OOOE  -06 

3 . OOOE+OO 

1 . OOOE-05 

2  .  OOOE -0 6 

3 . OOOE+OO 

3  .  OOOE-O 3 

1 . OOOE -03 

3. OOOE+OO 

1 . OOOE-O 4 

2. 000E-03 

3 . 0  0  0  E  +  0  0 

3  .  OOOE-04 

1 .000E-04 

3 . OOOE+OO 

1 , OOOE-03 

2. 000E-04 

3. OOOE+OO 

2 .OOOE  -03 

4  ,  OOOE- 0 4 

3 . OOOE+OO 

4 .OOOE-03 

3 . OOOE-O  4 

3 . OOOE+OO 

6 . 000E-03 

1 . 200E-03 

3 , OOOE  100 

1  .  000E-02 

2. OOOE-03 

^-00  0 E  1 0 0 

*  0  0 OE  -  0 

1 .000E-02 

3 . OOOE+OO 

: . oooE-oi 

2.  OOOE-02 

3.00  0 E  +  0  0 

-  t  r>  Q  c:  -  r, 

1 . OOOE-O L 

3 . OOOE+OO 

1  .  OOOE  +  OO 

2. OOOE-Ol 

Emissivity 

Transmissivity 

Reflectivity 

6.254E-03 

9 . 933E-01 

4 . 361E-04 

3.079E-02 

9.671E-01 

2.092E-03 

6 . 035E-02 

9 . 337E-01 

3 . 9+7E-03 

2 . 594E-C 1 

7.267E-01 

1  .  3  9  4  E  -  0  2 

4 . 363E-01 

5 . 4  4  1  E  -  0  1 

1 . 966E-02 

3. 61  IE- 01 

1 . 095E-01 

2 . 640E-02 

? . 362E-01 

3.661E-02 

2 . 7  2 1 E  -  0  2 

9.631E-01 

9 . 350E-03 

2  .  •’  3  1 E  -  0  2 

9 . 71 4E-01 

1 . 0  6  0  E  -  " 

2  .  7  3  5  E  -  0  2 

9 . 7  2  3  E  -  0 1 

1 . 345E-04 

2 . 755E-02 

3 .724E-01 

2 . 361E-06 

2.  -'55E-02 

o  t  7  2  4  E  -  0  1 

1  .  “» 7 1 E  -  2  3 

2  ■  732E-02 

9.724E-01 

0 . OOOE+OO 

2  .  ?  5  5  E  -  0  2 

9 . 72 4 E -01 

0.  OOOE  +  OO 

2.  ’53E-02 

9.724E-01 

0  .  OOOE  +  OO 

2  •  7  5  5  E  -  0  2 

3.572E-03 

o  ♦  °  6  2  E  -  0  1 

2 . 3  £  1 E  -  0  4 

1 . 7  7 1 E  -  0  2 

9.811E-01 

1 . 155E-03 

3 . 504E-02 

9 . 6  2  7  E  -  0  1 

: . 2 1*2-03 

1 .612E-01 

3 . 2  °  6  E  -  0  1 

201E-03 

2 . 9  1 3  E  -  0  1 

6.935E-0: 

1  . a  7  7E - 02 

7 . 63SE-0 1 

2 . 053  E - 0  1 

.  5  4  1 E  -  0  2 

9 . 0 6 2 E  -01 

h , 7  o  3  f.  - :: 

2  •  7  5  E  -  0  2 

9 . 3S7E-0 1 

1  .  4 13E - 02 

:  .  7  1  4  £  -  0  2 

9.714E-01 

1  •  4 1'E -0  7 

*  .  ^  1  •- c  -  0  .* 

0  .  '7  2  £  E  -  0  1 

1 . 370E-04 

.  720E-o: 

9.72SE-01 

?  .  5  4  0  E 

.  .  '20E-0:. 

9, 72SE-01 

f. ,  1  4  *5  ff  -  2  7 

2  .  ■'  2  0  £  -  0  2 

'?  .  7  ?  9  F  -  G  1 

v  .  OOOE  A  ■' 

2 . '  'OE  -02 

o  t  "■  2  8  E  -  0  1 

o .  o  oo  £+•:••■> 

2  .  7 2 0 E  -02 

°*728E-C1 

0  ♦  •  0  E  + 0 

- .  “205  -02 

1  . 419E-03 

?  ,  °  9  j  E  -  '•  l 

~  .  )  1 7  E  -  0  5 

7 . 073E-03 

?  .  *  2  6  E  -  I 

e  -  o  •> 

1  .  4  0  9  E  -  0  2 

Q  f  C  ..  1 

7  .  •  6  0  E  -  0  4 

6 . 3  4  5  E  -  0  2 

-•  ■  281E-01 

7  .  n  3  £  E  -  0  3 

1 . 321E-01 

3.613E-M 

i  .  0  7  E  -  0  3 

3 . 0  2  7  E  -  0  1 

4  .  7  3  3  E  -  ■  1 

1  .  392E-02 

7 . 4  3  7  E  -  0 1 

2 . 3  2  IE -01 

•  313E-02 

9 . 1  ?  1  E  -  0 1 

*  .  7  4  2  E  -  *. » 2 

: .  i :  3  e  -  o  2 

9 . 703E-01 

4 . 3  6  3  E  -  0  3 

2  .  4  £  3  E  -  0  2 

9 . 749E-01 

4 .63  IE -04 

.  4  6  F  -  0  j 

9 . 1 jIE-01 

7  .  310E- :•£ 

2 . 4 iiE-02 

9  .  ?  3  3  E  -  0 1 

-•  .  -?  i  r-  c  _  p 

2 . 1 6  -  0  2 

°.,53E-01 

0  .  V  0  E  r  0  0 

.  *.  i  i  E  -  0  2 

■-  .  7  2  3  E  -  0  1 

0 . OOOE+OC 

7  .  ■'  -  6  E  -  0  2 

•-  . 7S3E-0i 

0  .  0  0  C  +  0  0 

;■  .  4  £  £  -  o  2 
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Table  C-4. 


Broad-band  1R  properties  of  sooty  water  clouds  (continued). 


(3HJt/<>S00t 

So 

(qm  cm*  ) 

U500t 
(qm  cm*  ) 

Emissivity 

Transmissivity 

Reflectivity 

1 . OOOE+Ol 

1 . 000E-06 

1 . 0  0  0  E  -  0  7 

1  .  1  19E-02 

?  «  9  8  3  E  -  0  \ 

5  .  •) 3°E  -05 

1  .  0006+01 

5 . 000E-06 

3 . 000E-07 

3. 733E-03 

9 . 9406-0  1 

2  •  7  0  7  E  -  0  4 

1 . 000E+01 

1 . 00CE-05 

i  .  0  0  0  £  -  0  6 

1 . 1436-02 

9-8806-01 

5 . 3  315-04 

i  .  OOOE+Ol 

5 , 0  0  0  E  -  0  5 

3 . 0006-06 

5,5846-02 

9. 4166-01 

2 .5*56-03 

1  .  0  0  0  E  +  0  1 

1  .OOOE-C4 

i .oooe-05 

1 .0856-01 

3.3666-01 

’  .3  +  06-03 

1 .OOOE+Ol 

5.00 0  £ - 0  4 

5. 0006-05 

4 . 34<?E-01 

5 . 4926-01 

1 . 38BE-02 

1  .  OOOE  +  Ol 

1 . 0  0  0  E  -  0  2 

1 .0006-0  4 

6 . 7606-01 

3 . 0  3  3  E  -  0 1 

,  0  ^  E>  £  —  0  2 

1 .OOOE+Ol 

2 . 000E-03 

2 . 0006-04 

8.3336-01 

9 . 3936-02 

2  .  2756-02 

1. OOOE+Ol 

1 . 000E-O3 

4 .  0006-04 

9. 6 7 4E -01 

°  .  E  7  4  E  -  0  3 

2 . 30 26 -02 

1 .OOOE+Ol 

*6 . 0  0  0  E  -  0  3 

6 . 0006-04 

9.7596-01 

1 .05  46-03 

2 . 3036-02 

1  -  OOOE  +  Ol 

1 .000E-02 

1 .0006-03 

e . 7706-01 

l . 5  486-05 

2 .3036-02 

1 . OOOE+Ol 

5 . 0  0  0  E  -  0  2 

5.0006-03 

9. 7706—01 

4.136E-22 

2 . 303E-02 

1  .OOOE  +  Ol 

1  .  0  0  0  E  -  0  1 

1 .0006-02 

9.7706-01 

0.  OOOE+OO 

2 . 3O3E-02 

i - OOOE+Ol 

5.000E-01 

5.0006-02 

9.7706-0! 

O.OOOE+OO 

2 . 303E-02 

1.  .  ■}  0  0  F  +  0 1 

1 . 000E+00 

1 . 00 OE -0 1 

9 . 7706-01 

0 . 0006+00 

2  •  3 0 36-02 

3. OOOE+Ol 

1 . OOOE-06 

2.0006-09 

9.33 56-04 

9. 9906-01 

3 .3166-05 

3. OOOE+Ol 

5 . 000E-06 

1  .  00 OE  -0  ■* 

4  .  6336  -0  1 

9.9526-01 

1 . ° 006 -04 

3. OOOE+Ol 

1 . 0  COE- 0  5 

2  .  0  0  0  E  -  0 

9 . 29-16-  0  7 

J  .  ?  a  3  c  A  t 

2.7526-04 

3 . OOOE+Ol 

S  .000E  - 03 

i 

.  3  3  ;  E  -  '  2 

9  2  2  i  E  -  O'  1 

1 . 3.196-0  3 

3 . OOOE+Ol 

1  .  0  0  0  E  -  O'  •> 

2  e  -  -- 

-  A  *> 

-  .  >  -  E  E  •  1 

3 .  1 696-0. i 

j  .  0  0  0  E  +  0  i 

5  .  0006-04 

«  t  ^  ^ 

1  .  1  •  _  £  -  ‘j  ’ 

«  7  -  IT  .»  • 

1  .  2  -  E  -  0 ,2 

3. OOOE+Ol 

1 . 0006-03 

^  ,  ;«.)  ■  f  . 

-  .  *  ‘  L  '  "•  ! 

2  2  m  5  r  -  1 

:  .  9  4  E  -  0  2 

: . oooE+oi 

2  ,  A  A  A  t:  _  A  7 

•i  ■  7 

-  .  *  E  ~  '  1 

'  /  ~  c  r  - .-  i 

1.9.  It 

j  ♦  0  0  0  E  x  0  1 

•>  •  0006  -  0  3 

3. : 

:  #  ;  c  .  '  • 

3  .  5  -  “2  E  **  2 

: . 0196-02 

o.oooe+o: 

r  .  vC-OE  -  '• 7 

:  2  ■:  ■' 

. 

r~  ■  -  /  > 

^  r‘  -  0  2 

3  .  )  0  0  E  •+  0  1 

i . oooe -o: 

'.£-••  4 

;  ■  r  r  .  ,m 

2  .  3  ^  r  „  >  .4 

-  .  A  -<  -  r  ..  a  -) 

r  a  a  r;,'  + 

r-  -  -  -  - 

-  -  s  _ 

-  *- 

*  .  .  . 

»-  .  .  w  ...  .  . 

..  »_  .  “* ' )  _ 

3  *  0  0  0  E  +  0  1 

l  ,  2  •>  •;  £  -  j  ; 

J  ,  ; 

.  -M  ■ 

1  .  '  c  2  r  -  2  ’ 

C-,11 

r*  •  • 

- 

,  ^ 

^  A  ..  ,  ^ 

-  -V  *■ 

*•**■'■-  -  * 

•  «. 

-  .  .  -  _ 

...  ■’  *.  .  Z-  ~ 

oooE+oi 

1 . 0006+ 0 ’ 

3  -  '  %  ~  2 

■  b  '• .-  a  r  »  • 

T  •  ■'  C  ^  .'  "> 

:  7  r  .  •  <  *.  ^  •  ••  l-  i  '  -.  r  ..  .  c  . 
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Table  C-4.  Broad-band  IR  properties  of  sooty  water  clouds  (continued) 


m 

fM 

$ 


pH20,/Bsoot 

uh20 

(gm  cm-  ) 

Usoot 
(gm  cm*  ) 

Emissivity 

2 . 000E+02 

1 . OOOE-Oa 

5 . 000E-09 

8  .  ?  3  0  £  -  0  4 

2 . 000E+02 

5 . 000E-06 

2 . S00E-08 

4 . 437E-03 

2 . 000E+02 

1 . 0Q0E-05 

3 . OOOE-03 

9  *  992F-0"7 

2. OOOE+02 

5  .  OOOE-05 

2 . 5  0  0  E  -  0  7 

4 . 364E-02 

2  .  OOOE+02 

1 .000E-04 

5.000E-07 

8 . 3  2  9  E  -  0  2 

2 . 0002+02 

5  .  OOC'E-0  4 

2 . 5  0  0  E  -  0  3 

3 . 569E--0  1 

2 . 000E  +  02 

1  .  OOOE-03 

5 . 0  0  0  E  -  0  6 

5 . 3  0  9  E  -  0  1 

2 . 000E+02 

2  .  OOOE-03 

1 . 000E-05 

8. 1 1 7  E - 0 1 

2 . 000E+02 

4 .000E-03 

2 . 000E-05 

9 . 462E-01 

2. 000E+02 

6.  OOOE-03 

3. 000E-05 

9.724E-01 

2 . 000E  +  02 

1 . 000E-02 

5. OOOE-05 

°  .  800E-0 1 

2.00CE+02 

5 . 000E-02 

2 . 500E-04 

9.307E-01 

2. 000E+02 

1  .OOOE-Ol 

5 . 0  0  0  £  -  0  4 

9 . 307E-0 1 

2 . 000E+02 

5 . 0O0E-01 

2.500E-03 

9. 907!“ -  01 

2 . OOOE+02 

l .OOOE+OO 

5. 000E-03 

9.307E-01 

3 . 0QQE+02 

1 . 000E-06 

2 . 0  0  0  E  -  0  7 

8 . 3  4  9  £  -  o  4 

3. 000E+02 

5.000E-0A 

1 . 000E-08 

4  .  4 1 6  E  -  0  3 

5, 000E+02 

1 .000E-05 

2 . 000E-08 

3 . 312E-03 

3. 000E  +  02 

5.  OOOE-05 

1 . 0  0  0  E  -  0  7 

4 . 3  2  3  E  -  0  2 

3 . 0002+02 

1 . 0  0  0  £  -  0  4 

2 . 0  0  0  E  -  0  7 

8 . 4  8  3  E  -  0  2 

3. 000E  +  02 

5 .000E-04 

1 . 0  0  0  E  -  0  6 

3 . 3  4  2  E  -  0  1 

3. 000E  +  02 

1 .000E-03 

2 . OOOE-Oa 

3 . 7  7  3  £  -  0 1 

3, OOOE+02 

2. 000E-03 

4  .OOOE-Oo 

q  a;  v  r  _  a i 

5 . 0  0  0  E  +  0  2 

4  ♦  0  0  0  £  -  0  3 

8 . 0  0  0  E  -  0  i 

9 . 4  4  c-  E  -  0  1 

3 . OOOE+02 

6 . 000E-03 

1 . 200E-05 

3 . 713E-01 

3 . 000E+02 

1 . 000E-02 

2 . 0  0  0  £  -  0  5 

- . 30  IE -01 

3 , 000E+02 

3 . 000E-02 

1 . 000E-04 

•=>  .0,09  c  ~  a  i 

3 , 000E+02 

1 . OOOE-Ol 

2.00  0  E - 0  1 

■?  „  O  A  9  C  _  A  ! 

3 . 000E+02 

5. 000E-01 

1. 000E-03 

° . 8  0  9  E -01 

3 . 000E+02 

1  .OOOE  +  OO 

2 . 000E-03 

3.809£-v: 

3 . 0G0E+03 

1 .000E-06 

1 . 0  0  0  E  -  0  ° 

3 . 322E-0  4 

1 . OOOE+03 

3 . 0  0  0  E  -  0  6 

5 . 000E-0? 

4  .  4  0  3  E  -  0 .3 

1  .  0  00E  +  03 

1 . 000E-05 

1 . 0  0  0  E  -  0  8 

•3.7  23E-03 

1 . OOOE+03 

3 . 000E-03 

5 . 000E-08 

4.31 2  E - 0  2 

*  .  J  0  0  E  +  0  i 

1 . 000 E -04 

1  -  0  0  0  E  -  0  7 

3. 4  3  0  E  -  0  2 

1  *  0  0  0 17  r  •'“*  7 

3.000E-04 

3 . 0  0  0  E  -  0  " 

3 . 3  3  3  E  -  0  1 

\  '}AAC  i.’ji 

1 . 000E-03 

1 . 000E-06 

8 . 3 i OE -0 1 

1 . 000E+03 

2 . OOOE-03 

2 . 0  0  0  8  -  0  s 

.  0  '’l  c.  .-.  i 

1 . 000E+03 

4 . 000E-03 

4 . OOOE-Oo 

3  .  4  4  1  E  -  >  1 

1 .00  0  E  f  A  7 

£  .  0  0  0  £  --  0  3 

i . OOOE-Oa 

v  .  7  ■  ,4.,p  i 

1  .  0  0  0  E  t  0  3 

1 . 000E-02 

1 . 0  0  0  £ -03 

3 . 30  IE -01 

1 . 0C0E-03 

7 . 0  0  0  E  -  0  2 

Z  .  v  0  O' c  -  •}  r 

'J.yAJu  .  ^  • 

1  v  0  0  £  -f  0  3 

1 . OOOE-Ol 

1 . 0  0  0  E  -  0  4 

-  ,  <:•  1 

1 . 000 Etc-  ’ 

3. OOOE-Ol 

3 . OOOE-Ol 

:  .  j.iOl'-'i 

1  .  O'OOE  r O'  3 

1 . 000E+0 0 

1 • 000E-03 

- .  9oc'i>o : 

Transmissivity  Reflectivity 


'*•'  *■-  ■*.  ■■". 
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3.Q91C-01 
■f .  ?  5  4  E  -  0  ' 
9  .  ?  0  8  E  -  0  1 
?  .  E  4  7  E  -  0  1 
9 . 1 15E-0  1 
•i  .  3  1  '  E  -  ()  1 
4. 032E-01 

3 . 4  4  3  E  -  0  2 
2 . 2  , '  S  E  - .'  3 
6.592E-04 

4  .  3  6  3  E  -  2  0 

0 . 0  "  0  c-  *  a  a 

O.OOOE-O > 

?  .  o.j  iE:  ..0  • 

3 . 954E-01 

9.351E-0L 

:••  34  iE  -01 

4  .  o  ’ :  3  -  •  : 

1  .  3  3  3  E  -  ■ 1 

i  .  -  3  ■'■  3  -  A  ' 

5  .  *■  4  2  E  -  0  3 

4  •  03'E-l 1 

t  ■■  r,  ,j  r  ^  -,  a 

9  _ov  tr_,M 

9  -i  •,  9 1.  -  '.  • 


.  •  9  .*  r  _  . .  < 

*  -  1  «  C  ..  .  . 

»  •  1  -  . 

2  •  :  3  3  E  -  0  2 

-  .  3 .:  4  E  -  :  : 

? .  23-::-  ■  i 


T  -  1  '  (T  _  " 

— '  *  — ‘  1  ~  *-  V  fcJ 

t  .  3  4  9  E  -  0  4 
7  .  4  3  1  E  -  0  4 
1 . 676E-03 
3 . 200E-03 
1  .  1332-02 
1  .  2  2  4  E  -  0  2 

c-  O  c  —  A  ■' 

l  .  2  S  E  -  0  2 
1  -  ?  2  4  E  -  0  2 
1  .  3  3  3  E  -  0  2 
1  .  935E-02 
1 .  °3SE-02 
1  .  y  222-02 
1  .  9  3  5  E  -  0  2 

*  .  4  T  '■  r‘  -  0  L~ 

1  ,  713S-04 

2  .  4  2 1  E  -  0  4 
:  .  :■  4  7  E  -  0  7 
. '  ■  1  •’  2  2  -  0  2 
1  .  1  V9E-02 

1  .  -  ^  rr  ■■ 

1 . 2  3  5  E  -  0  2 
I.  .  9  0  3  E  -  0  2 
1  .  '5 1  42-02 
1  •  9  1  5  2  -  0  2 
i  .  317,E-02 

*  -  !.  l! 

: . ,J  1 5  E  -  0  2 

-■  ■  •■-.  C  ..A  '-■ 


3 . 4  :■ :  h  - ; 


I  ;  '  3  2  -  o  - 
1  3  Z  3  3  -  ■ '  2 

1 . 32  7£  -v  7 

1  .  V  0  1  E  -  0  2 

1  •  •  v  3  E  -  :• .. 

1  .  '  ■  ■ .  c  -  ■■;  2 

:  . -  2  32-02 


*.\v 


-*.-■74 


;  r\* 

ii-jia 


\  \%  .V 


'  j"7 


«  \  « 

-  •V  *. 

. . .  .  % 

•  /  •/  -.y/_ *.v»  -■,v;2.v^lv.v.>'l%vSVtJVJ^ 
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Table  C-4. 

0Hjf/B$OOt 

Broad-band 

'jh2o 

(am  cm'  ) 

IR  properties 

u 

soot 

( am  cm-i ) 

of  sooty 

Emissivity 

*  ■  %*. 

water  clouds 

Transmissivity 

(continued) 

Reflectivity 

.3  .  v  0  v  E  4-  0  5 

1  •  0  v  0  E  -  0  d 

?  .  i  0  E  -  0  1 

0  (  0  9  i  r  „  <\  ^ 

3 . 4  1  4  E  -  0  5 

':  .  000E+03 

5 . 000E-06 

1 . 000E-09 

4 . 392E-03 

9 . 0  S  4  E  -  0  1 

!.  .  7  0  0  E  -  0  4 

3  .  ■“•  0  0 c  A  •'  7 

1  .OOOE-OS 

2 . 0  COE- 0  4 

8 . 7  6  4  E  -  0  3 

0 . 9  0  9  E  -  0  1 

3 . 3  8  5  E  -  0  4 

i  . .  :•  •. 

0  .  OOOE-OS 

1 . OOOE-09 

4 . 2  0  2  E - 03 

9.553E-01 

1 . 630E-03 

i  Z  ■  '  0  }  ~  ■*■  •}  7 

: .  ;ooe-C4 

2 . OOOE -03 

S. 4 i 0  E - 0  2 

9 . 1 2  8  E  -  0 1 

3  .  1 12E-03 

-  •  -  w  "  »'  J 

5  .  j  0  0  E  -  0  4 

I  .  0  0  0  E  -  0  7 

3 . 3  2  8  E  -  0  1 

6 . 3  8  4  £  -  0  1 

•1  t  <  0  8  E  “  0  2 

’*<■'*  •'  V  —  ■'  T 

3 .  •:  •:  -j  e  -  0  “* 
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Table  C-4.  Broad-band  IR  properties  of  sooty  water  clouds  (Concluded) 


D  /  0 

H2O  soot 

Vo 

(gm  cm"  ) 

Usoot 
(gm  cm"  ) 

Emissivity 

Transmissivity 

Reflectivity 

i .OOOE+IO 

1 . 0  0  0  E  -  0  6 

1  .OOOE-16 

3.795E-04 

?  .  Q  9  1 E  -  0 1 

3 

. 410E-05 

t . OOOE  +  IO 

5 . 0  0  0  E  -  0  6 

5 . 0  0  0  E  -  1  6 

4 . 389E-03 

9 . 954E-01 

1 

.  6  9  8  E  -  0  4 

1 . OOOE  +  IO 

1 . 000E-05 

1 .000E-15 

3 . 759E-03 

9 . 909E-0 1 

3 

.  3  9 1  £  -  0  4 

1 . OOOE+IO 

5 . 000E-05 

5 . 000E-15 

4 . 299E-02 

9 . 554E-0 1 

1 

.  6  2  3  E  -  0  3 

1  »  OOOE  +  IO 

1 . 000E-04 

1 .000E-14 

3 . 405E-02 

9 . 1 23E-0 1 

3 
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1 . 0  0  0  E  +  1 0 

5 . 000E-04 

5 . 0  0  0  E  -  1  4 
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1 

. 550E-02 
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.  9  0  0  E  -  0  2 

1 , OOOE+IO 
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1 . 0  0  0  E  - 1 2 
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3 . 3  5 1 E  -  0  4 
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Table  C-5.  Optical  properties  in  the  visible  (X=0.55  wn)  of  sooty 
water  clouds. 


a  /o  T  *  'soot  ^H20  Q  T  £ 

H20'  soot  /not'cal  Oeoth)  'qm  an-2)  1  qm  (Pef lectivity)  (Transmissivity)  ( Emissi vi ty) 

l . 000E-03  4.422E+01  1.000E-03  1.000E-06  4.424E-02  3.986E-09  9.558E-01 

1 . 000E-03  2 . 21 1E  +  02  5.000E-03  5.000E-06  +.42+E-02  O.OOOE+OO  9.558E-01 

1 . 000E-03  4 . 422E+02  1. 000E-02  1.000E-05  4.424E-02  O.OOOE+OO  9.558E-01 

1 . 000E-03  2.21 1E+03  5.000E-02  5.000E-05  4.424E-02  O.OOOE+OO  9.558E-01 

1 . 000E-03  4 . 422E+03  l.OOOE-Ol  1.000E-04  4.424E-02  O.OOOE+OO  9.558E-01 

1 . 000E-03  2.21 1E+04  5.000E-01  5.000E-04  4.424E-02  O.OOOE+OO  9.558E-01 

1 . 000E-03  4 . 422E+04  l.OOOE+OO  1.000E-03  4.424E-02  O.OOOE+OO  9.558E-01 

l.OOOE-03  8 . 844E+04  2.000E+00  2.000E-03  4.424E-02  O.OOOE+OO  9.553E-01 

1.000E-03  1 . 769E+05  4.000E+00  4.000E-03  4.424E-02  O.OOOE+OO  9.558E-01 

1 . 000E-03  2 . 653E+05  6.000E+00  6.000E-03  4.424E-02  O.OOOE+OO  9.558E-01 

1 . 000E-03  4 . 422E+05  l.OOOE+Ol  1.000E-02  4.424E-02  O.OOOE+OO  9.558E-01 

1.000E-03  2.21 1E+06  5.000E+01  5.000E-02  4.424E-02  O.OOOE+OO  9.558E-01 

1.000E-03  4 . 422E+06  1.000E+02  l.OOOE-Ol  4.424E-02  O.OOOE+OO  9.55SE-01 

1.000E-03  2.211E+07  5.000E+02  5.000E-01  4.421E-02  O.OOOE+OO  9.558E-01 

1 . 000E-03  4 . 422E+07  1.000E+03  l.OOOE+OO  4.424E-02  O.OOOE+OO  9.558E-01 


1.000E-02  4 . 424E+00  1.000E-04  1.000E-06  4.333E-02  1.443E-01  8.124E-01 
1 . 000E-02  2.212E+01  5.000E-04  5.000E-06  4.425E-02  6.305E-05  9.557E-01 
1 . 000E-02  4 . 424E+0 1  1.000E-03  1.000E-05  4.425E-02  3.934E-09  9.557E-01 
1 . 000E-02  2.212E+02  5.000E-03  5.000E-05  4.425E-02  O.OOOE+OO  9.5S7E-01 
1 . 000E-02  4 . 424E+02  1.000E-02  1.000E-04  4.425E-02  O.OOOE+OO  9.557E-01 
1 . 000E-02  2 . 212E+03  5.000E-02  5.000E-04  4.425E-02  O.OOOE+OO  9.557E-01 
1 . 000E-02  4 • 424E+03  l.OOOE-Ol  1.000E-03  4.425E-02  O.OOOE+OO  9.55'E-01 
1 . 000E-02  8.847E+03  2.000E-01  2.000E-03  4.425E-02  O.OOOE+OO  9.557E-01 
1.000E-02  1 . 769E+04  4.000E-01  4.000E-03  4.425E-02  O.OOOE+OO  9.557E-01 
1 . OOOE-02  2 . 654E+04  6.000E-01  6.000E-03  4.425E-02  O.OOOE+OO  9.557E-01 
1 . 000E-02  4 . 424E+04  l.OOOE+OO  1.000E-02  4.425E-02  C.OOOE+OO  9.557E-01 
1. OOOE-02  2.21 2E+05  5.000E+00  5. 000E-02  4.425E-02  O.OOOE+OO  9.557E-01 
1.000E-02  4 . 424E+05  l.OOOE+Ol  l.OOOE-Ol  4.425E-02  O.OOOE+OO  9.3+J7E-01 
1. 000E-02  2.212E+06  5.000E+01  5.000E-01  4.425E-02  O.OOOE+OO  9.557E-01 
1.000E-02  4.424E+06  1.000E+02  l.OOOE+OO  4.425E-02  O.OOOE+OO  9.557E-01 


l.OOOE-Ol  4 . 439E-01  1.000E-05  1.000E-06  1.427E-02  8.236E-01  1.421E-01 
l.OOOE-Ol  2 . 220E+00  5.000E-05  5.000E-06  3.802E-02  3.795E-01  5.825E-OI 
l.OOOE-Ol  4.439E+00  1.000E-04  1.000E-05  4.350E-02  1.442E-01  3.123E-01 
l.OOOE-Ol  2 . 220E+01  5.000E-04  5.000E-05  4.442E-02  S.294E-03  9.5S5E-01 
l.OOOE-Ol  4.439E+01  1.000E-03  1.000E-04  4.442E-02  3.957E-09  P.556E-01 
l.OOOE-Ol  2 . 220E+02  5.000E-03  5.000E-04  4.442E-02  O.OOOE+OO  9.556E-01 
l.OOOE-Ol  4 . 439E+02  1. OOOE-02  1.000E-03  4.4a2E-02  O.OOOE+OO  9.556E-01 
l.OOOE-Ol  8 . 879E+02  2.000E-02  2.000E-03  4.442E-02  O.OOOE+OO  9.3E4E-01 
l.OOOE-Ol  1 . 776E+03  4.000E-02  4.000E-03  4.442E-02  O.OOOE+OO  9.55AE-01 
l.OOOE-Ol  2 . 664E+03  6. OOOE-02  6.000E-03  4.442E-02  O.OOOE+OO  9.356E-OI 
l.OOOE-Ol  4.439E+03  l.OOOE-Ol  1. OOOE-02  4.442E-02  O.OOOE+OO  9.53&E-01 
l.OOOE-Ol  2.220E+04  5.000E-01  5. OOOE-02  4.442E-02  O.OOOE+OO  9.55-4E-0I 
l.OOOE-Ol  4 .  439E  +  04  l.OOOE  +  OO  1.OO0E-01  4.442E-02  O.OOOE  +  OO  9.354E-'/: 
l.OOOE-Ol  2 . 220E+05  5.000E+00  3.000E-0 1  4.4-+2E-02  O.OOOE+OO  9.53*:£-ol 
l.OOOE-Ol  4 . 439E+05  l.OOOE+Ol  l.OOOE+OO  4.442E-02  O.OOOE+OO  9.3SoE-01 
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Table  C-5.  Optical  properties  in  the  visible  (A=0.55  pm)  of 
sooty  water  clouds  (continued). 

P  .  T  *  USOOt  UH20  R  T  £ 

H2n'  soot  (Optical  Depth)  (gm  cm-2)  (gm  cm-2 )  (Reflectivity) (Transmissivity) (Emissivity) 


5 . OOOE-Ol  9.013E-02  2.000E-Q6  1.000E-06  3.376E-03  9.618E-01  3.478E-02 
5 . OOOE-0 1  4 . 509E-01  1.000E-05  5.000E-06  1.454E-02  3.233E-01  1.621E-01 
5. OOOE-Ol  9 . 013E-01  2.000E-05  1.000E-05  2.439E-02  6.730E-01  2.976E-01 
5. OOOE-Ol  4 . 509E  +  00  1.000E-04  5.000E-05  4.425E-02  1.43SE-01  8.120E-01 
5. OOOE-Ol  9 . 0 1 8E  +  00  2.000E-04  1.000E-04  4.517E-02  2.071E-02  9.341E-01 
5. OOOE-Ol  4.509E+01  1.000E-03  5.000E-04  4.51SE-02  3.342E-09  9.543E-01 
5. OOOE-Ol  9.013E+01  2.000E-03  1.000E-03  4.518E-02  1.479E-17  9.543E-01 
5. OOOE-Ol  1 . 804E+02  4.000E-03  2.000E-03  4.518E-02  2.192E-34  9.548E-01 
5. OOOE-Ol  3.607E+02  8.000E-03  4.000E-03  4.513E-02  O.OOOE+OO  9.548E-01 
5. OOOE-Ol  5.411E+02  1.200E-02  6.000E-03  4.513E-02  O.OOOE+OO  9.548E-01 
5. OOOE-Ol  9 . 013E+02  2.000E-02  1.000E-02  4.518E-02  O.OOOE+OO  9.548E-01 
5. OOOE-Ol  4 . 509E+03  1. OOOE-Ol  5.000E-02  4.518E-02  O.OOOE+OO  9.548E-01 
5. OOOE-Ol  9.018E+03  2. OOOE-Ol  1. OOOE-Ol  4.518E-02  O.OOOE+OO  9.543E-01 
5. OOOE-Ol  4 . 509E+04  l.OOOE+OO  5. OOOE-Ol  4.518E-02  O.OOOE+OO  9.548E-01 
5. OOOE-Ol  9 . 018E+04  2.000E+00  l.OOOE+OO  4.518E-02  O.OOOE+OO  9.543E-01 


l.OOOE+OO  4 . 596E-02  l.OOOE-OA  1.000E-06  1.7&0E-03  9.307E-01  1.754E-02 
l.OOOE+OO  2 . 298E-01  5.000E-06  5.000E-06  8.154E-03  9.072E-01  8.469E-02 
l.OOOE+OO  4 • 596E-01  1.000E-05  1.000E-05  1.486E-02  8.230E-01  1.621E-01 
l.OOOE+OO  2 . 298E+00  5.000E-05  5.000E-05  3.952E-02  3.732E-01  5.823E-01 
l.OOOE+OO  4 . 594E+00  1.000E-04  1.000E-04  4.518E-02  1.432E-01  8.116E-01 
l.OOOE+OO  2 . 293E+01  5.000E-04  5.000E-04  4.413E-02  6.079E-05  9.533E-01 
l.OOOE  +  OO  4  <  596E  +  01  1.000E-03  1.000E-03  4.613E-02  3.703E-09  9.539E-01 
l.OOOE+OO  9 . 192E+01  2.000E-03  2.000E-03  4.613E-02  1.374E-17  9.539E-01 
l.OOOE+OO  1 . 838E+02  4.000E-03  4.000E-03  4.613E-02  1.392E-34  9.539E-01 
l.OOOE+OO  2 . 758E+02  6.000E-03  A.000E-03  4.A13E-02  O.OOOE+OO  9.539E-01 
l.OOOE+OO  4.59AE+02  1.000E-02  1.000E-02  4.613E-02  O.OOOE+OO  9.539E-01 
l.OOOE+OO  2 . 298E+03  5.000E-02  5.000E-02  4.613E-02  O.OOOE+OO  9.539E-01 
l.OOOE  +  OO  4 . 596E  +  03  1. OOOE-Ol  1. OOOE-Ol  4.613E-02  O.OOOE  +  OO  9.539E-01 
l.OOOE+OO  2 . 298E+04  5. OOOE-Ol  5. OOOE-Ol  4.A13E-02  O.OOOE+OO  9.539E-01 
l.OOOE+OO  4 . 596E+04  l.OOOE+OO  l.OOOE+OO  4.A13E-02  O.OOOE+OO  9.539E-01 


5 . OOOE+OO  1 . 058E-02  2.000E-07  1.000E-06  4.217E-04  9.960E-01  3.534E-03 
5 . OOOE+OO  5 . 292E-02  1.000E-06  5.000E-06  2.075E-03  9.304E-01  1.754E-02 
5. OOOE+OO  1 . 058E-01  2.000E-06  1.000E-05  4.070E-03  9.A11E-01  3.478E-02 
5. OOOE+OO  5 . 292E-0 1  1.000E-05  5.000E-05  1.748E-02  8.204E-01  1.621E-01 
5. OOOE+OO  1 . 058E+00  2.000E-05  1.000E-04  2.925E-02  6.733E-01  2.975E-01 
5. OOOE+OO  5.293E+00  1.000E-04  5.000E-04  5.252E-02  1.390E-01  3.085E-01 
5. OOOE+OO  1 . 059E+01  2.000E-04  1.000E-03  5.354E-02  1.933E-02  9.271E-01 
5. OOOE+OO  2.117E+01  4.000E-04  2.000E-03  5.356E-02  3.766E-04  9.461E-01 
5. OOOE+OO  4 . 234E  +  0 1  3.000E-04  4.000E-03  5.35AE-02  1.423E-07  O.  464E-01 
5. OOOE+OO  6.351E+01  1.200E-03  S.000E-03  5.356E-02  5.374E-U  9.464E-01 
5. OOOE+OO  1 . 0S8E+02  2.000E-03  1.000E-02  5.356E-02  7.667E-18  9.4&4E-0t 
5. OOOE+OO  S.293E+02  1.000E-02  5.000E-02  5.356E-02  O.OOOE+OO  9.464E-01 
5. OOOE+OO  1 . 059E+03  2.000E-02  1. OOOE-Ol  5.356E-02  O.OOOE+OO  9.444E-01 
5. OOOE+OO  5 . 293E  +  03  1. OOOE-Ol  5. OOOE-Ol  5.356E-02  O.OOOE  +  OO  9.4A4E-01 
5. OOOE+OO  1 . 059E  +  04  2. OOOE-Ol  l.OOOE  +  OO  5.356E-02  O.OOOE+OO  9.464E-01 


Table  C-5.  Optical  properties  in  the  visible  (A=0.55  ym)  of 
sooty  water  clouds  (continued). 

.  r*  Usoot  %0  R  T 

°H20/  soot  [Optical  Depth)  (gw  cm-2)  (gm  cm  )  (Reflectivity)  (Transmissivity)  (Emissivity) 


£ 

i 


l.OOOE+Ol  6.163E-03  1.000E-07  1.000E-06  2.514E-04  9,990E-0i  1.768E-03 
1 . OOOE+Ol  3.081E-02  5.000E-07  5.000E-06  1.247E-03  9.S99E-01  8.811E-03 
1 . OOOE+O 1  6.1A3E-02  1.000E-04  1.000E-05  2.468E-03  9.300E-01  1 . 754E-02 
1. OOOE+Ol  3 . 081E-01  5.000E-06  5.000E-05  1.140E-02  9.039E-01  8.468E-02 
1. OOOE+Ol  6.163E-01  1.000E-05  1.000E-04  2.072E-02  8.172E-01  1.621E-01 
1. OOOE+Ol  3 . 082E+00  5.000E-05  5.000E-04  5.415E-02  3.655E-01  5.803E-01 
1. OOOE+Ol  6 . 163E+00  1.000E-04  1.000E-03  6.140E-02  1.340E-01  8.046E-01 
1. OOOE+Ol  1 . 233E+01  2.000E-04  2.000E-03  A.251E-02  1.802E-02  9.195E-01 
1. OOOE+Ol  2.465E+01  4.000E-04  4.000E-03  6.253E-02  3.261E-04  9.371E-01 
1. OOOE+Ol  3.698E+01  6.000E-04  6.000E-03  6.253E-02  5.900E-06  9.375E-01 
1. OOOE+Ol  6.163E+01  1.000E-03  1.000E-02  A.253E-02  1.931E-09  9.375E-01 
1. OOOE+Ol  3.081E+02  5.000E-03  5.O00E-02  6.253E-02  O.OOOE+OO  9.375E-01 
1. OOOE+Ol  6 . 1 A3E+02  1.000E-02  l.OOOE-Ol  6.253E-02  O.OOOE+OO  9.375E-01 
l.OOOE+Ol  3 . 082E+03  5.000E-02  5.000E-01  A.253E-Q2  O.OOOE+OO  9.375E-01 
1. OOOE+Ol  6 . 163E+03  l.OOOE-Ol  l.OOOE+OO  6.253E-02  O.OOOE+OO  9.375E-01 


5. OOOE+Ol  2 . A25E-03  2.000E-08  1.000E-0&  1.145E-04  9.995E-01  3.539E-04 
5. OOOE+Ol  1.313E-02  1.000E-07  5.000E-06  5.714E-04  9.977E-01  1.769E-03 
5. OOOE+Ol  2 . A25E-02  2.000E-07  1.000E-05  1.140E-03  9.953E-01  3.534E-03 
5. OOOE+Ol  1.313E-01  1.000E-06  5.000E-05  5.595E-03  9.769E-01  1.754E-02 
5. OOOE+Ol  2 . 625E-01  2.000E-06  1.000E-04  1.093E-02  9.543E-01  3.478E-02 
5. OOOE+Ol  1.31 3E+00  1.000E-05  5.000E-04  4.577E-02  7.923E-01  1.620E-01 
5. OOOE+Ol  2 . 625E+00  2.000E-05  1.000E-03  7.455E-02  6.290E--01  2.964E-01 
5. OOOE+Ol  5 . 251E+00  4.000E-05  2.000E-03  1.042E-01  3.979E-01  4.979E-01 
5. OOOE+Ol  1 . 050E+01  8 . OOOE-OS  4.000E-03  1.209E-01  1.600E-C1  7.191E-01 
5. OOOE+Ol  1 . 575E+0I  1.200E-04  6.000E-03  1.236E-01  A.449E-02  9.119E-01 
5. OOOE+Ol  2 . 625E+01  2.000E-04  1.000E-02  1.241E-01  1.047E-02  8.654E-01 
5. OOOE+Ol  1.313E+02  1.000E-03  5.000E-02  1.241E-01  1.342E-10  8.759E-01 
5. OOOE+Ol  2 . 625E+02  2.000E-03  l.OOOE-Ol  1.241E-01  1.328E-20  8.759E-01 
5. OOOE+Ol  1.31 3E+03  1.000E-02  5.000E-01  1.241E-01  O.OOOE+OO  8.759E-01 
5. OOOE+Ol  2 . 625E+03  2.000E-02  l.OOOE+OO  1.241E-01  O.OOOE+OO  8.759E-01 


1 . 000E+02  2.183E-03  1.000E-08  1.000E-06  9.733E-05  9.997E-01  1.770E-04 
1 . 000E+02  1 . 092E-02  5.000E-08  5.000F-06  4.861E-04  9.986E-01  8.847E-04 
1 . OOOE+02  2 . 183E-02  1.000E-07  1.000E-05  9.709E-04  9.973E-01  t . 748E-03 
1 . 000E+02  1 .092E-01  5.000E-07  5.000E-05  4.302E-03  9.364E-01  3.9UE-03 
1. 000E+02  2.183E-01  1.000E-06  1.000E-04  9.473E-03  9.730E-01  I . 754E-02 
1. 000E+02  1 . 092E+00  5.000E-06  5.000E-04  4.260E-02  8.728E-01  8.464E-02 
1. OOOE+02  2 • 183E+00  1.000E-05  1.000E-03  7.511E-02  7.631E-01  l.ol9E-01 
1. 000E+02  4.366E+00  2.000E-05  2.000E-03  1.191E-01  5.856E-01  2.°53E-01 
1. 000E+02  3 . 733E+00  4.000E-05  4.000E-03  1.605E-01  3.479E-01  4.9li>E-01 
1. 000E+02  1 . 310E+01  6.000E-05  6.000E-03  1.752E-01  2.077E-01  6.171E-01 
1. OOOE  +  02  2 . 133E  +  01  1.000E-04  1.000E-02  1.82  3E-01  7.435E-02  7.433E-01 
1. 000E+02  1 . 092E+02  5.000E-04  5.000E-02  1.831E-01  2.602E-C6  3.166E-01 
1. 000E+02  2.183E+02  1.000E-03  l.OOOE-Ol  1.834E-01  7.006E-12  3.166E-01 
1. 000E+02  1.092E+03  3.000E-03  3.000E-01  1.334E-01  O.OOOE+OO  9.166E-01 
1. 000E+02  2 . 1 83E+03  1.000E-02  l.OOOE+OO  1.834E-01  O.OOOEtOO  8.166E-01 
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Table  C-5.  Optical  properties  in  the  visible  (x=0.55  um)  of 
sooty  water  clouds  (continued). 


aH2n/,pSOOt 

T  * 

fOotical  Depth) 

IJsoot 
(an  cm*  ) 

VO 

(gm  cm*  )  1 

R  T  e 

(Reflectivity) (Transmissivity) (Emissivity) 

2 . 000E+02 

1 .962E-03 

5.000E-09 

1 . 000E-06 

9.874E-05 

9.998E-01 

8.S51E-05 

2 . 000E+02 

9.B10E-03 

2 . 500E-08 

5 .  OOOE-OA 

4 . 434E-04 

9 .991E-01 

4.424E-04 

2 . 000E+02 

1 .9A2E-02 

5.000E-08 

1 .000E-05 

8 . 8A0E-04 

9.982E-01 

8.846E-04 

2 . 000E+02 

9 . 81 0E-02 

2 . 500E-07 

5 . OOOE-05 

4 . 399E-03 

9.912E-01 

4 . 4 1 5E-03 

2 . 000E+02 

1 .9A2E-01 

5 . OOOE-07 

1 .000E-04 

8 . 720E-03 

9.825E-01 

8.311E-03 

2 . 000E+02 

9.811E-01 

2 . 500E-06 

5 . 000E-04 

4.069E-02 

9. 160E-01 

4 . 327E-02 

2 . 000E+02 

1 .942E+00 

5 . 000E-06 

1. OOOE-03 

7 , 489E-02 

8.405E-01 

8.459E-02 

2 . 000E+02 

3.924E+00 

1 .000E-05 

2. OOOE-03 

1 . 281E-01 

7. 105E-01 

1 .614E-01 

2 . 000E+02 

7 . 848E+00 

2 . OOOE-05 

4. OOOE-03 

1 .938E-01 

5. 132E-01 

2.930E-01 

2 . 000E+02 

1.  177E+01 

3 . 000E-05 

b.  OOOE-03 

2 . 284E-01 

3.739E-01 

3 . 9?"?E-01 

2 . 000E+02 

1 .9A2E+01 

5 • 000E-05 

1 . 000E-02 

2.568E-01 

2.007E-01 

5.425E-01 

2 . OOOE+02 

9 . 8 1 1 E+O 1 

2 . 500E-04 

5.000E-02 

2.684E-01 

4.323E-04 

7 . 312E-01 

2 . 000E+02 

1 .962E+02 

5 . 000E-04 

1  .OOOE-Ol 

2.684E-01 

2.014E-07 

7.314E-01 

2. 000E+02 

9 . 310E  +  02 

2 . 500E-03 

5 . 000E-01 

2 . A84E-01 

4 . 468E-34 

7.31AE-01 

2.000E+02 

1 .9A2E+03 

5 .  OOOE-03 

1 .OOOE+OO 

2.684E-01 

0. OOOE  +  OO 

7 . 3 1 6E -0 1 

. 000E+O2 

1 .929E-03 

2 • 000E-09 

1 .000E-06 

8 . 357E-05 

9.999E-01 

3.541 E -05 

. 000E+02 

9. 147E-03 

1 .OOOE-08 

5 . OOOE-Ob 

4. 177E-04 

9. 994E-01 

1 . 7A9E-04 

.  000E+02 

1 .829E-02 

2.000E-08 

1 .000E-05 

8.349E-04 

9.988E-01 

3 . 539E-04 

. OOOE+02 

9. 147E-02 

1  .OOOE-07 

5. 000E-05 

4. 155E-03 

9.941E-01 

1 . 7A8E-03 

. 000E+02 

1 .829E-01 

2. 000E-07 

1 .OOOE-04 

8 . 261E-03 

9.882E-01 

3.534E-03 

. OOOE+02 

9.147E-01 

1 .000E-06 

5. OOOE-04 

3.943E-02 

9 . 430E-C1 

1  .  7  5  4  E  -  0  2 

. 000E+02 

1 .829E+00 

2 . 000E-06 

1 .000E-03 

7 . 454E-02 

8.907E-01 

3.476E-02 

. 000E+02 

3 . A59E+00 

4 . OOOE-OA 

2. OOOE-03 

1 . 340E-01 

7 . 978E-0 1 

6.822E-02 

. 000E+02 

7 , 318E+00 

8 . 000E-06 

4 , 000E-03 

2 . 209E-0 1 

6. 481E-01 

1 .31  IE-01 

. 000E+02 

1 .098E+01 

1 .200E-05 

6 . 000E-03 

2.739E-01 

5.328E-01 

1 .334E-0 1 

. 000E+02 

1 .829E  +  01 

2. 000E-05 

1 .000E-02 

3 . 457E-01 

3. A79E-01 

2.8A4E-01 

.OOOE+02 

9, 147E+01 

1  .OOOE-04 

5.000E-02 

4  ,  103E-01 

1 . 205E-C2 

5.776E-03 

. 000E+02 

1 .829E  +  02 

2 . 000E-04 

1 .OOOE-Ol 

4 . 104E-01 

1 . 74©E -04 

5.394E-01 

. 000E+02 

9. 147E  +  02 

1 .OOOE-03 

5 . OOOE-Ol 

4.  10'iE-01 

3.396E-19 

5 . 896E-01 

. OOOE+02 

1 .829E  +  03 

2 . OOOE-03 

1  .OOOE  +  OO 

4  .  104E-01 

0. OOOE  +  OO 

5 . 89AE-0 1 

1 .OOOE+03 

1 .785E-03 

1  .  OOOE-09 

1 .000E-O6 

8  .  134E-05 

9.999E-01 

1 . 77  2E-05 

1 . OOOE+03 

8 . 92AE-03 

5 . 000E-09 

5 . 000E-06 

4 . 09 1 E -04 

9 .995E-0  l 

8.947E-05 

1 .000E+03 

t . 785E-02 

1 . 000E-08 

1 . OOOE-05 

8.  179E-04 

9.990E-01 

1  . 77CE-04 

1 .000E+03 

3 . 926E-02 

5 . 000E-08 

5. 000E-05 

4 . 073E-03 

9.950E-01 

8 . 847E-04 

1 .000E+03 

1 .785E-01 

1  .OOOE-07 

1  .  000E-04 

3. 106E-03 

9.901E-C1 

1  .  •’  6  8  E  -  0  2 

1 .000E+03 

8.926E-01 

5. 000E-07 

5. OOOE-04 

3.898E-02 

9.522E-01 

8.91  IE-03 

1 .OOOE+03 

1 . 785E  +  00 

1  .  000E-O6 

1 .000E-03 

7. 43SE-02 

9 . 08 1 E -  0  1 

1 . 7  4  E  -  0  2 

1 .OOOE+03 

3 . 570E+00 

2 . 000E-06 

2. 000E-03 

1 .361E-01 

8.292E-01 

3 . 4‘75E-02 

1 .000E+03 

7 . 1 41E+00 

4 .000E-06 

4 . 000E-03 

2 . 31 1E-0 1 

7.005E-01 

6.310E-02 

1 .000E+03 

1 .071E  +  01 

6 . 000E-06 

6.  OOOE-03 

3 . 003E-0 1 

5.993E-01 

9 . 99  IE -02 

1 .000E+03 

1 .735E+01 

1 .000E-05 

1  .  000E-02 

3 . 898E-01 

4 . 515E-01 

1 . 58  7E-01 

1 . 000E+03 

8.926E+01 

5. 000E-05 

5 . 000E-02 

5.227E-01 

4.270E-02 

4 . 346E-01 

1 .OOOE+03 

1 .7S5E  +  02 

1 . OOOE-04 

1 .OOOE-Ol 

5 . 240E-01 

2.509E-03 

4  .  9 35E - 0 l 

1 .000E+03 

9.926E+02 

5 . 000E-04 

5. OOOE-Ol 

5 . 240E-01 

3.5S8E-13 

4 . -JOE-Ol 

1 .000E+03 

1 . 735E+03 

1  .  OOOE-03 

1 . OOOE+OO 

5. 240E-01 

1  .  775E-2S 

'  J 

l> 

O 

rn 

i 

o 
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Table  C-5.  Optical  properties  in  the  visible  (A=0.55  um)  of 
sooty  water  clouds  (continued). 


a  / o  T  soot 

^2n  soot  'Hoti'cal  lepth)  'am  cm-  ) 


jH20  3  T  e 

'  qm  cm'2)  f°ef leetivity)  .'Transmissivity)  fEmissivity) 


5 . 000E+03 
5 . 000E+03 
5 . 000E+03 
5 . 000E+03 
5 . 000E+03 
5 . 000E  +  03 
5 . OOOE  +  03 
5 . 000E  +  03 
j . 000E+03 
5 . 000E  +  03 
5 . 000E+03 
5 . 000E+03 
5 . 000E+03 
5 . 000E+03 
5 . 00OE  +  03 


1 .750E-03 
8 . 749E-03 
1 .750E-02 
8.749E-02 
1 .750E-01 
8.749E-01 
1 .7S0E  +  00 
3 . 500E+00 
6 . 999E+00 
1 . 050E+0  1 
1 .750E+01 
8 . 749E+0 1 
1 .750E+02 
9 . 7  49E  +  02 
1 . 750E  +  03 


2.000E- 
1 .000E- 
2.000E- 
1 .000E- 
2.000E- 
1 .OOOE- 
2.000E- 
4.000E- 
9 . OOOE- 
1 .200E- 
2.000E- 
1 . OOOE- 
2 . OOOE- 
1 .OOOE- 
2.000E- 


.  .OOOE-06 
.OOOE-06 
.  OOOE-OS 
.  OOOE-05 
.  OOOE-04 
.  000E-04 
.  000E-03 
.  000E-03 
1.000E-03 
> . 000E-O3 
.  .000E-O2 
;.000E-02 
.  .OOOE-Ol 
i . OOOE-Ol 
.  . OOOE+OO 


3.051E- 
4.023E- 
8 . 042E- 
4 . 008E- 
7.932E- 
3.S62E- 
7.424E- 
1 .377£- 
2.404E- 
3 . 196E- 
4 . 328E- 
7. 129E- 
7.415E- 
7.441E- 
7.441E- 


►.999E- 
>.996E- 
r ,  992E- 
’  .958E- 
».917E- 
f . 596E- 
’ . 222E- 
1.552E- 
’  .455E- 
>.594E- 
i . 325E- 
.  .405E- 
I.012E- 
!.564E- 
,  .  473E- 


I.S11E-06 
. 769E-05 
I.543E-05 
.  .  771E-04 
I.539E-04 
.  769E-03 
524E-03 
’ .  054E-03 
. 405E-02 
! . 098E-02 
: .  462E-02 
. 467E-0J 
t .  184E-01 
I.558E-01 
I.559E-C1 


1 .000E+04 
1 .000E+04 
1 .000E+04 
1 . 000E  +  04 
1 .000E+04 
1 .000E+04 
1 .000E+04 
1 .000E  +  04 
t . 000E+04 
1 . 000E  +  04 
1 . 0008+04 
1 .000E+04 
1 . 0008+04 
1 . 000E+04 
1 . 000E+04 


1 .715E-03 
8 . 727E-03 
1 . 745E-02 
8. 727E-02 
1 . 745E-01 
8. 727E-01 
l . 745E+00 
3.491E+00 
6 . 982E+00 
1 . 047E+01 
1 . 745E+01 
8. 727E+01 
1 .745E+02 
8 . 727E  +  02 
1  .  745E+03 


1 .OOOE- 
5. OOOE- 
1 .OOOE- 
5.000E- 
1 .OOOE- 
5.000E- 
1 . OOOE - 
2  .  OOO'. 
4.000E- 
6.000E- 
1 .OOOE- 
S.OOOE- 
1 . OOOE - 
S.OOOE- 
1 .OOOE- 


1 .000E-06 
5. OOOE-06 
l  . OOOE-05 
5. OOOE-OS 
1 .000E-04 
5.  OOOE-04 

1  .  OOOE -03 
! .  000E-03 
I.000E-03 
> .  0OOE-O3 
L  .  000E-02 
5 .  OOOE-02 
t  .OOOE-Ol 
5. OOOE-Ol 
1  .OOOE  +  OO 


8.032E- 
4.014E- 
8.026E- 
3.999E- 
7 . 9.J4E - 
3.858E- 
7.422E- 
1 . 379E- 
2 . 416E- 
3. 221E- 
4 . 388E- 
7.530E- 
8.014E- 
8  .  10<»E- 
8. 109E- 


.999E-01 
.996E-01 
.992E-01 
.959E-01 
.919E-01 
.605E-01 
.240E-01 
. 585E-0 1 
.514E-01 
.  6  7  3  E  -  0  1 
.  4  3  7E-01 
. ob 3E -0 1 
. 425E-02 
.  122E-05 
. 48 1 E -08 


1 .821E-0.*. 
3.817E-06 
1  .■’63E-05 
3 . 953E-05 
1 . 770E-04 
3.848E-04 
1 . 769E-03 
3.534E-03 
7 .  0S2E-O3 
1 .056E-02 
1 . .'51C-0 2 
3.021E-02 
l . 3448-01 
1 . 39 OE-01 
1  .  3  9  l  F  -  O  I 


1 . OOOE+OS  1. 
1 . OOOE+OS  8. 
1 . OOOE+OS  1. 
1  .  OOOE  +  OS  8. 
1 . 000E  +  05  1. 
1 . OOOE  +05  3. 
1 . 000E+05  1. 
1 . OOOE+OS  3. 
I . OOOE+OS  6. 
1 . OOOE+OS  1, 
1 . OOOE+OS  1 , 
1  . OOOE+OS  3. 
1.000E+05  1, 
1. OOOE  +  OS  9. 
1. OOOE  +  OS  1. 


74 1 E -03 
707E-03 
74 1 E -02 
707E-02 
74  IE-01 
707E-01 
741E+00 
433E+00 
966E+00 
045E+01 
741E+01 
707E+01 
741E+02 
707E+02 
’4 1E+03 


1  .  OOOE- 1 1 
5 . OOOE - 1 1 
1 .00OE-10 
5 • OOOE-l 0 
1 . OOOE -09 
S-OOOE-09 
1 . 000E-08 
2 . OOOE -08 
4 . 000E-08 
6 . OOOE -08 
1 . 000E-07 
S  .  OOOE-O7 
1 . OOOE -  06 
S  .0008-06 
1 . OOOE-OS 


1 .OOOE-06 
5 . OOOE -06 
1 .OOOE-OS 
5 . OOOE-OS 
1 .OOOE-04 
5. OOOE-04 
1 . 000E-03 
2 . 000E-03 
4 . 000E-03 
6.000E-03 
1 . 000E-O2 
S . OOOE-02 
1 . OOOE-Ol 
S . OOOE-Ol 
1 . OOOE+OO 


8.017E-05 
4 . 00  7E -0  4 
9 . 01  IE-04 
3 . 992E-03 
7.9S2E-03 
3. 853E-02 
7 . 420E-02 
1 . 381E-01 
2 . 426E-C1 
3 . 24SE-01 
4  .  4  4  3  E  -  0  1 
7.952E-0  1 
9. 783E-01 
9.251F-01 
9. 357E-0  1 


o.99ac. 
9 . 9-?6E 
9 . 992E - 
9 . 940E - 
9 . 920E- 
9.614E- 
9.256E- 
3.615E- 
7  .  S  6  6  E  - 
6  .  7  4  5  E  - 
5  .S39E- 
1  .9.S1E- 
l .04SE- 
B  .  6  8^8- 
6 . 008E- 


01  1.S20E-07 
01  8.3S'»E-07 
01  1.327E-06 
01  9 . 3  \  -  0  6 
01  I.''7 'E-OS 
01  9.366E-0S 
01  1 . 7  ~  3  E -  0  4 
01  2.5448-01 

01  ■’.0888-04 
01  1  .0638-0  ,.i 
0  1  1  .  7  7 1 E - 0  3 
0  1  2  .  r 2E  -0  3 
01  1 . 7 1 7  E - 0  2 
0  7  S..120E  -02 
04  .5.  373E-02 


Table  C-5.  Optical  properties  in  the  visible  ( x  =0 . 5 5  um)  of 
sooty  water  clouds  (continued). 


/q  n2u  ,  iout 

Hjfr  soot  f Opt i cal  Oepth)  (gm  cm*2)  (qm  cm*2)  (Reflectivity) (Transmissivity)  (Emissivity) 


1. OOOE+OA  1.741E-03  1.000E-12  l.OOOE-OA  8.024E-05  9.999E-01  O.OOOE+OO 
1.000E+06  8 . 705E-03  5.000E-12  5.000E-06  4.011E-04  9.996E-01  5.971E-08 
1 • OOOE+OA  1 . 741E-02  1.000E-U  1.000E-05  8.001E-04  9.992E-01  1.790E-07 
1 . OOOE  +  OA  8.705E-02  5.000E-U  5.000E-05  3.991E-03  9.9A0E-01  8.924E-07 
t . 0OOE+O6  1.741E-01  1 . OOOE-IO  1.000E-04  7.951E-03  9.921E-01  1.778E-06 
1. OOOE+OA  8 . 705E-01  5.000E-10  5.000E-04  3.853E-02  9.A15E-01  3.96AE-06 
1. OOOE+OA  1.741E+00  1.000E-09  1.000E-03  7.420E-02  9.253E-01  1.800E-05 
1. OOOE+OA  3 . 482E+00  2.000E-09  2.000E-03  1.381E-01  8.A19E-01  3.590E-05 
1. OOOE+OA  A . 9A4E+00  4.000E-09  4.000E-03  2.427E-01  7.572E-01  7.131E-05 
1. OOOE+OA  1 . 045E+01  A.0O0E-09  A.OO0E-O3  3.247E-01  6.752E-01  1.077E-04 
X. OOOE+OA  1 . 741E+01  1.000E-08  1.000E-02  4.448E-01  5.550E-01  X.796E-04 
1. OOOE+OA  8 . 705E+01  5.000E-08  5.000E-02  7.998E-01  1.993E-01  8.969E-04 
1. OOOE+OA  1.741E+02  1.000E-07  1.000E-01  8.880E-01  1.103E-01  1.790E-03 
1. OOOE  +  OA  8 « 705E+02  5.000E-07  3.000E-01  9.700E-01  2.153E-02  8.440E-03 
1. OOOE+OA  1.741E+03  l.OOOE-OA  l.OOOE+OO  9.77AE-01  7.850E-03  1.450E-02 


1 . 000E+07  1.741E-03  1.000E-13  l.OOOE-OA  7.961E-03  9.999E-01  8.889E-08 
1 . 000E+07  8.705E-03  S.OOOE-13  3.000E-0A  4.021E-04  9.99AE-01  O.OOOE+OO 
1.000E+07  1 . 741E-02  l.OOOE-12  1.000E-05  7.998E-04  9.992E-01  O.OOOE+OO 
1 . OOOE+07  8 . 705E-02  5.000E-12  5.000E-0S  3.990E-03  9.9A0E-01  3.354E-0S 
1 . 000E+07  X.741E-0X  l.OOOE-11  1.000E-04  7.949E-03  9.921E-01  2.A46E-07 
X. 000E+07  8 . 705E-0 1  5.000E-11  5.000E-04  3.853E-02  9.A15E-01  1.11XE-0A 
X. 000E+07  1.741E+00  1.000E-10  1.000E-03  7.420E-02  9.258E-01  2.059E-06 
1. OOOE+07  3 . 482E+00  2.000E-10  2.000E-03  1.332E-01  8.A13E-01  4.142E-0A 
X. 000E+07  A . 9A4E+00  4.000E-10  4.000E-03  2.428E-01  7.572E-01  8.096E-06 
1. 000E+07  1.045E+01  6.000E-10  6.000E-03  3.247E-0X  6.753E-01  1.217E-03 
1. 000E+07  1.74XE+0X  1.000E-09  1.000E-02  4.449E-0X  5.551E-01  2.025E-05 
1. 000E+07  8 . 703E+01  5.000E-09  5.000E-02  3.002E-01  X.997E-01  1.013E-04 
1. 000E+07  1.741E+02  1.000E-08  1.000E-01  8.389E-01  1.109E-01  2.02AE-04 
1. 000E+07  8 . 705E+02  5.000E-08  3.000E-01  9.750E-01  2.400E-02  1.00AE-03 
1. 000E+07  1 . 741E+03  1.000E-07  l.OOOE+OO  9.864E-01  1.1A6E-02  1.972E-03 


1 . 000E+08  1 . 741E-03  1.000E-14  l.OOOE-OA  8.151E-05  9 . 99<»E-01 -9 . 60 7E-0U 
1 . 000E+08  8 . 705E-03  5.000E-14  5.000E-06  3.984E-04  9.99AE-01  9.&04E-03 
1 . 000E+08  1.74XE-02  l.OOOE-13  1.000E-05  7.965E-04  9.99’E-Ol  O.OOOE+OO 
1 . 000E  +  08  8 . 705E-02  5.000E-13  5.000E-05  3.983E.-03  9.960E-01  O.OOOE  +  OO 
1 . 000E+O8  1.741E-01  l.OOOE-12  1.000E-04  7.953E-03  9.920E-01  O.OOOE+OO 
1.000E+08  8 . 703E-01  3.000E-12  3.000E-04  3.833E-02  9.615E-01  1.348E-07 
1 . 000E  +  08  1  .741 E  +00  l.OOOE-ll  1.000E-03  7.420E-02  9.258E-01  5.333E-0" 
1 . 000E  +  08  3 . 482E  +  00  2.000E-11  2.000E-03  1.3B1E-01  8.613E-01  9.U3E-07 
1 . 000E+08  6.9A4E+00  4.000E-11  4.000E-03  2.428E-01  7.572E-01  1.748E-0A 
1 . 000E+08  1.043E+01  6.000E-11  6.000E-03  3.247E-01  6.753E-01  2.599E-0A 
1 . 000E+08  1.741E+01  1.000E-10  1.000E-02  4.449E-01  5.551E-C-I  4.385E-0A 
1 . 000E+08  3.705E+01  3.000E-10  5.000E-02  8.003E-01  1.997E-01  2.17OE-05 
1 . 000E+08  1 . 741E+02  1.000E-09  l.OOOE-Ol  8.890E-01  1.109E-01  4.338E-03 
1.000E+08  8 . 705E+02  3.000E-09  3.000E-01  9.755E-01  2.427E-02  2.1AAE-04 
1 . 000E+08  1 . 741E+03  1.000E-08  l.OOOE+OO  9.374E-01  1.21SE-02  4.312E-04 


Table  C-5. 

Optical 

properties 

in  the  visible  (A=0. 

55  urn)  of 

sooty  water  clouds 

(Concluded) . 

★ 

Usoot 

Uh20  . 

R 

T 

£ 

h 2 n  soot  (Ootical  OeotM 

( am  cur* ) 

(qm  cnr'M  (Ref lectivity)  (Transmissivity )  (fmissivity) 

1 .000E+09 

1 .741E-03 

1 .000E-15 

1 .000E-0A 

7.970E-05 

9 .999E-01 

6 • 038E-08 

1 .000E+09 

8.705E-03 

5 . 000E-15 

5 . OOOE-OA 

3 . 984E-04 

9.99AE-01 

A.03AE-03 

1 .00OE+09 

1 .741E-02 

1 .OOOE-14 

1 . 000E-05 

7.9A4E-04 

9 .992E-01 

A.034E-08 

1 .000E+09 

8 . 705E-02 

5.000E-14 

3. OOOE-OS 

3.992E-03 

9 .9A0E-01 

A.015E-03 

1 .000E+09 

1 .741E-01 

1  .OOOE-13 

1 .OOOE-04 

7 . 953E-03 

9.920E-01 

0. OOOE+OO 

1 .000E+09 

3.70SE-01 

5. OOOE-13 

5 . 000E-04 

3 . 853E-02 

9.615E-01 

1 .742E-07 

1 .000E+09 

1 .741E+00 

1.000E-12 

1 .000E-03 

7. 420E-02 

9 . 258E-01 

2.79AE-0? 

1 .0O0E+09 

3 . 482E+00 

2.000E-12 

2 .  OOOE-03 

1 . 331E-01 

8.618E-01 

5.727E-0? 

1 .000E+09 

A . 944E+00 

4. 000E-12 

4 • 000E-03 

2.428E-01 

7 . 572E-01 

l  .098E-06 

1 .000E+09 

1 .045E+01 

A.000E-12 

A . 000E-03 

3 . 247E-01 

6.753E-01 

1 . A33E-06 

1 .000E+09 

1 .741E+01 

l.OOOE-11 

1 .000E-02 

4.449E-01 

5.5S1E-01 

2.755E-06 

1 .000E+09 

8.705E+01 

5.000E-11 

5.000E-02 

8 . 003E-0 1 

1.997E-01 

1  » 373E-05 

1 .OOOE+09 

1 .741E+02 

1  .OOOE-IO 

1 .OOOE-Ol 

8. 891E-01 

1 .  109E-C1 

2.745E-05 

1 .000E+09 

8 . 705E+02 

5 .  OOOE-IO 

5.000E-01 

9 . 75AE-01 

2.430E-02 

1 .371E-04 

1 .000E+09 

1.741E+03 

1 .000E-09 

1 .OOOE+OO 

9 . 875E-01 

1 . 223E-02 

2 • 735E-04 

1 .000E+10 

1 .741E-03 

1 .000E-1A 

1.000E-06 

3.212E-05 

9.999E-01 

0. OOOE  +  OO 

X .OOOE+IO 

8.705E-03 

5.000E-1A 

5 . 000E-06 

3 . 987E-04 

9 . 99AE-0 1 

0. OOOE  +  OO 

1 .OOOE+XO 

1 .741E-02 

1  .OOOE-13 

1 .000E-05 

7 . 972E-04 

9.992E-01 

A .  21 7E -08 

1 .OOOE+IO 

8.705E-02 

5 . 000E-15 

5  •  000E-05 

3 . 997E-03 

9.960E-01 

A. 197E-08 

X .OOOE+XO 

1 . 741E-01 

1 .OOOE-14 

1 .000E-04 

7 . 950E-03 

9.921E-01 

1 .235E-07 

1 .OOOE+IO 

8.705E-01 

5. OOOE-14 

5. 000E-04 

3 . 853E-02 

9.A15E- 01 

1 . 795E-07 

1 .OOOE+IO 

1 . 741E+00 

1 .OOOE-13 

1  .OOOE-03 

7 , 420E-02 

9.238E-01 

2.305E-9? 

1 .OOOE+IO 

3  .  482E  +  00 

2. OOOE-13 

2. OOOE-03 

1 .3B2E-01 

8.  A19E-01 

3.3A3E-07 

1 .OOOE+IO 

A.9A4E+00 

4. OOOE-13 

4. OOOE-03 

2 . 429E-01 

7.572E-01 

1 . 095E-OA 

1 .OOOE+IO 

1 .045E+01 

A. OOOE-13 

6. 000E-03 

3 . 247E-01 

A .  753E-01 

1 .599E-06 

1 .OOOE+IO 

1 .74XE+01 

1 rOOOE-12 

1 .OOOE-02 

4.449E-01 

5.531E-01 

2.59AE-0A 

1 .OOOE+XO 

8 . 705E  +  0 1 

5.000E-12 

5  .  OOOE-02 

8 . 003E-01 

1 .997E-01 

1 .294E-05 

1 .OOOE+IO 

1 .741E+02 

X  .OOOE-ll 

1 .OOOE-Ol 

3 . 89  IE -01 

1 . 109C-01 

2  •  53aE-0!j 

1 .OOOE+IO 

8 . 705E+02 

5.000E-11 

5. OOOE-Ol 

9. 73AE-01 

2 . 430E-02 

1 . 292E-04 

1 .OOOE+IO 

1 . 74XE+03 

1  .OOOF.-IO 

1  .OOOE  +  OO 

9 .873E-01 

1 . 224E-02 

2.57AE-04 

*  Note  that  the  optical  depth  is  defined  by  Equation  2-36  and  includes 
effects  due  to  both  scattering  and  absorption. 
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APPENDIX  D 

SOLUTIONS  TO  THE  NUCLEAR  WINTER  MODELS 


THE  PROBLEM  OF  FIGURE  3-1 


The  solution  to  the  problem  of  Figure  3-1  is  as  follows: 


I  n  1  ayer  1 : 


S  =  (1-Aj)  x 


U  =  \  ( 1  -Aj )  (1-A2)  l  (2  +  £j  2j) 
D  =  1  (1-Ai)  (1-A2)  £ 


x  --  L  (1-Aj)  ( 1  -  A2  )  Z  (1*6^2!  ) 


where  z,  measures  the  distance  from  the  top  of  layer  1. 


In  1 ayer  ? : 

-e  z2 

S  =  (1-Aj)  (1-A2)  l  e  52 


( D- 1 ) 


M  =  i  (1-Ai)  (1  -A2 )  Z  [lecx  +  uj  f(l-u)  e'ts. 


^  ^  ( 1  -  A  ■  )  (1-A2)  l  Ui  +  ii  |  + 

0  1  l 


(D-2) 


l  ( i _a -  ?  I -A  )  :  ♦  ( L  - 1 


wh'-re  /  is  th»*  listari  .r  r*M,i  tne  r  >o  of  lriy>>r  ?. 


PREVIOUS  MOI 
IS  BLANK 


In  layer  3: 


S  =  (1-AJ  (1-A2)  e  e  “s2 

U  =  L  ( 1 -Ax )  (1-A2)  E  [l+a+u^  + ( 1 -a )  e  Us2  +  e  Ms2  £j  z3j 

( D-3 ) 

D  =  L  (1-Ai )  (1-A2)  E  [1  +  a  +  Uj  -(1+a)  e  US2  +  e  Us2  Ej  z3] 

t  =  L  ( 1  -  A  x )  (1-A2)  e  [l  +  +  a(l-e  U$2)  +  e  P$2  £j  z3J 

where  z3  is  the  distance  from  the  top  of  layer  3. 

For  the  earth: 

x  =  L  (1-a3 )  (1-A2)  e  [l+u  +  (l+p  )e'Us2  +  a(l-e’Ws2)]  (D-4) 

2  1 1  1 3 

THE  PROBLEM  OF  FIGURE  3-4 

An  approximate  solution  to  the  problem  of  Figure  3-4  in  the  case 
Uj^O  and  defining  A  =  ( 1  -Aj )  (1-A2)  is  as  follows: 

In  the  convection  zone: 

-e  z2 
S  =  A  E  e  S2 

U  =  o  [(T0  +  +  4  JL-  (T0  +  6  z2 ) 3  +  12  ( _§_)2  (T0  +  6z2)2 

CI2  eI2 

(D-5 ) 

+  24  (JL)3  (T0  +  Bz2)  +  24  -  "i2-  E  A  e’^ 

'iz  I2  I2 
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0  =  0  [  (T0  +  Bz2)4  -  4  {T0  +  Sz2)3  +  12  (_§_)  (T„  +  Bz2)2 


ei2  £I2 


\'fP\ 


-  24  (_§_)  (T0  +  Bz2)  +  24  (J_)  ]  -  [To" 

£  *  £  t 


4  (J_)  T03  .  ;2  (-i-)2T02  -  24  (-§-)3T0  +  24  (-i_)4]  e"^*2 


CI2  CI2 


& 


T  =  T0  +  6  z2 


where  z2  is  the  distance  from  the  top  of  the  convection  zone. 


In  the  convection- free  region  of  layer  2 


S  =  A  E  e 


"ec„  z2 


-e  d  -e  z2 

11  =  1  [AEe  2  |  ( 1-a ) e  2  +  o  -  1}  +  2o(T0  +  Bd)4J 

2  a 


-es  d  _es  z2 

D  =  L  [AEe  2  [-(l+a)e  2  +  a-i}+2a(T0+6d)4] 

2  a 


-e,  d  -e  Z2 

t  =  I  [AEe  2  (a  -  I)  (1-e  2  )  +  2a(T0  +  Bd)4] 


(0-6) 


where  z2  is  the  distance  from  the  top  of  the  convection-free  region  of 


S  =  A 


U  =  I 
2 


D  =  L 
2 


t  =  I 
2 


where 


r  =  I 

2 


In  1 ayer  3 : 


£  e 


[2o(T0  +  iid)1*  +  An  (a  +  1)  e  ^  -  a  e  Us2  +  e  Us2(l  +  e,  z3)}j 

a  1 3 


% 

t'r*2 


[2o(To  +  yd)4  +  Ai.{(a  +  L)  e  S2  -  a  e  S2  +  e  S2(-l  +  z3)jj 


[2a(T0  f  yd)4  +  Az{ (a  -  1)  e  a  e  ^  ♦  e  ^  e.  z3}J 

a  1 3 


z3  is  the  distance  from  the  top  of  layer  3. 


(D-7) 


>>5 

& 

Sl 

Wr 


>1 

;5>>l 


For  the  earth 


L 2 o ( T o  +  yd)4  +  Aii{  (a  -  L)  e  S2  -  a  e  S2  +  e  S2  (1  +  y  ) }  j  (D-8) 

a  1 3 


In  the  above  equations: 


THE  CASE  OF  HIGH  SURFACE  ALBEDO 


For  the  solar  fluxes  in  the  case  where  the  earth  has  an  albedo, 
Ag,  in  the  solar  we  have: 

Above  layer  1: 


S 


D 


=  E 


Sy  =  S  [Aj.  +  A2 


AjA2  + 


(l~Ai )2 (1-A2 )2  Ap  e'2"52 

-2us 

1  -  A2  Ae  -  A!  Ae(l-A2)  e  S2 


(D-9) 


In  layer  1: 

SQ  =  ( 1  -Aj )  I 


(l-AO 


+  (1-Aj)2(l-A2)2  Ag  e  ^ 

-2u 

1  -  A2  Ag  -  Aj  Ae(l-A2)  e 


(D-10) 


In  layer  2: 


2m, 


sn  .  !  Id-*.)  (1-A;)  *  A,  (1-A,);(1-A;);  Ae  e 


e  z. 


-2m, 


-j  e  s2 


1  .  A2  Ag  -  Aj  Ae(l-A2)  e  s2 

c  _  r  r  (1-Ai)(l-A2)2  a  e  S2  e  Z2 

‘  1  L - “ - j  e  s2 

U  -2uc 

1  -  A2  Ae  -  Aj  Ae(l-A2)  e  S2 

where  z2  is  the  distance  from  the  top  of  layer  2. 


( D- 11) 
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APPENDIX  E 

CONVECTIVE  EQUILIBRIUM  WITH  MULTIPLE  IR  GROUPS 


In  this  appendix  we  shall  demonstrate  that  a  region  as  might  be 
formed  by  the  NRC  baseline,  soot  and  water-vapor  layer  (Reference  E-l)  can 
be  represented  by  a  set  of  IR  groups  with  constant  opacities  such  that  the 
resulting  temperature  profile  is  significantly  different  from  that  for  a 
single  IR  group.  In  particular,  if  the  range  of  IR  opacities  needed  to 
represent  the  water-vapor  emissivity  includes  opacities  both  less  than  and 
greater  than  the  sol ar  opacity,  the  tenperature  profile  will  always  have  a 
peak  at  about  one  solar  optical  depth.  The  temperature  gradient,  in  the 
absence  of  convection,  on  the  top  of  the  layer  will  normally  exceed  the 
critical  lapse  rate. 

Without  using  a  dynamical  model  for  the  profile  changes  which 
result  from  convection,  it  is  difficult  to  be  confident  that  a  complete 
understanding  has  been  gained.  However,  it  is  possible  to  include  a  con¬ 
vective  flux  component  with  transfer  equations  in  order  to  find  a  new 
equilibrium  condition  which  has  a  convective  layer  which  holds  the  tem¬ 
perature  gradient  at  the  critical  lapse  rate  through  a  modest  mixing  of 
the  unstable  region  with  the  optically  inert  atmosphere  above  it.  This 
results  in  a  modified  temperature  profile  including  an  altered  asymptotic 
temperature  at  the  base  of  the  soot  layer. 

OPTICAL  PROPERTIES  OF  THE  LAYER 

The  NRC  baseline,  30°-70°  N  case  shows  that  after  ?5  days  a 
minimum  surface  temperature  is  reached  and  at  this  time  the  optical  depth 
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for  solar  extinction  is  2.5  due  to  the  presence  of  soot  which  is  taken  to 
be  uniformly  distributed  from  0  to  9  km  above  the  surface.  The  optical 
depth  is  the  product  of  the  physical  depth,  z,  and  the  solar  opacity,  ^  , 
so  we  can  take  xs  “  0.25  km-1.  The  broad-band  infrared  emissivity,  7,  is 
mainly  from  water  vapor  and  its  behavior  with  path  length,  U,  is  given  by 
Figure  E-l. 


Figure  E-l.  Broad-band  IR  emissivity  for  water  vapor 
(from  Reference  2-10). 


In  the  absence  of  scattering  and  for  a  parallel  beam,  the  emis¬ 
sivity  is  related  to  the  path  length  by 


=  1 


-kU 


(E-l) 


m 


where  k  is  a  mean  mass  absorption  coefficient  with  units  of  cm2  g-1.  Here 
kU  is  the  optical  depth. 

The  NRC  baseline  case  provides  for  a  water-vapor  pathlength  of 
1.45  gm  cm-2  giving  an  emissivity  of  0.63  from  Figure  E-l.  Thus,  the 
optical  depth  for  water  vapor  will  be  very  nearly  one  for  a  mean  opacity 
over  a  9  km  depth  of  x^q  *  0.1  km-1. 

Use  of  this  mean  value  will  give  the  correct  emissivity  at 
U  =  1.45  gm/cm2  but  it  would  provide  a  very  poor  representation  for  the 

F(U)  curve  of  Figure  E-l.  In  fact,  F(U)  results  from  an  averaging  over  a 

complex  absorption  structure.  It  is  possible  to  represent  this  curve  as 
resulting  from  a  set  of  IR  groups  with  signif icantly  different  optical 

properties.  For  example,  a  reasonable  fit  to  e*(U)  is  given  by 

e  -  .2(1  -  e'1300,J)  +  .2(1  -  '50U) 

(E-2) 

+  .31.1  -  e  J  +  .3(1  -  e  ) 

over  the  range  of  U  from  10"1*  to  102.  This  form  for  F  would  result  if 

there  were  a  set  of  four  distinct  IR  frequency  bands  which  contributed 

with  uniform  opacities  within  each  band.  The  opacities  x  =  kp  of  the  four 
groups  in  this  fit  correspond  to,  for  this  mean  density  if  p  *10"6 
gm/cm3,  values  of  130  km*1,  5  km-1,  0.1  km"1,  and  0.012  km"1.  Recall  that 

the  solar  opacity  was  x$  “  0.1  km"1  so  that  there  are  IR  groups  with  sig¬ 

nificant  activities  which  are  both  more  and  less  active  than  the  solar 
absorption. 

TRANSFER  EQUATIONS 

The  equation  for  the  steady  state  radiative  transfer  for  a 
planar  atmosphere  can  be  written  as  (Reference  E-?) 
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where  the  specific  intensity  I(z,n,v)  is  a  function  of  position  z,  direc¬ 
tion  and  frequency  and  has  units  of  ergs  cm-2  sec"1  Hz"1  sr"1.  Here  y  is 
the  cosine  of  the  angle  relative  to  vertical  while  n  is  the  emission  co¬ 
efficient  and  x  is  the  opacity  with  units  of  cm"1.  Generally  x  is  the  sum 
of  absorption  and  scattering  contributions  but,  for  the  present  discus¬ 
sion,  scattering  will  be  ignored.  Given  n  and  x»  both  possibly  functions 
of  frequency,  direction,  and  position,  it  is  possible  to  solve  for  l. 
However,  in  general,  o  and  x  are  functions  of  1,  providing  a  nonlinear 
transfer  equation.  In  a  uniform  medium  at  steady  state  there  can  be  no 
net  energy  accumulation  or  loss  so  the  emission  coefficient  and  opacity 
are  related  by 

n(v)  =  x(v)  I(n,v)  (E-4 ) 

which  is  Kirchoff's  law.  For  an  enclosure  in  strict  thermodynamic  equili¬ 
brium  the  radiation  field  is  given  by  the  Planck  function  at  temperature  T 
so 

n(v)  =  x(v)  Bv(T)  .  (E-5) 

If  a  medium  has  only  small  gradient  lengths  relative  to  photon 
mean  free  paths  the  approximation  of  local  thermodynamic  equilibrium  (LTE) 
is  often  useful  and  this  implies 

n(z,v)  =  x(z,v)  BV(T)  (E-6) 

so  that  the  emissivity  can  be  expressed  in  terms  of  the  opacity  and  the 
known  Planck  function.  Here  B^T)  has  the  property  that  the  flux  from  a 
blackbody  is 


0 


which  defines  the  Stefan-Boltzmann  constant,  o. 


Consider  now  an  idealized  atmospheric  layer  which  has  the  pro¬ 
perties  that  there  are  a  finite  number  of  infrared  frequency  ranges  of 
importance  and  each  of  these  has  constant  opacity.  We  assume  that  all 
thermal  radiation  takes  place  in  these  groups  and  that  the  solar  spectrum 
does  not  significantly  overlap  the  IR.  Thus  there  is  no  emission  into  the 
solar  spectrum.  The  layer  then  serves  to  absorb  the  solar  radiation  and 
convert  the  energy  to  thermal . 


THREE  GROUPS 


In  order  to  get  an  analytic  handle  on  the  problem,  we  will  now 
develop  a  set  of  transfer  equations  appropriate  for  a  three-group 
two-stream  approximation  where  the  three  groups  consist  of  two  IR  groups 
and  one  solar  group.  The  transfer  equation  under  LTE  is 


u  =  x(v)(B  (T)  -  I) 

o  Z  v 


Integrating  over  the  frequency  domains  gives 


(E-8) 


M  II) 


(E-9a) 


Here  the  are  the  opacities  of  IR  groups  and  o.  is  the  fraction  of  the 
thermal  emission  which  takes  place  into  the  ith  group.  Both  and  ai 

will  be  taken  as  known.  (The  opacities  x-  may  be  considered  to  be  the 
product  of  the  density  and  mass  absorption  coefficient.) 

TWO  STREAMS 

The  full  solution  to  the  transfer  equations  includes  the  angular 
dependence  of  the  specific  intensity.  The  details  of  this  dependence  are 
not  expected  to  have  a  significant  impact  on  the  result  and  it  is  conve¬ 
nient  to  make  an  approximation  which  allows  closure  of  the  set  of  moment 
equations  generated  by  Equation  E-3.  Eddington's  approximation  provides 
that  the  upward  and  downward  components  of  I  are  each  independent  of  u  but 
different  from  each  other. 

The  flux  associated  with  I  is 


Flux  =/  ul(u)  dp  d4>  (E-10) 

and  defining  I+  i  I(u  >  0)  and  1“  =  I(y  <  0)  the  upward  flux  is 

l 

U  =  /  u  I+dy  /  d*  =  ttI+  (E-ll) 

o 

and  the  downward  flux  is 

o 

D  =  -  /  u  1“  dp  /  d<t>  =  Tip  .  (E-12) 

-l 

Thus,  from  Equation  E-7  upon  integration  over  angle 

~  [U.J  =  2  xppT4  -  2  xiUi  ( E - 1 3 a ) 

l-Di)  =  2  aT4  -2XjD.  ( E -1 3 b ) 

a  7  1  i  i 
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for  i  =  1,2.  Since  the  solar  flux  is  only  down,  irl$  =  S  and 


liSl  *  2  v  s 

3Z 


(E-14) 


Finally  taking  x  =  -2z  as  a  convenient  unit  of  depth  and  expressing  the  1R 
fluxes  in  terms  of  plus,  P  =  U+D,  and  minus,  M=U-D  components,  we  arrive 


(E-15a) 


( E  - 1 5b  > 


(S)  = 
3x 


(E-15c ) 


which  give  the  three-group,  two-stream  transfer  equations.  The  steady 
state  condition  implies 


(y.M.-  S)  =  0  . 

3x  1 


(E-16 ) 


These  results  are  correct  for  any  number  of  IR  groups.  However,  the 
essential  features  appear  with  two  groups  and  this  case  can  be  solved 
analytically,  providing  a  solution  as  follows: 


S  =  S0  exp(-x  t) 


(E-17a ) 


M!  =  A  exp(/x  t)  +  B  exp(/y  x)  +  X  exp(/7  x) 


=  S  -  M, 


(E-17b) 

(E-17c) 


pi  *  Pio  +  x i (exp  (/y  x)  -1)  -  -8~  (exp (/y  x)  -1) 

/y  /y 


(E-17d) 


+  (exp  (“x^  x)  -  1)] 
*s 


-V.*-  A.NV.V. 


P2  »  P20  +  x i [ “  (exp (/y  t)  -1) 

/V 

-  S-0-2*  (exp(-x  t )  -  1)] 

*s 


+ 


-B  (exp(/y  t)  -1) 
/y 


(E-17e) 


where  S0  is  the  incident  solar  flux  at  the  upper  surface,  t=0. 


=  C 01  i  x2  +  a2xi)xix2 

(°lXl  +  CI2X2) 


and 


2  .  2 


X  =  s, 


“1X1  x2  -X 


«ixi  +  «2X2  y  -  X21 


(E-18) 


( E - 19 ) 


Finite  fluxes  in  an  infinitely  deep  layer  requires  A=0  so  that 


—  f (S0-X)(l  *  5.)  ♦  2  D20l  -  L  [*(1*  *i)  *  2  Dl0) 

B  =  52. _ X2 _  al _ Xi _ 

/y  r  — l —  +  — 1 — 1  +  I —  +  .I — 

v  CJI  x  1  <*2X2  a!  a2 

(E-20) 

where  DiQ  provides  for  possible  non-zero  downward  fluxes  at  t=0.  The 
boundary  conditions  at  t=0  are 

Mio  *  B  ♦  X  =  P10  -  2  D  10  (E-21a) 

and 

^20  =  Sq  "  ^10  =  P20  ”  ^  D2o  •  ( E -2 lb ) 

These  conditions  provide  the  full  solution  for  an  infinitely  deep  region 
consisting  of  two  IR  groups  and  a  solar  group. 


To  illustrate  the  behavior  exhibited  by  these  conditions  for  a 
single  layer,  consider  an  example  of  equal  strength  IR  groups  («.=  a2= 


'wimmrwwimrmrwB  wmmsm  wwi 


0.5)  with  arbitrary  normalization  of  xi  =  1  and  x2  =  10"2.  The  tempera¬ 
ture  profiles  which  result  are  shown  for  several  values  of  xs  in  Figure 
E-2.  If  xs  <  X2  the  temperature  increases  with  depth  to  an  asymptotic 
value.  If  X2  <  x$<  Xi  the  temperature  reaches  a  maximum  at  a  depth  com¬ 
parable  to  x”^»  or  at  a  solar  optical  depth  near  1.  If  x$>  Xi  the  tempera 
ture  monotonical ly  decreases  while  if  x  »  Xi  one  sees  structure  at  both  of 
the  depth  scales  of  x^1  and  (/x’iX^)"1  ^  but  the  temperature  still 

decreases  monotonical ly. 


As  illustrated  by  the  previous  discussion  of  soot  layer  optical 
properties,  the  activities  of  the  equivalent  IR  groups  span  that  for  the 
solar  radiation.  In  order  to  demonstrate  the  basic  features  expected  in  a 
realistic  temperature  profile,  let  us  approximate  the  water-vapor  emis- 
sivity  using  just  two  IR  groups.  We  will  take 


e(z)  -  .4  (1  -  e'12>5z)  +  .3(1  -  e'*0625z)  (E-22) 

where  z  is  in  kilometers.  Thus  we  are  using  aj  =  .4,  Xi  =  12.5  km-1,  a2  = 
.3,  and  X2  =  t/16  km*1,  while  taking  x$=  1/4  km-1.  The  above  fit  to  e"(U) 
is  not  very  good  unless  U*  1  but  the  essential  feature  is  preserved-- 
namely,  there  are  IR  groups  both  more  and  less  active  than  the  solar.  The 
exact  values  for  a1  and  xi  are  not  well  determined  by  such  a  fit  but  it 
turns  out  that  the  essential  features  of  solution  remain  largely  indepen¬ 
dent  of  the  selection  of  parameter  values. 


The  temperature  profile  which  results  from  the  two  IR  groups 
plus  one  solar  group  approximation  for  two  streams  is  shown  in  Figure 
E-3.  Since  x2  <  x ,/  Xi,  there  is  a  maximum  temperature  which  occurs  and 
it  is  at  a  depth  of  t  =  2z  =  2.2  km,  which  corresponds  to  a  solar  optical 
depth  of  *1/2.  The  temperature  gradient  at  the  top  of  the  layer  greatly 
exceeds  the  critical  lapse  rate  of  about  7  degrees  per  kilometer  so  that 
the  profile  will  be  unstable  against  convection.  The  resulting  convection 


will  alter  the  profile  so  that  the  gradient  will  not  exceed  the  critical 
lapse  rate  anywhere.  It  is  not  clear  that  there  will  exist  a  steady  state 
profile  maintained  by  convection.  Without  entering  into  a  dynamical 
calculation  we  now  want  to  consider  the  character  of  a  possible  steady 
state  profile  which  is  maintained  by  convection,  allowing  for  the 
possibility  that  a  portion  of  the  layer  may  mix  with  air  above  to  reduce 
the  opacities  and  effectively  expand  the  layer. 

CONVECTIVE  LIMIT 

In  a  region  in  which  convection  enforces  a  maximum  gradient 
equal  to  some  critical  lapse  rate,  8 ,  our  previous  solution  for  2  IR 
groups  is  inappropriate.  In  such  a  region  the  gradient  will  be  constant 
at  S  and  there  will  exist  an  additional  flux  from  the  convection. 

The  transfer  equations  for  the  IR  and  solar  groups  are  the  same 
as  in  the  absence  of  convection.  Namely, 


fr(pi)"xiMi  <E-23a> 

^  I.M.)  =  XjPj  -  2  x(  (E-23b) 

-X  T 

s  =  So  e  (E-23c) 

but  the  steady  state  condition  becomes 

d—  QMr  S  +  F)  =  0  (E-24) 


where  F  is  the  convective  flux, 
restricted  to  the  form 


Furthermore,  the  temperature  now  is 


where  8  is  the  critical  gradient  and  T0  is  an  unknown  temperature  at  the 
top  of  the  layer. 


Eliminating  using  Equations  E-23a  and  b,  while  making  use  of 
Equation  E-25,  yields 


32P, 


L  -r  *  *iPi  *  2Xl“iTo4(l  *  (f-  )t) 

Xi  d-c  'o 


(E-26) 


which  has  a  general  solution  of  the  form 


P,  =  A.e 


-x,t 


xiT  4  i 

+  Bie  +  l  aijT 
j=0 


(E-27) 


The  coefficients  a-,  can  be  expressed  in  terms  of  T0  as 

*  si 


a)4  ■  2  «,«  V  (e/T„)“ 

a,3  -  2  V  To'^CB/T,,)3) 

at2  ■  2  v  T(,4(6(S/T0)2)  *  12  a,,/*,2  (E-28) 

a,i  ■  2  T0*(4(b/T01)  4  6  a,3/x,2 

a10  ’  2  V  T»'  *  2  a12'*12  ' 

Using  Equation  (E-23a)  the  can  be  written  as 

"xiT  xi T  4  i-l 

M  =  -A.e  1  +  B.e  1  +  I  j  a..TJ  1  .  (E-29) 
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For  no  downward  flux  at  the  boundary  at  t=0 


Pi0  =  Mi0 


Ai  *  (ail/xi  ‘  ai0^2  * 

Furthermore,  fixing  the  convective  flux  to  be  zero  at  t=0  forces 


which  requires 


l  pio  Mi0  •  S< 


l  Bi  -  S0  -  l  (ail/Xi  ♦  ai0)/2 


(E-30) 


(E-31) 


(E-32) 


(E— 33 ) 


Thus  the  remaining  unknowns  are  T0  and  Bj.  (From  Equation  E-33  B  j+  B2 
can  be  written  in  terms  of  T0.  So  it  is  only  the  split  between  the  Bi 
which  is  unknown.) 

The  convective  flux  can  be  found  by  integrating  Equation 
E-24,  giving 

F  •  S0{e'Xs’  -  1)  -  (J  A.(e'XiT  -1).[  B1(eX,T  -  1) 


♦  5  \  J(^) 

i  j=2  Xi 


( E  -  34  ) 


where  the  condition  F=0  at  t=0  is  imposed  giving  no  convective  flux  at  the 
boundary. 

Figure  E-3  indicates  that  in  the  absence  of  convection  there 
will  be  a  finite  region  with  a  temperature  gradient  exceeding  the  critical 
lapse  rate.  Since  this  region  is  unstable  to  convection,  it  will  be 


altered  through  convection,  perhaps  including  a  mixing  with  the  inactive 
region  above  and  also  perhaps  dynamically  changing  the  previously  stable 
region  below.  Without  exploring  the  dynamics  of  such  processes,  we  shall 
now  consider  the  possibility  that  there  may  exist  a  different  steady  state 
equilibrium  in  which  convection  holds  the  upper  region  at  a  gradient  3 
while  the  semi-infinite  lower  region  is  unaltered  except  for  possible 
changes  in  boundary  conditions  at  the  interface.  This  new  configuration 
is  to  allow  the  possibility  that,  through  mixing,  the  upper  region  is  uni¬ 
formly  stretched,  if  necessary,  to  provide  a  solution  which  matches  the 
boundary  conditions  at  the  interface.  This  stretching  will  be  carried  out 
by  a  uniform  reduction  of  the  x's  in  the  upper  region  such  that  their 
ratios  are  fixed. 

The  conditions  to  be  met  are  as  follows: 

The  convective  flux,  F,  is  to  vanish  at  the  top  (t=0) 
and  at  the  interface,  t0. 

The  gradient  of  the  convective  flux,  3F/3t,  is  the 
vanish  at  the  interface--this  forces  the  temperature  to 
be  continuous. 

All  IR  fluxes  M.j  and  are  to  be  continuous  at  the 
interface--actual  1y  only  one  of  the  four  conditions  is 
independent. 

The  maximum  temperature  gradient  in  the  lower  non- 
convective  region  must  be  no  greater  than  b. 

In  general,  it  is  found  that  parameters  t0.  T0,  ancl  Bi  (°r*  equivalently, 
the  mix  between  Bj  and  B2)  can  be  found  which  satisfy  the  above  conditions 
without  any  stretching  or  mixing  of  the  upper  convective  region. 
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It  is  further  found  that  there  exist  solutions  with  a  stretching 
up  to  some  maximum  value  beyond  which  the  temperature  gradient  at  the  top 
of  the  nonconvecti ve  lower  region  exceeds  8.  The  maximum  allowed  stretch¬ 
ing  is  generally  no  more  than  about  50  percent  for  typical  parameters. 
Any  degree  of  stretching  less  than  this  maximum  results  in  a  discontinuity 
in  the  temperature  gradient  at  the  interface.  While  there  is  no  physical 
reason  to  rule  out  such  a  possibility,  it  seems  likely  that  the  dynamics 
of  the  change  from  an  unstable  profile  to  one  maintained  by  convection 
would  not  leave  such  a  discontinuity.  This  reasoning  is  further  supported 
by  the  example  in  Figure  E-4,  which  indicates  that  the  final  convective 
region  is  much  larger  than  the  original  region  for  which  dT/dT  >  B* 

Figure  E-4  shows  three  temperature  profiles  which  result  from 
opacities  which  approximate  those  for  a  soot  layer.  The  profiles  are  for 
a  fully  nonconvecti ve  equilibrium  (ORIGINAL),  for  an  unstretched  con¬ 
vective  layer  above  a  stable  layer  (UNSTRETCHED),  and  for  a  stretched 
convective  layer  above  a  stable  layer  such  that  the  temperature  gradient 
is  continuous  across  the  interface  (STRETCHED). 

The  effect  of  introduction  of  a  convectively  stabilized  layer  is 
generally  a  reduction  in  temperatures  at  the  peak  and  at  the  bottom  of  the 
region  with  an  attendant  sharp  increase  in  the  temperature  at  the  top  of 
the  layer.  The  temperature  peak  of  the  original  region  is  at  a  solar 
optical  depth,  x$t,  of  about  1/2  while  the  convection  induced  peak  is  at 
an  optical  depth  of  about  2. 

REMARKS 

The  goal  of  the,  currently  incomplete,  work  described  in  this 
appendix  is  to  understand  the  development  of  the  instability  at  the  top  of 
the  soot  layer  when  no  condensation  cloud  is  present.  The  most  obvious 
conclusion  of  the  work  so  far  is  that  it  is  necessary  to  resolve  the  IR 
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opacity  of  H20  into  more  than  one  band  unlike  was  done  in  Chapter  3.  If 
there  are  portions  of  the  IR  band  where  the  opacity  is  higher  than  is  the 

solar  opacity  and  portions  where  it  is  lower  than  the  solar  one  obtains 

behavior--the  details  of  which  cannot  be  simulated  by  any  single  IR  band. 
These  details  are  not  of  great  importance  for  estimating  the  maximum 
temperature  or  the  temperature  level  in  the  cloud,  such  as  at  the  earth's 
surface,  but  they  are  very  important  for  understanding  stability  proper¬ 
ties.  The  opacity  of  water-vapor  is  so  high  in  the  ten  to  twenty  percent 

of  the  IR  band  which  has  the  highest  opacity  and  so  low  in  the  thirty 

percent  which  is  lower  that  we  believe  there  will  always  be  substantial 
parts  of  the  IR  band  when  the  opacity  is  on  either  side  of  the  solar.  For 
these  reasons  amounts  of  water-vapor  which  are  too  small  to  affect  the 
main  character  of  the  temperature  profile  as  calculated  in  Chapter  3  can 
be  important  to  the  development  of  the  instability.  Note  that  the  IR 
opacity  of  condensed  water  is  quite  constant  as  a  function  of  frequency  so 
the  problem  under  discussion  here  must  affect  the  stability  calculations 
of  Chapter  3  which  were  for  condensed  water. 
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